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Max von Laue 
9. 10. 1879 bis 24. 4. 1960 


Mitten aus intensiver Tatigkeit wurde uns Laue entrissen: Auf der Fahrt 
zu wissenschaftlichen Besprechungen erlitt er den Unfall, dem er zwei Wochen 
später erlag. 


Im Fritz-Haber-Institut war ihm, auch als er die Leitung desselben auf- 
gegeben hatte, noch jede Arbeitsmöglichkeit gegeben, wovon er regen Ge- 
brauch machte. Wer Laue am 9. Oktober vorigen Jahres bei der Feier seines 
80. Geburtstages in voller geistiger Frische sprechen hörte, zweifelte nicht 
daran, daß er sein Lebenswerk, entlastet von Verwaltungsarbeit, noch weiter- 
führen werde. Nun ist diese Zuversicht plötzlich zerstört! Für diejenigen, 
die mit ihm im mündlichen oder schriftlichen wissenschaftlichem Austausch 
standen, ist ein schwerer persönlicher wissenschaftlicher Verlust entstanden. 
Aber Laues Lebensarbeit bleibt weiter wirksam: Obwohl seit der Laueschen 
Idee, durch Beugung der Röntgenstrahlen am Kristallgitter Interferenzen zu 
erzeugen, bald 50 Jahre verflossen sind, ist die Anwendung der Röntgen- 
spektroskopie in Physik, Chemie und Technik noch immer im Steigen. Und 
auch viele andere Arbeiten Laues sind noch lebendig. Wie sehr seine wissen- 
schaftlichen Leistungen bewundert werden, kam wieder bei der Feier von 
Laues 80. Geburtstag in Ansprachen von Vertretern der Physik, der Chemie, 
der Mineralogie und der Technik zum Ausdruck. 


Aber ebenso unvergessen wie als Wissenschaftler wird Laue als Mensch 
sein: Seine menschliche Größe zeigte sich besonders in der schweren national- 
sozialistischen Zeit, in der er allen Anstürmen gegenüber standhaft blieb, 
ja sogar vielen half, dem Unheil zu entgehen und trotz großer Schwierig- 
keiten Verbindung mit ihnen aufrecht hielt. In den letzten Jahren gehörte er 
wieder zu den „Göttingern‘“, die mit Otto Hahn dagegen auftraten, daß die 
Kernenergie statt nur für nützliche Zwecke auch zur Zerstörung bereitgestellt 
wird. Durch seine hohen seelischen Eigenschaften erwarb sich Laue nicht nur 
bei seinen Kollegen, sondern überall eine solche Achtung, daß er auch im Aus- 
land, wie sein Duzfreund Einstein, als einer der ganz Großen angesehen 
wird. Diese umfassende Anerkennung fand auch bei der Feier seines 80. Geburts- 
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tages lebendigen Ausdruck, indem zu den vielen früheren Ehrungen wieder neue 
dazukamen: Das große Verdienstkreuz, nun mit Stern und Schulterband, die 
Ernst-Reuter-Medaille des Senats von Berlin, die Helmholtz-Medaille, die 
höchste Auszeichnung der Deutschen Akademie, der Laue seit langem ange- 
hörte. 200 Kristallographen aus USA hatten ein Glückwunschschreiben der 
American Crystallographic Association unterschrieben. 


Die Annalen der Physik haben noch besonderen Grund, um Laue zu trau- 
ern. War er doch ein besonders aktives Mitglied ihres Kuratoriums. Laues 
engere Freunde aber, von denen ja nur noch ganz wenige leben, verlieren in 
ihm wie seine Gattin, die ihm immer, auch auf seinen Dienstreisen, treu um- 
sorgte, wie sein Sohn und seine Tochter sowie deren beide Familien, einen 
verehrten Menschen, an den sie immer mit bewunderndem Stolz und inniger 
Freude zurückdenken werden. 


W. Meißner, München. 


On 
St 
condt 
The < 
as we 
that 
condı 
mitti 
dispe 
the ı 
peak 
flecti 
T 
depo 
Son 
depe 
quen 
(1)8 
tion 
relay 
7 
and 
of tl 
This 
For 
integ 
disp 
(1953 


On the Electrical Conductivity of Polar Semiconductors 
at High Frequencies 


By P. H. Fang 
With 2 Figures 


Abstract 


Stolz’ expression for the complex electrical conductivity of polar semi- 
conductors as a function of the applied frequency is numerically evaluated. 
The conductivity is represented in the usual form of the frequency spectrum, 
as well as in the form of an Argand diagram. The latter representation shows 
that this dispersion function of the complex conductivity of polar semi- 
conductors is almost identical to the dispersion function of the dielectric per- 
mittivity given by Cole and Cole. From the approximate symmetry of the 
dispersion function in logarithmic relaxationtime space, it is concluded that 
the relaxation time can be determined with reasonable accuracy from the 
peak of the imaginary part of the complex conductivity, or from the in- 
flection point of the real part. 


I. Introduction 


The apparent electrical conductivity of semiconductors is influenced by 
depolarizing fields. This influence has been discussed by Howarth and 
» Sondheimer?) and Stolz?). Howarth and Sondheimer gave the energy 
dependence of the relaxation time. Stolz derived the formula for the fre- 
quency dependence of the complex conductivity. In the present paper: 
(1) Stolz’ formula will be calculated numerically; (2) a method of representa- 
tion of numerical values of the conductivity will be introduced; and, (3) the 
relaxation time distribution will be discussed. 


II. Complex conductivity 


The dispersion functions for the complex conductivity (ce. e.) of non-polar 
and polar semiconductors were derived by Stolz?) under the usual assumption 
of the existence of a relaxation time, t, for the electron-lattice interaction. 
3 This relaxation time is dependent on the energy distribution of the electrons. 
For degenerated non-polar semiconductors, the evaluation of the conductivity 
# integral is quite simple, and the dispersion for the c. c. leads to a resonance 
dispersion similar to that of an harmonic oscillator. 


ass J. Howarth and E.A. Sondheimer, Proc. Roy. Soc. London A 219, 53 


*) H. Stolz, Ann. Physik 19, 394 (1957); 1, 334 (1958). 
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In the case of polar semiconductors, Stolz derived the formula for the 
c. c. by using the result of Howarth and Sondheimer for the relaxation 
time T, i.e. 
@M(ha,)* 
V2m* 
where E and m* the kinetic energy and the effective mass of electrons respec. 
tively, a the lattice constant and ®, the limiting frequency of the longitudinal 
optical lattice oscillations. M is the reduced mass of ions of the polar crystal 
An abbreviation symbol X, will be introduced such that Eq. (1) is represented 

Ya 

as T= K, . 

The ce. c. of the electron is, according to Stolz, given by: 


t(E) E'l: (1) 


E 4e? 
where e= —, K, = -—-_n,, with n, the number density of the electrons. 
kT m* 

Substituting (1) into (2), we obtain, 

K 
= K de. 3 
” 

0 


The quantity @ is connected with the external frequency, w, through the 
relation, 


=0(1-%) (4) 
where w, is the frequency parameter associated with the depolarizing field’). 
Stolz evaluated Eq. (3) in terms of polynomials which involve the exponential 
integral and the error integral. For low frequencies, simple asymptotic expres- 
sions were developed. We have also developed asymptotic expressions appli- 
cable at high frequencies. 

Moreover, since the relaxation frequency, from which important characte. 
ristics of the relaxation process are to be found, cannot be determined from 
the asymptotic forms, it is necessary to evaluate Eq. (3) directly. In the 
following, we give the procedure for numerical evaluation. 

In the region near the relaxation frequency, Eq. (3) can be evaluated by 
using the numerical table of Dingle et al?). In their notation: 


1 
(x) = p!, (3) 
Therefore, (3) becomes, after the separation of the real and imaginary parts, 
o, = K, Ky [2a — i(2.5!) x's 5 (zx)] (6) 


with x= (@ K,)-*. Numerical values of W, for the above values of indes 
with x from x = 0 to 1 at 0.1 interval and 1 to 20 at 1.0 interval can be found 


3) R.B. Dingle, D. Arndt and S.K. Roy, Appl. Sci. Res. B6, 1441 (1957). 
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in Dingle’s table. This range of x covers the relaxation frequency well. To 
extend the range, for x > 20, Stolz’ asymptotic form can be used. For 
x < 0.1, asymptotic forms are given by Dingle which involve some tran- 
scendental functions. But in the range of x, values of o, accurate to the third 
decimal place can be obtained by using following form: 


1 x x? 
x 
p-l1 
for p= 2. Both expressions for A, are simplified from those of Dingle. 


for p = 2.5 and 


=—[1- (7b) 


III. Representation 


In order to analyze experimental data and to compare them with theore- 
tical values, graphical representations are sometimes helpful. The usual way 
of representing the complex con- 
ductivity data is to plot the fre- 
quency spectra of and o;’. 
This is shown in Fig. 1. 

In the literature on dielectric 
dispersion’) and paramagnetic 
dispersion®), the Argand dia- 
gram of the complex quantity is 
sometimes used. The Argand 
diagram is the frequency trace 
in the space of the complex . 


quantities. The use ofthismethod 106 (X) 


has the Fig. 1. The frequency spectrum of the complex 
one to determine the characte- conductivity of polar semiconductors from 
ristic relaxation time and also the Eq. (6) 


nature of the distribution of the 
relaxation times. In Fig. 2 we have shown the Argand diagram of the com- 
plex conductivity given by Eq. (6). In this way, we find that the are in the 
Argand diagram is almost 


identical with at most a St 
deviation of + 0.5%, to the Ri 

curve given by the disper- 1 
sion function of Cole and "st a 
Cole‘), with (1—a)= 3°. af 
The parameter x is a mea- Fr R 
sure of the depression of the ol 

center of a semicircle from 3 


4K.S. Cole and R.H. 
Cole, J. Chem. Phys. 19, 1484 
(1941). 

5) H. B. C. Casimir and 
F.K. DuPré, Physica 5, 507 Fig. 2. The Argand diagram of the complex conduc- 
(1938). tivity of polar semiconductors from Eq. (6) 
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the real axis. It is interesting to observe that in the dispersion of the 
conductivity of plasmas to microwaves, for the particular case in which the 
collision cross-section is inversely proportional to the velocity, the Argand 
diagram also is very similar to the arc of the Cole and Cole function, with 


the parameter 5 (1 — x) = 4°), 


IV. Distribution function of the Relaxation times 


From the fact that Fig. 2 does not exhibit a perfect semi-circle, it can be 
shown’) that there is a distribution of the relaxation times for the electrical 
conductivity of polar semiconductors. Now according to the formal definition 
of relaxation time distribution’), the complex dispersion quantity o, can be 
written as: 

0,(0) = (w) + to’ (w) 


g(t) (8) 


Here g(r) is the distribution function of relaxation times. If there is no distri- 
bution of the relaxation times g(t) becomes a Dirac delta function. 

In the present case of the conductivity of polar semiconductors, we find, 
by comparing Eq. (3) with Eq. (8): 


g(t) = 2K, (x) 9 


In the case of a non-zero distribution, i. e. there are many relaxation times 
where characteristic relaxation time r, has been defined in three different 
ways: 

i. Tg), the value of t where g(t) has a maximum; 

ii. T, the value of t such that WoT). = 1 where wo. is determined from 

(w) = 1/2 (0) — 0’(o0)]; 

iii. T)3, the value of t such that wo3 T93 = 1, where a3 is determined from 
the maximum of o” (w). 

If g(t) is symmetric in the sense that 9(=) = 6) i. e., there is 4 

01 
symmetry of g(t) with respect to a value Tp, in the log t-space, Tg; = Ty. = Typ 
and no distinction in the determination is necessary. On the other hand, 
if g(t) is not symmetrical, 1); =# Ty = T 3°). Now, 7,, has a direct connection 
with the distribution function of relaxation times, on the other hand, 7,, and 
T3 are derived quantities whose relations to the distribution of relaxation 
times depend on the specific distribution function involved. From this point 
of view, the use of 7,, is preferable for defining relaxation time characteristics. 


6) P.H. Fang, Physic. Rev. 118, 13 (1959). 

7) P.H. Fang, to be submitted to J. chem. Physic. 

8) See, for example, J.R. Macdonald and M. K. Brachman, Rev. mod. Physic. 28, 
392 (1956). 

%) L.C. Van der Marel, J. Van der Broek, and C. J. Gorter, Physica 23, 361 
(1957); P.H. Fang, Physica 24, 970 (1958). 
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By maximizing Eq. (9), 75, is obtained: 


Tor = (2.500) Ky. (10) 
To examine the symmetry property, we have, with y = log ( = ) 
01 
01 
(=) g = exp (10y — 10 Sinh y) 


(11) 
mw exp (F y) for small values of y. 


Therefore, g(t) is not symmetrical in the over-all log T-space, but it is almost 
symmetrical for small values of y, i.e. near r=r,. From graphical inter- 
polation, we find: 
= (2.570)' Ky (12) 
To3 = (2.561)'/: K,. (13) 
Since the differences are very small from the numerical point of view, one 
can ignore them. 

In conclusion, we have calculated numerically the complex conductivity of 
polar semiconductors as a function of frequency parameter (@ K,). The 
Argand diagram of this conductivity is very close to that produced by the 
Cole and Cole function with zu —«) = 3°. Finally, the distribution 
function of relaxation times of this conductivity is discussed. In this way, 
we have shown that this distribution function is almost symmetrical in the 
log t-space, and therefore an approximately accurate relaxation time can be 
determined directly from the data of the real or the imaginary part of con- 
ductivity. 
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Elektronenenergieniveaus in Kristallen. I 


Thalliumaktivierte Alkalihalogenide’) 
Von M. Balarin?) 
Mit 2 Abbildungen 


Inhaltsiibersicht 


Mit Hilfe des Zentralionenmodells wird der Einfluß verschiedener Alkali- 
halogenid-Grundkristalle auf die Elektronenterme betrachtet. Für Tl*-ak- 
tivierte Kristallphosphore wird nach Seitz angenommen, daß ein Aktivatorion 
ein normales Kation substituiert. Es wird gezeigt, daß Gesetzmäßigkeiten, 
die aus experimentellen Arbeiten seit langem bekannt sind, sich auch aus der 
Theorie ergeben, wenn die Gitterionen nicht als punktförmige Ladungen an- 
gesehen werden, sondern die räumliche Verteilung der Elektronenwolken be- 
rücksichtigt wird. Sich anbietende Schlußfolgerungen werden diskutiert und 
anderen Arbeiten gegenübergestellt. Im Anhang wird eine benutzte analytische 


Approximation für die „‚nichtpunktförmige‘‘ Ladung von Halogen- und Alkali- 
ionen angegeben. 


Bekanntlich lassen sich Berechnungen von Kristallsystemen nur nach sehr 
weitgehenden mathematischen Vereinfachungen, die nicht immer physikalisch 
ausreichend gerechtfertigt werden können, näherungsweise auswerten. Da- 
bei treten sowohl praktisch rechnerische Schwierigkeiten auf, wie z. B. der 
Umfang der sukzessiven Näherung eines ,,self-consistent‘‘ Hartree-Fock-Sy- 
stems mit an und für sich unendlich vielen Gleichungen, als auch prinzipielle 
Probleme, so z. B., in welcher Form kann die Austauschwechselwirkung ver- 
einfacht betrachtet werden, oder wie können die Korrelation oder die Eigen- 
heiten, die durch die magnetischen Quantenzahlen ausgedrückt werden und 
in bezug auf die die Hartree-Fock-Methode unempfindlich ist, Berücksichti- 
gung finden. Es ist daher für theoretische Untersuchungen solcher Vielelek- 
tronensysteme besonders wichtig, sich eine klare und dabei oftmals idealisierte 
Vorstellung über die physikalischen Vorgänge bei der zu betrachtenden Auf- 
gabe zu schaffen, wobei man die wesentlichen Einflüsse herausfinden und diese 
in möglichst übersichtlicher Form darstellen muß. 

Im Falle starker Bindung aller Elektronen, d. h. auch der Valenzelektronen, 
wie sie bei Ionenkristallen vorliegt, wahren die einzelnen Ionen im Gitterver- 


1) Diese Arbeit entstand im wesentlichen während eines Studienaufenthaltes des 
Verfassers an der Universität Leningrad. Der Verfasser möchte der Regierung der DDR 
sowie den zuständigen Stellen der Sowjetunion, die diesen Studienaufenthalt ermöglicht 
haben, seinen Dank zum Ausdruck bringen. 

2) Jetzt Zentralinstitut für Kernphysik Rossendorf b. Dresden. 
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band in weitem Maße ihre Individualität. Abweichungen von den entsprechen- 
den freien Ionen lassen sich hier mit Hilfe des Zentralionenmodells (ZIM)-$) 
bestimmen. 

Die hauptsächlichsten Voraussetzungen, die bei der Einführung des ZIM 
gemacht wurden, sowie einige Resultate für reine Ionenkristalle wurden in- 
zwischen auch durch eine experimentelle Untersuchung der Elektronenver- 
teilung in Kristallen”) als richtig bestätigt *). 


1. Der Einfluß eines Kristallfeldes auf ein Gitterion 


Hier soll das ZIM auf gestörte Ionenkristalle angewandt werden®), und 
zwar möge ein Tl+-Ion als Aktivator ein Kation des Grundkristallgitters eines 
Alkalihalogenides substituieren. Gefragt wird nach den möglichen Energie- 
zuständen für die zu diesem Ion gehörenden Elektronen. Aus Experimenten 
ist bekannt, daß bei Tl*-aktivierten Alkalihalogeniden typische Absorptions- 
und Emissionsspektren auftreten®), die einwandfrei dem Aktivator zuzu- 
schreiben sind. Sie können mit den wichtigsten optischen Übergängen (1S, — *P,, 
18, > !P,) im freien TI*-Ion verglichen werden. Wir gehen also von der von 
Seitz!) gegebenen Vorstellung über den Aufbau derartiger Kristallphos- 
phore aus, wobei wirannehmen, daß Übergänge desoptischen Elektronszwischen 
diskreten Niveaus des Aktivatorzentrums selbst möglich sind, im Unterschied 
zu dem entsprechenden Vorgang in aktivatorfreien Alkali- und Silberhalo- 

niden. 
= In der vorliegenden Arbeit wird versucht, den Einfluß verschiedener Alkali- 
halogenid-Grundkristalle zu vergleichen. Bei einer solchen Aufgabenstellung 
wird es verständlich sein, daß mehrere Vereinfachungen in Kauf genommen 
werden mußten und daß es vor allem darauf ankam, ein qualitativ richtiges 
Bild zu erhalten. Insbesondere werden wir die lokalen und Grundgitter- 
schwingungen nicht berücksichtigen, was aber für die Form der Spektren sowie 
andere Fragen, z. B. die Löschung, wichtig wäre. In einer von Kristoffel!)12) 
mit Hilfe des ZIM entwickelten Berechnungsmethode für Aktivatorzentren 
werden die Kernkonfigurationen mit berücksichtigt. Doch ist die hier gestellte 
Aufgabe auch dort nicht direkt lösbar. Eine wichtige Vereinfachung bestand 
darin, die Wechselwirkung der Elektronen verschiedener Ionen durch effek- 
tive approximierte Potentiale zu beschreiben, woraus sich erklärt daß im 
weiteren keine S-Energie auftritt!*), obwohl keine Orthogonalitätsforderung 


®) M. I. Petrashen, A.V. Iwanowau.H.Wolf, Nachr. Leningr. Univ. 10, 29 (1956). 
one M.I. Petrashen, T. L. Gutman u. M. Balarin, Nachr. Leningr. Univ. 22, 28 

). 

5) M.I. Petrashen u. T.L. Gutman, Nachr. Akad. Wiss. UdSSR, physik. Serie, 
22, 668 (1958). 

6) N.N. Kristoffel, Arb. Inst. Phys. Astr. ESSR 7, 112 (1958). 

?) G.Schoknecht, Z. Naturf. 12a, 983 (1957). 
ly nn beider Korrektur: siehe auch M. Balarin, Z. Physik 158, 

(1960). 

8) Berechnung von F-Zentren mit Hilfe des ZIM siehe I. V. Abarenkov, Nachr. 
Leningr. Univ. 22, 14 (1958). 

®) Siehe Landolt-Börnstein, Tabellen 1/4, $99ff. 

10) F. Seitz, J. Chem. Phys. 6, 150 (1938). 

11) N.N. Kristoffel, Arb. Inst. Phys. Abstr. 1ISSR 7, 85 (1958). 

12) N.N. Kristoffel, Arb. Inst. Phys. Abstr. ESSR 10, 3 (1959). 

18) P.O. Léwdin, Thesis Uppsala 1948. 
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gestellt wurden. Eine Forderung, die, wenn auch dort nur formal, so dog 
in *) enthalten ist. In vorhergehenden Arbeiten wurde gezeigt, daB man dj 
Wellengleichungen fiir die Elektronen des Zentralions durch ein Variation 
prinzip, ausgehend von der Gesamtenergie des Kristalls, erhalten kann. Wi 
wollen die Gleichungen, denen die Einelektronenfunktionen des Zentralion 
als Eigenfunktionen genügen müssen, in der Form schreiben 


[Hy + U,+AU + A] px (f) = Ex yx (F). (I 
H, = Operator, der alle Glieder enthält, die nur vom Zentralion abhängen; 


U, + AU = Operator des elektrostatischen Potentials, das durch die ib. 
rigen Ionen hervorgerufen wird, mit U, = punktförmiger Anteil, AU = nicht. 
punktförmiger Anteil; 


A = Operator der Austauschwechselwirkung mit den benachbarten Ionen. 


Da für kubische Kristalle die Winkelabhängigkeit der Eigenfunktione 
durch eine Reihe gruppentheoretischer Arbeiten gelöst ist 14-16), die auf Bethes 
klassische Arbeit zur Termaufspaltung in Kristallen!?) zurückgehen, bleibt 
lediglich übrig, die Gleichung (1) in bezug auf die Radialfunktionen zu lösen 
und die Eigenwerte zu bestimmen. Die Störungsrechnung wird hier so dureh. 
geführt, daß die einzelnen auftretenden Einwirkungen schrittweise beriick. 
sichtigt werden. Als nullte Näherung werden daher Funktionen der freien 
Ionen verwendet. In der ersten Näherung wird das elektrostatische Potential 
punktförmig angenommener Ionenladungen berücksichtigt, und in der zweiten 
Näherung wird die räumliche Verteilung der Elektronenladungen beachtet, 
das bedeutet, daß die von uns sogenannte Nichtpunktförmigkeitskorrektur‘) 
und die Austauschkorrektur auftreten 8). 

Die Eigenfunktionen der gestörten Aufgabe in der zweiten Näherung 
kann man auf dieselbe Art erhalten wie das für die erste Näherung ursprüng- 
lich angegeben wurde®). In der Störungsenergie, die den asymptotischen Ver- 
lauf der Radialfunktionen bestimmt, treten nur additiv die zwei obengenannten 
Korrekturen zusätzlich auf. Die Eigenwerte der gestörten Aufgabe sind somit j 
gleichzeitig bestimmt zu: 


Ex = EX — Ex,v, — Ex,av — Ex,a. (2) 


Es kann als bestätigt angesehen werden, daß die Hauptgröße E% auch für die 
dem Leitfähigkeitsband nahe liegenden Elektronenniveaus trotz der Störung 
der Eigenfunktionen praktisch nicht verändert wird!?). Für die Korrektur 
der ersten Näherung ist bekannt, daB Erw, =A+nr ist, mit 


nx = [ (2-z)fkdr; fx-Radialteil der Eigenfunktion des betrachteten 


a r 


a 
Zentralions. Die Madelungenergie A ist umgekehrt proportional der Gitter- 
konstanten und für das weitere ist entscheidend, daß sie für alle Elektronen 


14) Bouchardt, R. Smoluchowsky u. E. Wigner, Phys. Rev. 50, 58 (1936). 
15) F,van der Lage, H. Bethe, Phys. Rev. 71, 612 (1947). 

16) D. Bell, Rev. Mod. Phys. 26, 311 (1954). 

17) H. Bethe, Ann. d. Phys. 8, 133 (1929). 

18) M. Balarin, Z. phys. Chem. 213, 44 (1960). 

19) J. V. Kossenkowa, Diplomarbeit Univ. Leningrad 1958. 
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eines Ions gleich groß ist und mit demselben Vorzeichen auftritt. In Ex, 47 
und Ex, a gehen die benachbarten Ionen additiv ein. Diese Korrekturen sind 
stark vom Quantenzustand des betrachteten Elektrons abhängig. Es soll gleich 
gezeigt werden, daß jeweils nur die erstnächsten Nachbarionen für diese 
Korrekturen von Bedeutung sind. Größenordnungsmäßig gilt meist 


E% >Ex,v, > Exr,a = Ex,av- (3) 


Folglich gleichen die inneren Elektronen jedes Ions den entsprechenden Elek- 
tronen des freien Ions. Für das Valenzelektron jedoch wird. EX so klein, daß 
es mit Ex,v, vergleichbar wird. Für ein optisches Elektron schwerer Akti- 
vatoren wachsen die beiden letzten Korrekturen wegen seines großen Radius 
stark an. 


2. Eine grobe Abschätzung der Punktförmigkeit der Ionenladungen 


Es soll gezeigt werden, daß man in fast allen auftretenden Fällen in dem 
Nichtpunktförmigkeitsglied nur die erstnächsten Nachbarn des zu betrach- 
tenden Ions zu berücksichtigen braucht. Dazu wollen wir das Verhältnis der 
durch die einzelnen Nachbargürtel aufgebrachten Anteile abschätzen. 


Für ein Elektron des Zentralions mit dem Zentrum in r = 0 ergibt der 
nichtpunktförmige Anteil der elektrostatischen Ladung eines Ions aus einem 


Nachbargürtel im Abstand r =a Vm ein Potential, das näherungsweise dar- 
gestellt werden kann in der Form (siehe Anhang): 


N ir—aVm| _ L Alr—aVm| _ ‚-Alr+aYm| 


(4) 


2vaym r 2haVm r 
Die sphärische Symmetrie in bezug auf das Zentrum des Zentralions folgt aus 


der Zerlegung nach der nullten kubischen Funktion, der winkelabhängigen 
Eigenfunktion der vorliegenden Aufgabe. 


Aus dem Potential (4) und der bekannten Elektronen-Radialfunktion fx (r) 
kann die zugehörige Energiekorrektur mit Hilfe der Störungstheorie numerisch 
berechnet werden. Eine Abschätzung jedoch kann einfacher gehalten werden. 
Wir wollen insbegondere die Abhängigkeit von a Vm fiir folgende zwei Extrem- 
fälle prüfen: 1. las Zentralion sei ein normales Kation; 2. es sei ein Aktivator. 


Im ersten Fall können wir ausnutzen, daß praktisch für alle Alkalihalo- 
oo a 
genide [ fk(r) dr<[ fk (r) dr ist. 
a 0 


Daraus folgt, daß die Nichtpunktförmigkeitskorrektur beschränkt wird 
auf den Wert 
N aym va\m 
- fk (r)dr = I, (5) 
v 


aym aym 


wobei J von a Vm, der oberen Grenze faktisch nicht abhängt, d.h. für alle 
Gürtel gleichgesetzt werden darf. y,, ist die Besetzungszahl des Gürtels mit 


‚dem Index m im Abstand r =a|/m vom Zentralion. 
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In diesem ersten Fall ist Ex sy immer klein, besitzt aber prinzipielle Be- 
deutung wegen des Beitrages in das Abstoßpotential. Für ein bestimmtes 
System ergibt sich, daß der Einfluß des ersten Anionengürtels größer ist als der 
Einfluß aller weiteren Kationen- und Anionengürtel zusammengenommen. Für 
ein zentrales Anion spielt neben dem ersten Kationengürtel der nächstfolgende 
Anionengürtel noch eine, wenn auch geringe Rolle. 


Im zweiten Fall jedoch erstreckt sich die Verteilung des optischen Elek- 
trons des Aktivatorions bis in den Bereich der benachbarten Ionen. Hier 


wird daher entscheidend, daß das Potential (4) in a Vm ein stets sehr spitzes 
Maximum besitzt. Das veranschaulicht die Abbildung der Hilfsfunktion @, (r), 
wobei gesetzt wurde 


eralm @,(r) = © (6) 


r 


Insbesondere gilt 


za 
®,(a Vm +e) aVm—e’ (ie <aym). 


Der in (6) und (7) nicht berücksichtigte, 
von A abhängige Teil bewirkt eine weitere 
200 | Anhebung der Spitze von (4). 

Aus diesem Grund wird für die Energie- 
korrektur der Verlauf der Eigenfunktionen 
00+ im Gebiet der benachbarten Ionen maBgeb- 
lich. Wenn mit E (Cl-, a) die Nichtpunkt- 
P förmigkeitskorrektur, bedingt durch einen 

eon ae ersten Chlornachbarn, bezeichnet wird, 

usw., so setzt sich die Nichtpunktförmig- 
Abb. 1 Schematischer Verlauf von keitskorrektur für ein Ti*-Ion in KCl 

w=1,4;a=5; m=1) 
insgesamt folgendermaßen zusammen: 


Ex, av = 6E(CH, a) + 8 E(Kt+, a/2)+... (8) 


Beschränken wir uns wieder wie in (6) auf den größeren Teil des Potentials (4), 
so kann man wie folgt abschätzen: 


-v|r—a¥m| _ -v|r +aYm| 
E(<CH,K+ >aym)=,,., : fk(r) dr 
vaym 
N or | _ 


und da fx(r)/r sich in (e’, 2a — e’) gegenüber dem anderen Faktor langsam 
ändert 


> vey, Fila Vm) (e’, »). (9) 


| 
En 
qu 
da 
un 
so 
Co 
Ve 
B 
H 
Tis 
| 
ei 
Al 
dı 
zu 
ke 
de 
io 
W 
K 
se 
| 
| : 
u 


M. Balarin: Thalliumaktivierte Alkalihalogenide 125 


Mit Hilfe der Formel (9) wurden dann die folgenden Verhältnisse erhalten?) 


Für die Zustände 

6E (CT, a) | | 163 | 
Ym E (CT, a Vm) | | 2840 


6E (CH, a) | 


== 4 _ 17,5 - 178 


3. Optische Übergänge 


Bekanntlich erhält man für den Übergang eines Valenzelektrons von einem 
Energieniveau auf ein anderes die Frequenz des dabei ausgestrahlten Licht- 
quants bzw. des dabei absorbierten, aus der Differenz der Terme. Für ein 
freies TI+-Ion erhält man z. B. 6s = — wobei zu beachten ist, 
daß es mehrere sich in der magnetischen Quantenzahl sowie der Symmetrie 
unterscheidende 6 P-Zustände gibt. Nimmt man jetzt ein Tl*-Aktivator, 
so ergibt sich bei analogem Vorgehen und unter Berücksichtigung des Frank- 
Condon-Prinzips, daß die Übergangsenergie sich um folgende Größe, die wir 
Verschiebungsgröße nennen wollen, von der im freien Ion unterscheidet: 


oo 
B(H,M; Ti+) = (E$, — E$,) — (Lop — Ess) =3 f — er’Ir+al) 


(e-älr-al — dr + 6 (ney — Nos) 
+6AE, (10) 


H und M stehen hier an Stelle der die Grundanionen bzw. -kationen charakte- 
risierenden Parameter. Das Feld punktförmiger Ladungen geht, abgesehen 
von den vernachlässigbar kleinen 4, gar nicht in diese Verschiebungsgröße 
ein. Für die Austauschkorrekturglieder ist es schwierig, eine übersichtliche 
Abschätzung durchzuführen. Anscheinend kommen für sie ebenfalls nur die 
durch die ersten Nachbarn bedingten Glieder in Betracht. Diese sind in (9) 
zudem kleiner als die entsprechenden Glieder der Nichtpunktförmigkeits- 
korrektur. 

Die Größe a kann man als Aktivatorkonstante im gegebenen Kristall be- 
zeichnen. Sie ist nicht identisch mit der Gitterkonstanten, sondern stellt 
den Abstand der Gleichgewichtslage der erstnächsten Nachbarn zum Zentral- 
ion dar. Da sich bei Anregungsprozessen die Ladungsverteilung ändert, 
wird folglich auch dieser Abstand geändert. Daher fallen bei aktivierten 
Kristallphosphoren Absorptions- und Emissionsfrequenzen für ein und den- 
selben Elektronenübergang nicht zusammen, wie das für freie Ionen der Fall 
ist. 

20) Die Radialfunktionen von Tl+* unter Berücksichtigung der Gittereinwirkung und 


für die angegebenen Zustände wurden mir dankenswerterweise von M. J. Petrashen 
und T.L. Gutman zur Verfügung gestellt. 
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Im einzelnen ergeben sich aus der Kenntnis der die Halogen- und Alkali- 
ionen kennzeichnenden Parameter sowie der Kenntnis der Radialfunktionen 
des optischen Elektrons für die wichtigsten Übergänge 1S, > ®P,; 18, > 1P, 
folgende Regeln: (entsprechende experimentell bestimmte Werte sind zum 
Vergleich in Tabelle 1 angeführt ®)) 


Tabelle 1 
Verschiebungsgröße für Tl+ in Alkalihalogeniden (in eV) 
18, > ®P, 0 1 18, > *P, 
Br- | B(CI-): B(Br-): 
| | B(J-) 
55:1,39 
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1. Die Verschiebungsgröße und damit die Lage der Absorptionsmaxima 
wird durch die ersten Nachbarn des Aktivatorions bestimmt. 

2. Dabei wächst innerhalb der Halogenreihe Cl- > Br~ — J~ die Ver- 
schiebungsgröße monoton an, unabhängig von den Kationen. 

3. Für den Übergang aus dem Grundzustand in das !P,-Niveau ist die 
Verschiebungsgröße größer als für den Übergang in das 3P,-Niveau. 

4. Daum r =a ftp, langsamer abfällt als f?p,, ergibt sich eine Verringerung 
des Verhältnisses der Verschiebungsgröße für verschiedene Halogene für den 
Übergang in das höhere Niveau. 

5. Für Isostrukturphosphore ist zu erwarten, daß eine heterogene Kationen- 
zusammensetzung keinen Einfluß hat, eine heterogene Anionenzusammenset- 
zung aber zu einer Aufspaltung und proportionalen Verschiebung der Maxima 
führt. 

6. Da nach dem Frank-Condon-Prinzip die Gleichgewichtslage der an- 
geregten Niveaus für die Emission verantwortlich ist, diese aber stärker vom 
zweiten Nachbarionengürtel, d. h. dem ersten Kationengürtel abhängen, ist 
für die Lage der Emissionsmaxima eine gewisse Abhängigkeit von den Kat- 
ionen zu erwarten. Außerdem wächst die Verschiebungsgröße für die Emission 
an.r 


4. Diskussion 


Da gerade die Alkalihalogenidphosphore schon seit über zwei Jahrzehnten 
experimentell sehr intensiv untersucht wurden und auch seit langem prak- 
tische Bedeutung besitzen, ist es verwunderlich, daß die hier angeführte Ver- 
schiebungsgröße, die so eindeutig von dem Grundkristall abhängt, bisher 
praktisch keine Beachtung gefunden hat. Es wurden zwar empirisch alle 
Feinheiten der Banden, Konzentrations- und Temperaturabhängigkeiten u.a.m. 
aufgenommen, eine systematische Gegenüberstellung der einzelnen Ergebnisse 
aber fehlte. Ein Zeichen dieses rein empirischen Vorgehens ist das Auftragen 
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von Kurven gegenüber Wellenlängeneinheiten. Man interessierte sich nur 
für die optischen Eigenschaften eines Kristalls, ohne dahinter wesentliche 
energetische Zusammenhänge bestimmter Bausteine gefunden zu haben. Die 
meisten Deutungsversuche liefen im Prinzip auf eine Beschreibung durch 
klassische Wechselwirkung hinaus. Da aber noch nicht einmal eindeutig ge- 
klärt war, worauf das Auftreten der Spektren, ihre Lage und Form beruht, 
konnten solche Fragen wie der Einfluß verschiedener Grundkristalle nicht 
gelöst werden. Auch die weitsichtigsten Erklärungen, z. B. von Schali- 
mowa”!) und Klement”) konnten mathematisch nicht fundiert werden. 
Man vergleiche nur moderne Bücher über Festkörperphysik ®), die alle ein 
Kapitel über Lumineszenz und aktivierte Kristallphosphore enthalten, meist 
sogar über Tl*-aktivierte Alkalihalogenide. 


Die Erfolge der Halbleiterphysik andererseits führten dazu, auch die 
Alkalihalogenidphosphore durch Vorstellungen und Begriffe, die der Halb- 
leiterphysik entlehnt waren, zu erfassen, wodurch die Verhältnisse noch mehr 
verwirrt 


In der Theorie wurden seit Seitz’ grundlegender Arbeit!) Berechnungs- 
methoden für Aktivatorzentren entwickelt. Bedingt durch den Umfang 
der Rechenarbeiten wurden bisher nur einzelne Systeme berechnet, und zwar 
mehrfach KCl (T1)12)2)%), sowie NaCl (Ag)?”) und KCl (In)®). 

Verständlicherweise sind alle bekannten Methoden bei genauer numerischer 
Durchrechnung so kompliziert, daß der hier so elementar erscheinende Ver- 
gleich von ihnen nicht geliefert werden kann. Die gute Übereinstimmung un- 
serer qualitativen Ergebnisse mit den experimentell bekannten Daten deutet 
aber darauf hin, daß auch nach den gemachten Idealisierungen die Aufgabe 
an und für sich noch ziemlich gut erfaßt wurde. Ist das aber der Fall, so er- 
kennt man auch hier einige Unzulänglichkeiten der Methode von Williams”), 
gegen die teilweise schon früher Bedenken erhoben wurden. Es kommt auf 
den spezifischen Verlauf des Störpotentials an. Das Abstoßungspotential 
der Born-Mayerschen Theorie ist zu grob in bezug auf die Störung des 
optischen Elektrons, vor allem für r >a. Die darin eingehenden Konstanten 
müßten dem System — ein Aktivatorion in Grundkristallumgebung — an- 
gepaßt werden. Dafür gibt es jedoch kein eindeutiges Kriterium. Noch pro- 
blematischer ist das für den angeregten Zustand. Als Wellenfunktionen des 
Aktivators dürfen nicht die des freien Ions verwendet werden und außerdem 
benötigt man Funktionen für die einzelnen möglichen angeregten Zustände, 
was bedeutet, daß die Hartree-Fock-Methode durch die halbempirische Me- 
thode ergänzt werden muß. Es liegt daher nahe, daß die verhältnismäßig gute 


2) K.V.Schalimowa, Ber. Akad. Wiss. UdSSR 97, 437 (1954). 

2) F.D. Klement, Arb. Inst. Phys. Astr. ESSR 1, 15 (1955). 

2) F. Seitz, Modern Theory of Solids, N. Y. 1940, 151; J. Kittel, Introduction to 
Solid State Physics, N. Y., London 1956, Ch. 18; A. J. Dekker, Solid State Physics, 
1958, Ch. 16; N. F. Mott, R. W. Gurney, Electronic Processes in Ionic Crystals; Ox- 
ford 1957, Ch. 6, 4. 

*) D. I. Blochinzew, Ber. Akad. Wiss. UdSSR 2, 76 (1934). Ausführlich behandelt 
bei E. I. Adirowitsch, Einige Fragen zur Theorie der Lumineszenz der Kristalle, M., L. 
1950 (Dtsch. Berlin 1953). 

5) F. Williams, J. Opt. Soc. Am. 47, 869 (1957), dort weitere Literaturhinweise. 

*) R.S. Knox u. D. J. Dexter, Phys. Rev. 104, 1245 (1956). 

®) N. D. Potechina, Opt. Spektroskopija 2, 388 (1957). 


\ 

kali- 
1p, 
zum 
r): 
37 
40 
39 
cima 
Ver- 
die 
ung 
den 
nen- 
set- | 
‘ima 
an- 
vom 
, ist 
<at- J 
sion 
iten 
rak- 
Ver- 
sher 
alle 
1.m. 
lisse 
igen 


128 Annalen der Physik. 7. Folge. Band 6. 1960 


Übereinstimmung zwischen den Resultaten von Williams für den bisher 
einzigen berechneten Übergang in KCl (Tl):1S, > 3P, und dem experimen. 
tellen Ergebnis auf die ausgewählte Anpassung der Parameter zurückzuführen 
ist. Die Berechnung von NaCl (Ag) ergab eine schlechtere Übereinstimmung, 

Es ist also festzustellen, daß die Gesamtheit aller Gitterionen die An. 
hebung der Elektronenniveaus von Kationen und die Absenkung bei Anionen 
bestimmt, und zwar betragsmäßig fast ausschließlich durch ihren punkt. 
förmigen Ladungsanteil. Die Annäherung von Niveaus eines Kations unter. 
einander hängt jedoch nur von der Nahordnung ab. Wesentlich für die opti. 
schen Eigenschaften ist also, daß das in Betracht kommende Elektron sich 
im angeregten und im Grundzustand bis in die Nähe des Kerns eines Nach- 
barions ausbreitet, wo die benachbarte Ionenladung nicht mehr punktförmig 
ist. Bedeutende Verschiebungen der Übergangsenergien ergeben sich daher 
für Schwermetallionen als Aktivatoren. Ihre optischen Elektronen besitzen 
einen großen effektiven Radius, ihre Verteilung verläuft bereits über die Nach- 
barkerne hinaus, aber sie sind auch in den ersten angeregten Zuständen noch 
um den Aktivator lokalisiert. 


Zusammenfassung 


Während bisherige theoretische Arbeiten über thalliumaktivierte Alkali- 
halogenide sich jeweils nur auf einen bestimmten Grundkristall bezogen, 
wird eine Abschätzung vorgeschlagen, die den Unterschied zwischen den ver- 
schiedenen Alkalihalogenidgrundkristallen aufzeigt. Dieser Unterschied kann 


mathematisch formuliert werden, wobei als Parameter Größen auftreten, 
die auf rein theoretischem Wege bestimmt wurden. Gegenüber Born -Mayer- 
Potentialen ist man nicht auf empirische Parameter angewiesen, die gar nicht 
der gestellten Aufgabe zu entsprechen brauchen, sondern aus andersgearteten 
Eigenschaften bestimmt werden. 

Es erweist sich als günstig, daß man nicht absolute Größen berechnen 
muß, z. B. sämtliche Eigenwerte, sondern, daß man für jeden Übergang nurf 
die Verschiebungsgrößen zu vergleichen braucht, d.h., daß alles auf ent 
sprechende Größen des Aktivatorions im freien Zustand bezogen wird. In 
den Verschiebungsgrößen wird das Kristallfeld, das in der Umgebung des Ak- 
tivators herrscht, berücksichtigt. Die sich daraus ergebenden, in Abschnitt 
angeführten Regelmäßigkeiten für die Abhängigkeit von den Grundkristallea 
befinden sich in guter qualitativer Übereinstimmung mit den experimentell 
ermittelten Daten, die bisher keine einwandfreie theoretische Deutung er 
fahren hatten. Dieses Ergebnis kann gleichzeitig als ein weiterer Beweis für 
die Gleichartigkeit der Natur der Phosphoreszenzzentren in aktivierten Alkali 
halogeniden dienen. 

Einige der angegebenen Regeln bedürfen jedoch noch der ausreichenden 
experimentellen Bestätigung. So kann aus dem experimentellen Befund fü 
Phosphore auf Mischkristallbasis nur indirekt auf eine Superposition mehrerer 
Banden geschlossen werden und für die Emission sind Angaben in der für 
einen Vergleich notwendigen Vollständigkeit nicht bekannt. 

Auch quantitativ erhält man ein zufriedenstellendes Resultat. So ergibt 
z.B. allein die Nichtpunktförmigkeitskorrektur mit dem Potential (4) fü 
KCI (Tl) und KJ (Tl) und die Übergänge 18, > 3P, und 18, > !P, rund dre 
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Viertel der experimentell ermittelten Verschiebungsgrößen, was sich mit der 
zu erwartenden Größenordnung für die übrigen Korrekturglieder und für die 
gemachten Vernachlässigungen verträgt. 


Für die Aufgabenstellung und aufmerksame Anleitung danke ich M. I. 
Petrashen. 
Anhang 


Analytische Näherung des effektiven Potentials einiger Ionen 


Gewöhnlich wird ein Ionenkristall durch die Art der vorkommenden An- 
ionen und Kationen bezeichnet und eine besonders wichtige Rolle spielt die 
Gitterkonstante. Ein punktförmiges Ionengitter wäre instabil und so wird 
üblicherweise ein Abstoßungspotential vorgegeben, dessen Parameter so be- 
stimmt werden, daß beobachtbare Größen, vor allem mechanische, sich aus 
der Theorie wieder ergeben. Für bestimmte Erscheinungen, wie z. B. optische 
Übergänge, scheinen solche phänomenologischen Potentiale nicht in allen 
Fällen ausreichend und genügend empfindlich zu sein. Andererseits ist klar, 
daß die vorhandene gegenseitige Abstoßung benachbarter Ionen durch die 
Wechselwirkung der räumlich verteilten Elektronenwolken hervorgerufen 
wird. Es lag deshalb nahe, die bisherigen Berechnungen von Ionen genauer 
auszuwerten. Insbesondere schien das effektive elektrostatische Potential 
für den Vergleich verschiedener Ionen geeignet. Dieses gibt nicht nur an, 
daß es sich um ein positives oder ein negatives Ion handelt, eine Bezeichnung 
übrigens, die erst für große Radien eindeutig ist, sondern daß im Zentrum eine 
quasipunktförmige Kernladung Z vorhanden ist, und wie um diese symmetrisch 


N°-Elektronen verteilt sind. Aus diesem Potential V (r) = > kann man 


einen punktförmigen Anteil aussondern. 
r Ir |? 


>, 


| 
Z = Kernladung; N° = Elektronenzahl; 9 = Elektronendichte. 
Als „nichtpunktförmige‘‘ Ladung kann dann 
Q,(r) = N°— dr’ =Z, + (N°—Z) 
vr? 
betrachtet werden. Offensichtlich ist Q,(r) > 0 für r > oo und Q,(r) > N° 
fürr > 0. Aus der Tabelle 2 ist zu ersehen, daß es für die Alkali- und Halo- 


genionen einige typische Größen (2,(0), 9 (0)/Q,(0), die sich 


für den jeweiligen Ionisationsgrad monoton mit der Ordnungszahl ändern. 
Das erlaubt auch eine gewisse Extrapolation für die noch nicht berechneten 
Ionen Cst, Br- und J-. 

Es wurde eine analytische Näherung gesucht, die, um vergleichende Ab- 
schätzungen durchführen zu können, für alle Ionen die gleiche Form haben 
mußte und die außerdem für die Anwendbarkeit im Rahmen des ZIM nach 
kubischen Funktionen zerlegbar sein mußte. Als brauchbar erwies sich fol- 
gende Approximation 

Q,(r) = + Let 


») Für den Einzelfall eines Chlorions war eine solche Näherung bereits früher von 
M.I. Petrashen und T.I. Gutman versucht worden. 
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Q,(0) =N+L, —Q,(0)=Nv+ Li. 


Selbstverständlich kann das Feld, das schon wegen der Struktur der inneren 
Elektronenwolken keinen solch glatten Verlauf hat, durch eine solche Nähe. 


— 125.67 5567" 
-14r 
55.673 r 


ee. 2 3r 
Abb. 2 Approximation der „nichtpunktför- 
migen‘‘ Ladung für Cl-. 


= 


rung nicht genau nachgebildet 
werden. 

Die Genauigkeit der gewählten 
Darstellung ist für den Fall eine 
Chlorions aus der Abbildung 2 m 
ersehen. Die von uns bestimmte 
Approximationsparameter findet 
man in Tabelle 2. 

Bei der Bestimmung der Ap. 
proximationsparameter wurden Er. 
gebnisse der Berechnungen freier 
Ionen verwendet. Die mit Hilfe 
des ZIM berechenbare Aufweitung 
der Kationen und Kompression 


der Anionen im Kristallfeld stören offenbar nicht das typisch monotone 
Verhalten der Parameter fiir die Ionen innerhalb jeder der Gruppen sowie 
deren Verhältnis untereinander. Dabei tendiert eine zahlenmäßig geringe 
Änderung in der Richtung der Verringerung des Verhältnisses vanaıı: YHalogen 
und das Verhältnis N:Z wächst für Kationen und fällt für Anionen. 


Tabelle 2 


Parameter für die analytische Approximation 
[extrapolierte Werte in ( )] 


|95(0) | —@(0) | —Q5(0)/2,(0) 


N r(Q,(0)/2) 


34,6 3,46 
74,5 4,14 
5,21 


5,65 


8,93 0,233 


2,65 
2,89 1,40 


(1,30) | (17,5) 
(1,24) | (22) 


12,5 


Bei der Bestimmung der Parameter erweist sich » als der empfindlichste. 
Interessant ist, daß für verschiedene Ionen sich die » (auch die A) fast umgekehrt 
proportional den entsprechenden Kristallradien verhalten. 


29) Kristallionenradius nach Pauling, The nature of the chemical bond, p. 350. 


Rossendorf b. Dresden, Zentralinstitut für Kernphysik, Bereich Werk- 


stoffe und Festkörper. 


Bei der Redaktion eingegangen am 31. Juli 1959. 
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Untersuchungen zur neuen Tamm-Dancoff- Methode*) 
Von Hans Jiirgen Kaiser 


Mit 7 Abbildungen 


Inhaltsiibersicht 


Das unendliche Gleichungssystem zwischen den Amplituden der neuen 
Tamm-Dancoff-Methode wird abgebrochen, indem man geeignete Linear- 
kombinationen (gy-Funktionen) dieser Amplituden nullsetzt. Die Vorschrift 
blieb bisher auf Tamm -Dancoff-Amplituden beschränkt, die den Vakuum- 
zustand enthalten. Eine Verallgemeinerung der y-Funktionen erlaubt die 
Berechnung von Differenzen benachbarter Eigenwerte, wodurch sich bessere 
Näherungen und einfachere Gleichungen ergeben. Das Verfahren wird an 
quantenmechanischen Beispielen untersucht und auf ein nichtlineares Bose- 
feld sowie auf Heisenbergs Elementarteilchentheorie angewandt. Es ergibt 
sich ein unendliches, diskretes Massenspektrum m, ~ 2? [3/(2 2 + 1)]"«. 


Einleitung 


Die neue Tamm-Dancoff-Methode (NTD-Methode) hat durch Heisen- 
bergs Elementarteilchentheorie ein besonderes Interesse gewonnen. Denn 
sie stellt bis jetzt die einzige Möglichkeit dar, konkrete Aussagen aus dieser 
Theorie abzuleiten. Die NTD-Gleichungen erlauben die Behandlung feld- 
theoretischer Probleme unabhängig von einer Entwicklung nach Potenzen 
der Kopplungskonstanten. Erste Erfolge hatte die Tamm-Dancoff-Methode 
bei der Untersuchung der Wechselwirkung von Mesonen mit Nukleonen, die 
wegen der Stärke der Kopplung keine Störungstheorie zuläßt. 

Die alte Tamm-Dancoff-Methode betrachtet Matrixelemente der freien 
Erzeugungsoperatoren zwischen Eigenzuständen der freien Felder. Dabei 
treten nichtrenomierbare Vakuumdivergenzen auf. Diese werden vermieden 
in der NTD-Methode, deren Matrixelemente (NTD-Amplituden) zwischen 
den unbekannten Zuständen der gekoppelten Felder genommen werden. Die 
NTD-Gleichungen stellen ein unendliches System von Integralgleichungen 
zwischen diesen Amplituden dar, das zur praktischen Behandlung abgebrochen 
werden muß. Die hierbei benutzte Vorschrift beeinflußt das Ergebnis ent- 
scheidend. 

Das war der Ausgangspunkt für die Untersuchungen in dieser Arbeit. 
In 1. wird gezeigt, daß die von Heisenberg benutzte Abschneidevorschrift 
für das NTD-Gleichungssystem der Annahme entspricht, daß die physikalischen 
Erzeugungsoperatoren genähert mit den freien Erzeugungsoperatoren über- 


*) Auszug aus der Dissertation des Verfassers, Dresden 1959. 
g* 
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einstimmen. Die freien Erzeugungsoperatoren miissen nicht notwendig zu 
demjenigen Hamiltonoperator gehören, der nach Streichen des Wechsel. 
wirkungsterms zurückbleibt. Bei starken Störungen erhält man — wie in 2, 
gezeigt wird — eine bessere Näherung oder überhaupt erst die Möglichkeit, 
das Problem zu behandeln, wenn man in einem System von freien Operatoren 
rechnet, das denselben Eigenwert hat wie die gekoppelten Felder in der be- 
trachteten NTD-Näherung. Aus einer Analyse der Abschn eidevorschrift für 
das NTD-Gleichungssystem ergibt sich eine neue Möglichkeit, höhere Energie. 
eigenwerte aus NTD-Gleichungen erster Ordnung zu gewinnen, während bis- 
her Gleichungen n-ter Ordnung erforderlich waren, um Aussagen über den 
n-ten Zustand zu erhalten. 

In 3. wird dieses Verfahren auf die mehrzeitige Darstellung übertragen. 
In 4. wird ein von Thirring angegebenes nichtlineares Bosefeld nach der 
NTD-Methode behandelt. Die Ergebnisse werden verglichen mit der Methode 
von Schiff, der das Problem in einen Raum mit diskreter Struktur trans- 
formiert hat. Die gewonnenen Erkenntnisse werden in 5. auf die Heisen- 
bergsche Elementarteilchentheorie angewandt. Es konnte ein unendliches, 
diskretes Massenspektrum berechnet werden, das an vielen Punkten dem 
empirischen nahekommt. Dabei mußte die Identifikation des niedrigsten 
Fermioneneigenwerts mit dem Nukleon fallengelassen und statt dessen die 
Masse des u-Mesons als Ausgangspunkt genommen werden. 


1. Die NTD-Gleichungen 


Das quantentheoretische Problem, zu einem gegebenen Hamiltonope- 
rator H Eigenfunktionen |%) und Energieeigenwerte E, zu bestimmen, 
führt auf die Untersuchung der Schrödingergleichung 

H |n> = E,, |. (1,1) 
Dabei seien die Lösungen des „ungestörten‘ Problems 


H, |n> = EA |n), (1,2) 


10), 


bekannt. 
Für beliebige A gilt die Beziehung 


LH, A]| m) = E,m R|A| m), Eun =E,-E,. (1,3) 


wie man durch Anwendung von H nach links bzw. nach rechts nachweist. 
(1,3) ist die Grundgleichung der NTD-Methode. Sie erlaubt nur die Bestim- 
mung von Energiedifferenzen. Das hat zugleich den Vorteil, daß Vakuum- 
divergenzen überhaupt nicht in Erscheinung treten. Aus (1,3) läßt sich ein 
(im allgemeinen unendliches) Gleichungssystem zwischen den o-Amplituden 


Onm (k\l) = a!| m) (1,4) 


ableiten. In der erforderlichen Vorschrift zum Abbrechen dieses Gleichungs- 
systems liegt die eigentliche Problematik der NTD-Methode. Sie wird zu- 
nächst für den Fall der Quantenmechanik untersucht. Zur Anwendung auf 
die Feldtheorie sind einige Modifikationen erforderlich, die in 3. angegeben 


werden. 
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Die Tamm-Dancoff-Methode ist eine Entwicklung nach der Zahl der 
virtuellen Teilchen. Da die Teilchen der Feldtheorie hier den Anregungs- 
stufen des Oszillators entsprechen, heißt das, die Eigenfunktionen des Hamil- 
tonoperators H sollen durch Zustände |m) endlicher Anregung genügend 
genau approximiert werden können: 


En |m> nm |m> fir N > n. (1,5) 


Unmittelbar benutzt wird (1,5) in der alten Tamm-Dancoff-Methode. NTD 
macht nur die weniger einschränkende Annahme, daß niedrige Potenzen der 
freien Erzeugungsoperatoren genügen, wenn es sich um die Untersuchung 
von Übergängen zwischen tiefen Zuständen handelt: 


Ono (k |) & 0 fir k,I>n. (1,6) 


In 1. Näherung berücksichtigt man hiernach nur Zustandsänderungen 
um ein Teilchen: 
(1,7) 


Für Übergänge zwischen benachbarten Zuständen gibt (1,7) einen ersten 
Aufschluß, während andere Matrixelemente und Energiedifferenzen, z. B. E,, 
hierdurch überhaupt nicht approximiert 
werden können. 3 

Die Verhältnisse werden illustriert durch 
Abb. 1. Die horizontalen Linien sind die , 
freien Eigenzustände |n). Die dicken 
Punkte liegen auf den Niveaus |n) und |0), 
wenn Übergänge 0,9 betrachtet werden. 
im darüberliegenden Punkte geben weitere 4 
m) an, die am Aufbau der betreffenden a 2 
physikalischen Zustände beteiligt sind. 
Wenn der Hamiltonoperator symmetrisch 
ist, H(q) = H(—q), so können das nur Zustände gleicher Parität sein; dies 
wurde in den Abbildungen vorausgesetzt. Die Linien zwischen den Punkten 
entsprechen den Übergängen, die in der untersuchten NTD-Näherung be- 
rücksichtigt sind. 

Als Beispiel wird ein schwach gestörter harmonischer Oszillator mit dem 
Hamiltonoperator 


A 
(1,8) 


behandelt. H läßt sich mit Hilfe von 
1 


1 


2 u (a* + a), p V? u (at a), 


y 
durch Erzeugungs- und Vernichtungsoperatoren darstellen: 


1 A 
H (at4 + 


1,10 
+ 6at? a? + data? + at + 6at? + 12ata + 6a? + 3). ae 
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In der Näherung (1,7) ist nur die Amplitude o,,(1|) von Null verschieden, 
Wegen 


[H,a*]} =ua* + (at? + 3a*2a + 3at a? + a? + 3at + 3a) (1,11) 


lautet die NTD-Gleichung 


3 
Ey O0(1|) = + 
Der Energieeigenwert ist also 
3 
Ey=u(1 (1,13) 


A 


(1,12) 


in Übereinstimmung mit der Störungsrechnung. 
In 2. NTD-Näherung berücksichtigt man die Beimischung von Zuständen 
nächsthöherer Ordnung in |%) und |o): 


|0> + |2>, | 1> = oy |1> + (1,14) 


Übergänge zwischen diesen physikalischen Eigenvektoren werden vermittelt 
durch die Operatoren 
at; a,a*?, at a2, (1,15) 
wie aus Abb. 2 hervorgeht. 
Mit 059(3|) ist im Operatorensystem (1,15) zugleich ein Matrixelement 
zwischen |3) und 10) vorhanden, das zur Bestimmung von E,, in erster 
Ordnung dienen kann (Abb. 3). 


2 


1 


0 


Abb. 2. NTD-Amplituden 9 in Abb. 3. 0, in zweiter 
zweiter Näherung NTD-Näherung 


Im Beispiel (1,10) treten in 2. Näherung zu (1,11) die H-Vertauschungen 
[H, a] [H, at]*, 


[H, a+?] = 3u at? + e [at5 + 3at4a 
+ 3at? a? + at?a? + Gar? + 9at?a + 3ata? + dat + 2a], 
[H, at = nat? + [-2atta — — Zatta (1,16) 


+ atat + a5 — 3at? — 12a? a — bat a? + 3a? — 3at], 
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hinzu. Rechts sind ebenfalls nur Ausdrücke (1,15) zu berücksichtigen. Damit 
lautet die Determinante der Gleichungen zwischen den g-Amplituden: 


e(1]) e(|1) o(3]) e(1|2) 


Sie hat die Eigenwerte 


58), Ey =-u(l +e), 


= —p(1 +32), (1,18) 


wenn man nur die ersten Glieder in der Entwicklung nach e ausrechnet. 
Ein Vergleich mit der Störungsrechnung 


(1,19) 


lehrt, daß der in ¢ quadratische Term der Energiedifferenz E,, durch (1,18) 
nicht richtig wiedergegeben wird. Die Erklärung hierfür liegt darin, daß 
(1,14) die Eigenfunktionen |1> und |5) nicht bis zur 1. Störungsordnung 
vollständig darstellt. Da H die Terme a** und a‘ enthält, treten in den Zu- 
ständen |%) bereits in der Ordnung & Beiträge |n + 4) auf, so daß die 
Übergänge 05) und 9,4 noch von der Größenordnung &? sind. Um bessere NTD- 
Eigenwerte zu erhalten, ohne von sehr hohen o-Funktionen Gebrauch machen 
zu müssen, benutzt man eine Methode, die in 2. bei der Behandlung stark ge- 
störter Systeme angegeben wird. Der Eigenwert E,, ist in der Näherung e 
richtig. 

Unter den Lösungen (1,18) sind noch die Übergänge E,, und E,, enthalten, 
da diese ebenfalls durch das Operatorensystem (1,15) dargestellt werden 
können. 

Die NTD-Gleichungen lassen sich auch in den Operatoren p und q dar- 
stellen t)2). Statt der o(k|l) benutzt man in diesem Falle die t-Amplituden 


Tno(k|l) = (q* |. (1,20) 
wobei man das symmetrische Produkt 
(1,21) 


durch Mitteilung über alle Permutationen der q und p erhält. Die q sind 
hier ebenso wie die p untereinander gleich. 

1) W. Heisenberg, Nachr. Gött. Akad. Wiss. 111 (1953). 

*) K. Symanzik, Dissertation, Göttingen 1954. 
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Das Gleichungssystem der t-Amplituden kann unter den gleichen Nahe. 
rungsannahmen wie in der at, a-Darstellung abgebrochen werden. Hierzu 
muß man die freien Erzeugungsoperatoren durch p und gq ausdrücken. Die 
freien Erzeugungsoperatoren lassen sich sofort angeben, wenn man nur Po- 
tenzen von q zu ihrem Aufbau benutzen will. Da die |n> Eigenfunktionen 
des harmonischen Oszillators sind, haben die Hermiteschen Polynome ge- 
rade die Eigenschaft 


Im = |o>. (1,28) 

Entsprechend dem Vorgehen in (1,6) kann man Matrixelemente der höheren 

H „ bei der Betrachtung von Übergängen zwischen niedrigen Niveaus vernach- 
lässigen: 

Hy (qVu) ~ 0 für N>n. (1,23) 


Man kann also das Gleichungssystem abbrechen, indem man die Hy ent- 
sprechende Kombination von t-Amplituden nullsetzt. Das Vernachlässigen 
der höheren t-Amplituden selbst steht im Widerspruch zu (1,23) und führt 
zu keiner brauchbaren Näherung. 

Da sich im allgemeinen keine NTD-Gleichungen mit ausschließlich p- 
freien t-Amplituden aufstellen lassen, genügt die Kenntnis von (1,23) nicht 
zum Abschneiden des Gleichungssystems. Es kann gezeigt werden, daß die 
von Heisenberg!) benutzten g-Funktionen eine allgemeine Darstellung 
von freien Erzeugungsoperatoren durch p und q geben. Die g-Funktionen 
gehen aus den t-Amplituden dadurch hervor, daß man zu ihnen Ausdrücke 
in niedrigeren t-Funktionen hinzufügt. Diese gewinnt man in folgender Weise: 
Aus den r(k|l) nehme man Paare von Operatoren q und p heraus und ersetze 
sie durch die entsprechenden Vakuumerwartungswerte mit negativem Vor- 
zeichen (mathematisches Vakuum, vgl. aber 2.). Wenn man die q und p nume- 
riert, heißt das 


Pno(4ı 92 - |Pı Po - +) = Tro (G1 Yo |Pı Pa — <O| (41 Tro (4s - - 

| Py —<O | (G1 92)s| - - |Pı P2--) — <O | (4s Pr)s| Tro (Ge %-- 
| Po Ps -) <O |(Py Pa)s| Tro - |Pa- — + 
<9 | (93 Ga)s| Tn0(95 - - |Pı Pe - - 

Die Vakuumerwartungswerte werden abgekiirzt durch 


A = <0| (492), |9 = <0| (Fas) 


während 
at ; 


V2 u 
verschwindet. 
Hier sind die qg und p untereinander gleich, und man kann Terme zusammen- 
fassen 


p), =<|(° *) =0 (1,26) 


(-1) 1 


+(-1) 
pres kil! 
= > =, 3) 2m)! (l—2n)! m! nit (k—2m 2m). 
(1,27) 
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Für / = 0 ist die rechte Seite von (1,27) gerade das Matrixelement der 
Hermiteschen Polynome (1,23), womit für 1] = 0 nachgewiesen ist, daß die 
g-Bedingung 

Pno(k\l) 0 für (1,28) 


der Benutzung freier Erzeugungsoperatoren zum Abbrechen des NTD-Glei- 
chungssystems äquivalent ist. Der allgemeine Beweis für beliebige k und / 
läßt sich mit Hilfe des Wickschen Theorems erbringen. 


2. NTD-Methode bei starker Störung 


In 1. wurde die NTD-Methode auf schwach gestörte Systeme angewandt, 
deren Hamiltonoperator nur um ein kleines Korrekturglied vom harmoni- 
schen Oszillator H, abweicht. Dann sind die physikalischen Erzeugungsope- 
ratoren zu H anzunähern durch die freien Erzeugungsoperatoren, die zu H, 
gehören. Bei starker Störung kann eine solche Annahme allenfalls in sehr 
hoher NTD-Näherung zu richtigen Eigenwerten führen, da in diesem Falle 
auch weit von |n> abliegende Zustände noch stark zu |%) beitragen. 

Gegeben sei der Hamiltonoperator 


>0 beliebig. (2,1) 


In diesem Falle muß ein geeignetes ,,ungestértes System“ 
H, = (2,2) 


gewahlt werden, um NTD-Methode oder Stérungsrechnung anwenden zu 
können. Dabei ist x zunächst beliebig. Aus (2,1) und (2,2) folgt als ‚‚Stör- 
glied‘“ 


Wendet man hierauf die Störungsrechnung an, so findet man 
1 34 
E, 5) + +22 +1) (2,4) 


in 1. Näherung. 

Zur Festlegung von x kann man ein Variationsverfahren durchführen. 

Unter diesem Gesichtspunkt dient die Aufspaltung (2,2), (2,3) nur dazu, ein 

Orthogonalsystem von Eigenfunktionen |nx) zu (2,2) zu erhalten. Mit diesen 

Funktionen werden die Eigenwerte 

E„(x) = <n, |H| n,> 2,5) 

von (2,1) angenähert. Da (2,5) für die richtige Eigenfunktion > von H 

ein Minimum annimmt, kann durch Variation der günstigste Wert für den 

freigelassenen Parameter x gefunden werden. Diese Variation darf für jeden 

Eigenwert einzeln durchgeführt werden, so daß man zu jedem E, ein opti- 

males x, erhält. Auf diese Weise kann gezeigt werden, daß die Annahme 
Enm 

n—- m 


(2,6) 


der Wahl eines optimalen x nahekommt. 


Nähe. 
lierzy | 
. Die 
r Po. 
ionen 
(1,22) | 
nach- 
(1,23) | 
ent. 
ssigen 
führt 
ch p- 
nicht 
| 
lung 
ionen 
rücke 
/eise: 
setze | 
Vor- 
ume- 

» 
1,24) 
1,25) 
1,26) 
men- 
-2n). 
1,2) | 


138 Annalen der Physik. 7. Folge. Band 6. 1960 


Um die NTD-Methode im Operatorensystem zu (2,2) mit (2,6) durchzu. 
fiihren, setzt man 


1 

a= p= (b+ — b) (2,7) 
und drückt (2,1) durch 6* und 6 aus. Für die Differenz der untersten Eigen. 
werte von (2,1) findet man mit E,,=x aus der Vertauschung [H, b+] die 
NTD-Gleichung 


E 8A 
09(1|) =" (1}) + gp (2,8) 
mit derselben Lésung 


V3 = 1,144, exakt: Ey = 1,087 (2,9) 


wie bei Heisenberg!). 

In p, q-Darstellung hat man bei starker Störung die Kontraktionen eben- 
falls im Operatorensystem (2,2) zu bilden, also in (1,25) u durch x zu ersetzen. 
Dieses Vorgehen stimmt in erster Näherung mit dem Vorschlag von Heisen- 
berg!) überein, die Kontraktionen mit dem physikalischen Vakuum zu 
bilden. 

In diesem Zusammenhang wird ein Verfahren zur genaueren Bestimmung 
höherer Eigenwerte vorgeschlagen. Als Beispiel wird die Energiedifferenz 
benachbarter Zustände für den Hamiltonoperator (2,1) in r-Darstellung 
berechnet. Aus den H-Vertauschungen von g und 2 findet man als Gleichungs- 
system zwischen den NTD-Amplituden 


Enn-1 Tn,n-ı (1 \) t Tn,n-ı (| 1) 

Eu.n- 1 Tn,n-ı (| 1) = id Tn,n-1 (3]). 
Die höhere t-Funktion r(3|) darf man dabei nicht durch Nullsetzen der Ampli- 
tude 9’(3|) annähern, da diese nur in Anwendung auf das Vakuum dem Er. 
zeugungsoperator b*+? entspricht. Es wäre schwierig, für beliebige n die ge- 
eignete Kombination von r-Funktionen aufzufinden, die man nullsetzen 
dürfte. Daher wird das Verhältnis von r(3|) zu r(1|) in erster Näherung 


mittels der freien Zustände zu (2,2) bestimmt. Dann läßt sich (3 |) ausdrücken 
als 


(2,10) 


7(3|) Po cn n -D + + Imn-ı Zu 3n 
t(1|) <n\q |~—1> In, n-1 2E,,n 
da 


‚ (2,11) 
-1 


at++a 


|n -)=-— (2,12) 


VeR,„.-ı 


ist. Für die Energiedifferenz benachbarter Zustände folgt dann aus (2,10) 


In,n-ı 


Ena (2,13) 


Summiert man diese Differenzen, so kann man eine viel genauere Näherung 
für E,, erhalten als durch direktes Aufstellen von NTD-Gleichungen hierfür. 
Denn das für jedes n neu gewählte Operatorensystem paßt sich dem wirk- 
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lichen Verlauf von H weit besser an als eines, das auf Grund des untersten 


Eigenwerts bestimmt wurde. Die Tabelle gibt einen Vergleich der Lösungen, 
die man auf diese Weise erhalten kann, mit den numerischen Lösungen!)3). 


Tabelle 1 


En,n-ı (2,13) E„o summiert | E,9 numerisch 


1,145 | 1,145 1,087 
1,442 2,587 2,538 
1,651 | 4,238 4,200 
1,819 | 6,057 


Die hier vorgeschlagene Methode ist auch zur Bestimmung von Eigenwert- 
differenzen E„,n-ı bei großen n brauchbar. Es ist möglich, die Untersuchung 
für den allgemeinen Fall 


+20 N gerade (2,14) 


durchzuführen. Die NTD-Methode ergibt die Beziehung 


9 (X —1]) 


Die t-Amplituden werden wieder durch Matrixelemente der ungestörten Eigen- 
funktionen des harmonischen Oszillators der Frequenz x = E„n-ı ange- 
nähert: 

Tun-ı (N — 1) & | (0 — (2,16) 


Diesen Ausdruck kann man durch eine einfache Uberlegung auswerten. Dazu 
betrachtet man in Abb. 4 alle möglichen Zwischenzustände und berücksichtigt 
(2.12). Bei n>N wird die Ausdehnung des 
Rechtecks nach unten nicht behindert. Dann 
enthält jeder Linienzug (N — 2)/2 nach oben lau- 
fende und N/2 nach unten laufende Verbindungen ra 


benachbarter Niveaus. Es gibt insgesamt (oa ) 
verschiedene Verbindungsméglichkeiten. Für 
n> N können außerdem die Unterschiede zwi- r 


schen n,n —1,...n— N vernachlässigt werden 
und man erhält 


( (2,17) Abb. 4. Zur Berechnung von 


Damit geht die Eigenwertgleichung über in | 


BE 


Enn-1 


*) L. I. Schiff, Phys. Rev. 92, 766 (1953). 
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Für große n kann andererseits die exakte Lösung auf einfache Weise mit 
Hilfe des Phasenintegrals ermittelt werden. Aus 


$ =4/2E, (“*) 


1 1 


2N 


“N+2 ya (+ 1)1(5): 


Abb. 5 bringt einen Vergleich der Koeffizienten Ky, und K, für alle X. 

Man erkennt daraus, daß das vorgeschlagene Verfahren auch bei sehr großen N 

brauchbar bleibt. Im Falle de 

2 Kastenpotentials tritt die größte 
= =247 Abweichung mit 23% auf. 


3. Mehrzeitige Theorie 


In 1. und 2. wurde die neue 
Tamm- Dancoff-Methode auf Pro- 
bleme der Quantenmechanik ange- 

- 20 -  Wandt. Die Operatoren wurden 


Abb. 5. Annäherung der Eigenwertdiffe. ‘bei als zeitunabhängig ange 
renzen En,n-ı,n > co durch NTD für ver- Sehen, die gesamte Zeitabhängigkeit 


A durch die Eigenfunktionen über- 
wird aus der kanonischen Koordi- 
nate q; mit dem diskreten Index i die Feldvariable q(r). (Die übliche Bezeich- 
nung wurde vermieden, um keine Verwechslung mit den p-Funktionen auf- 
kommen zu lassen.) In Schrödinger- und Wechselwirkungsdarstellung 
weisen die t-Amplituden für diese Feldvariable keine völlige Symmetrie 
zwischen Raum- und Zeitkoordinaten auf. Nur in Heisenbergdarstellung 
ist q(t, ¢) in der Form q(x) allein vom Vierervektor x, abhängig. Die rela- 
tivistische Kausalität ist gewährleistet, wenn die t-Amplituden auf Feyn- 
mansche Art von der Zeit abhängen. An die Stelle der symmetrisierten 
Produkte (1,20) muß daher das Wicksche 7-Produkt der Feldoperatoren 
treten. 

Nach diesen Vorüberlegungen kann die t-Amplitude für Bosefelder (für 
Fermifelder vgl. (5,2)) definiert werden als 


T„m(123...)=<n| T qg(2) g(x)... | m). (3,1) 
Diese Amplituden enthalten gegenüber (1,20) keine Impulsoperatoren; sie 


sind nicht erforderlich, denn es bestehen keine kovarianten Gleichungen zwi- 
schen t(x,) und den aus q- und p-Operatoren gemischten t-Funktionen. 
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Auch zu den neuen t-Amplituden (3,1) können g-Funktionen definiert 
werden. Der Zusammenhang bleibt der gleiche wie in (1,24), nur sind jetzt 
überall die 7T-Produkte an Stelle der symmetrisierten Produkte zu nehmen. 
Ain (1,25) geht über in die Vertauschungsfunktion der freien Felder 


Ap — 2) = <0| T |) = 4 (12), (3,2) 


wobei d(12) nur zur Abkürzung eingeführt wurde. Die ,,-Funktionen sind 
jetzt 


23...) (123...) — (12) ty9(3. . .) — d(13) Ty9(2...) —+--- 
(3,3) 
Sind die nichtlinearen Terme in der Feldgleichung groB, so kann man wie in 
2. vorgehen und statt (3,2) die Kontraktionen mit dem Vakuum zum Hamil- 
tonoperator 
2 

bilden. 

In Schrödingerdarstellung konnten Näherungslösungen der NTD- 
Gleichungen durch Nullsetzen der g-Funktionen mit höherer Variablenzahl 
gewonnen werden. Dagegen ist auch bei schwacher Störung das Nullsetzen 
der mehrzeitigen g-Funktionen nicht mehr gerechtfertigt, da es schon im 
Falle eines Freiheitsgrades falsch ist. Die p-Funktion dreier Variabler zum 
Hamiltonoperator (1,8) ist z. B.!) 


Be 
iu—(-t+t) 


sie ist also nur in der Nähe von ¢, = t, zu vernachlässigen. Die Brauchbarkeit 
der NTD-Gleichungen wird daher stark eingeschränkt, sobald g-Funktionen 
zu ungleichen Zeiten auftreten. Das ist der Fall bei der Berechnung höherer 
Eigenwerte mit Hilfe von t-Amplituden mehrerer Variabler und bei der Be- 
stimmung weiterer NTD-Näherungen für den untersten Eigenwert. Um mehr- 
zeitige 9-Funktionen zu vermeiden, muß man in diesen Fällen eine zusätzliche 
Vereinfachung dadurch einführen, daß man alle Ableitungen nach Differenz- 
koordinaten vernachlässigt und nur noch die Abhängigkeit von der Schwer- 
punktskoordinate betrachtet. Zur Bestimmung der höheren Eigenwerte wird 
es daher wie in (2,13) zweckmäßiger sein, Energiedifferenzen benachbarter 
Zustände aufzusummieren, da dann in der ersten Näherung nur t-Amplituden 
einer Variablen auftreten. 

Die p-Bedingung (3,3) ist nur gültig in Anwendung auf das Vakuum. Zur 
Bestimmung der Differenzen benachbarter höherer Eigenwerte muß man eine 
Verallgemeinerung durch die Forderung aufsuchen, daß die neue y-Bedingung 
für ein ungestörtes System exakt erfüllt ist. Das Vorgehen wird an einer g- 
Funktion von 3 Variablen erläutert. Für den harmonischen Oszillator gilt 


Tr,n-1 (123) = Vn [(m + 1) 4 


(2x 


bei ¢; >t, >t,; den Ausdruck für andere Zeiten findet man durch Vertau- 
schung der Indizes. Gesucht wird eine Darstellung von t„„-ı(1 2 3) durch r(1), 


(3,6) 
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t (2), 7(3), die für alle Zeiten gültig ist, also nicht auf t, >t, >t, beschränkt 
bleibt, d. h. die gesuchte Beziehung muß symmetrisch in allen 3 Zeiten sein, 
Der allgemeinste Ansatz ist unter dieser Voraussetzung 


Ta,n-ı (123) =a [d(12) + d(13) + d(23) t(1)] 
+ 5 [d-(12) + d-(13) + d-(23) r(1)]. 


Es erweist sich als erforderlich, neben 


(3,7) 


0 0 
d(12) = <0 =| 92 | 


= t. 
<0} 


die Vertauschungsfunktion mit umgekehrter Zeitordnung 


<0 
<0 
einzuführen, da sonst für n > 1 die Amplitude Ty, n-1 (1 2 3) nicht vollständig 
dargestellt werden kann. Für den Oszillator lautet (3,7) explizit beit, > t, >t 


Tn,n-1 (1 2 3) Vn 
(2x) 
Durch Vergleich mit (3,6) erhält man unmittelbar 
a+1= 2a, n=a-+b, n—1= 2b, (3,11) 
so daß die y-Bedingung lautet 


Gnn-1 (123) (123) fa (12) + + d(28) 


[d-(12) (3) + d-(13) 7(2) + d-(23) r(1)] = 0. 
(3,19) 


Auf die gleiche Weise findet man als p-Bedingung für 5 Variable, wenn 
man T,n-ı(12345) nur durch r(1),r(2)... auszudrücken wünscht, also 
t-Amplituden von 3 Variablen nicht vorkommen sollen: 


d-(12) = <0 TI, = 4h Zt, (3,9) 


a + eir + 


(3,10) 


Pan-1 (12345) 9-1 (12345) — 
— 2) + a(13) d(24) TB) +...) 
[d(12) d-(34) + d(13) d-(24) + ...] (3,13) 


12) d(B4) (5) + 4-13) d(24) 


ua (n — a. — 2) [d-(12) d-(34) t(5) + d~(13) d~(24) t(5)+...J= 0. 


In den eckigen Klammern stehen jeweils die 4) = - = 15 Terme, die durch 


Herausnahme von 2 Variablenpaaren zu Kontraktionen aus r(1 2 3 4 5) ent- 
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stehen. Mit Hilfe von (3,13) ist es möglich, das von Heisenberg‘)®) zur Be- 
stimmung des untersten Eigenwerts angewandte Verfahren zur Berechnung 
höherer Eigenwerte zu benutzen. 


4. Anwendung auf nichtlineares Bosefeld 


Dem anharmonischen Oszillator entspricht in der Feldtheorie das nicht- 
lineare Bosefeld 


das von Thirring®), Schiff?) und Lieb’) untersucht wurde. Schiff?) be- 
handelt das Problem (4,1), indem er eine diskrete Raumstruktur einführt. Da- 
durch werden die ö-Funktionen in der gleichzeitigen Vertauschung von p 
und q beseitigt, und es können keine Divergenzen mehr auftreten. Die Feld- 
variablen werden durch gemittelte Größen ersetzt, die nur an den Gitterpunk- 
ten definiert sind. Der Gradientterm (ög/ör)? geht in einen Ausdruck über, 
der eine Kopplung zwischen nichtlinearen Oszillatoren bewirkt. Obwohl durch 
die Einführung des Gitters der Hamiltonoperator nicht mehr invariant 
gegenüber Translationen des Raumes ist, konnte Schiff zeigen, daß Energie- 
eigenwerte mit der richtigen relativistischen Beziehung zwischen Energie 
und Impuls existieren. Diese Lösungen können daher als „Teilchen“ inter- 
pretiert werden. Vorausgesetzt wurde dabei, daß die Compton wellenlänge 
des Teilchens groß ist gegenüber der Gitterkonstanten. 

Für den Fall a = 0 werden die so gefundenen Masseneigenwerte umgekehrt 
proportional zur Gitterkonstanten. Das muß so sein, da erst durch die Einfüh- 
rung des Gitters eine charakteristische Länge in die Theorie hineingebracht 
wird. A ist dimensionslos, denn aus der Vertauschungsrelation 


[p,q] = —i d(x — 1’) (4,2) 
folgt wegen p = q, und da die 6-Funktion die Dimension 73 hat, als Dimension 
von q 

(4,3) 
J dr q* hat also bereits mit 2-1 die Dimension von H (es ist h = c = 1), so 
daß A ein Zahlenfaktor ist. 


Die Bewegungsgleichung der Operatoren folgt aus dem Hamiltonope- 
rator (4,1) zu 


(4,4) 
Deswegen lautet die Gleichung zwischen den niedrigsten t-Amplituden 
— =Ar(1 11). (4,5) 
Mit Hilfe der Greenschen Funktion der ungestérten Gleichung 
(Cy — a) @(1 2) = ö(1 2) (4.6) 
kann man (4,5) in eine Integralgleichung umformen: 
=Af @(1 2) r(2 2 2). (4,7) 


‘) W. Heisenberg, F. Kortel u. H. Mitter, Z. Naturforsch. 10a, 425 (1955). 
5) W. Heisenberg, Varenna-lectures 1958. 

6) W. E. Thirring, Helv. Physica Acta 26, 33 (1953). 

?) E.H. Lieb, Proc. Roy. Soc. London A 241, 339 (1957). 
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Da in der rechts stehenden t-Amplitude alle 3 Feldvariablen zur gleichen 
Koordinate x, gehören, darf man die zugehörige y-Funktion genähert vernach. 


lassigen : 
(2 2 2)= r(2 2 2) — 3d(2 2) 0. (4,8) 
Damit lautet die Integralgleichung fiir r(1) _— 
t(1) = 3A d(2 2) f dx, G(1 2) r(2). (4,9) 
Sucht man eine Lésung zum Viererimpuls p, so hat t(1) die Form At 
t(l)=reipz, (4,10) § Gitteı 
und man erhält aus (4,9) einfach 
34 d(1 1) G(p) = 1. (4,11) 

Für den untersten Eigenwert — das ist die Masse des leichtesten Teilchens, ee ei 
das durch die Quantisierung des Feldes (4,1) entsteht — findet man unter übers 
diesen Voraussetzungen die Beziehung aber | 

dk 1 Extra 
Dasselbe Resultat erhielt Lieb?) (2,21), mit Hilfe des Feynmanschen path- § kleine 
integral in erster Näherung. Für den niedrigsten Masseneigenwert ergibt sich f noch; 
also ein divergentes Integral. Eine Renomierung des Ausdrucks wird in dieser § Lösur 
Arbeit nicht versucht. Zumindest im Falle «4 = 0 könnte sich auch durch § larisie 
Renomierung kein sinnvoller Eigenwert ergeben, da A dimensionslos ist. Eine f bei g 
endliche Lösung kann statt dessen durch formale Regularisierung erhalten Zu 
werden. Hierdurch kommt man zu einem Vergleich mit der Methode von § Eigen 
Schiff’). Man kann nämlich die Regularisierung mit einer Masse durch- | zogen 
führen, deren Comptonwellenlänge dem Abstand zweier benachbarter Os- f 2n,n- 
zillatoren im Gitterraum gleich ist. Der Vorteil gegenüber *) besteht darin, § zu (4 
daß man die NTD-Rechnung in allen Schritten kovariant durchführen kann und 
nicht erst nachträglich die richtigen Transformationseigenschaften der Lösung Sn ho 
nachweisen muB. 

Die 1. NTD-Näherung liefert nur dann einen endlichen Wert für die Masse Auch 
des leichtesten Teilchens, wenn die Vertauschungsfunktion am Lichtkegel Zu 
nicht divergiert. Nennt man die Regularisierungsmasse M, so genügt es nicht, > 
im Impulsraum (4,18 

zu setzen, denn mit dieser Annahme bleibt d(1 1) bei k — co noch logarithmisch § wenn 
divergent. Erst nach doppelter Regularisierung mit derselben Hilfsmasse 
1 1 

d (k) (M + x2) + M?) (k2 + einfü 
(M? — (4,14) nimr 
8) Mi steh 

wird die Vertauschungsfunktion am Lichtkegel endlich. Mit der Wahl (4,14) 
wurden die beiden ersten Terme der Taylorentwicklung von d(k) um die Durc 
Stelle x? — M? gestrichen: n-teı 


d(k) = d(k)|yy + — M2) d(k)|y + (4,15) 


Aus (: 
Al 
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Aus (4,14) folgt 


d'(11) = = [ae — In]. (4,16) 
Hiermit geht die Gleichung . % für den ersten NTD-Eigenwert über in 


= 2+ — — x? In =]- (4,17) 


Abb. 6 gibt einen Vergleich mit den Ergebnissen von Schiff*). Seine 

Gitterkonstante wird durch 

l (4,18) 
mit der Regularisierungsmasse in Verbindung gebracht. Für kleine A besteht 
also eine recht gute Übereinstimmung zwischen beiden Methoden. Bei großen A 
überschreiten die Eigenwerte von Schiff die Regularisierungsmasse. Es muß 
aber bezweifelt werden, ob die £ 
Extrapolation auf Teilchen mit —2 
einer Compton wellenlänge 7 
kleiner als die Gitterkonstante 
noch gerechtfertigt ist. Die NTD- 
Lösungen erreichen die Regu- 
larisierungsmasse asymptotisch 
bei großen A. 

Zur Bestimmung der höheren 
Eigenwerte wird (3,12) herange- 
zogen. Die Eigenwertdifferenz 
E„n-ı findet man in Analogie 0 


zu (4,11) aus 10 : A 
4 Abb. 6. Niedrigster Masseneigenwert durch 
3nad’(1 1.(4,19) NTD und Gitterraumquantelung 


04 
—— ıD 


— — Gitferraum 


Auch in der Vertauschungsfunktion d’(1 1) benutzt man zweckmäßig statt 
x = Ey den gesuchten Eigenwert E,„,„-ı, da man hoffen kann, auch im feld- 
theoretischen Fall hierdurch die beste Näherung zu erzielen. Mit (4,16) geht 
(4,19) über in 


=P + (1 — n-1 (4,20) 


wenn man die Abkiirzungen 


16m’ 


einführt. Die Gleichung (4,20) wird nur für den Fall » = 0 ausgewertet. Sie 
nimmt dann dieselbe Form an wie für den untersten Eigenwert. Die Lösungen 
stehen daher in der Beziehung 

Nn,n-1 i’). (4,22) 
Durch Summation aller Eigenwertdifferenzen bis zur n-ten findet man den 
n-ten Term des Massenspektrums 


No (A) = A’). (4,23) 


(4,21) 


Ann, Physik. 7. Folge, Bd. 6 
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Es ist eine sehr bemerkenswerte Eigenschaft des Thirring-Modells mit 
a = 0, daß die Quotienten der Eigenwerte 
Eng __ "no 
Ey Me (429 
unabhängig von der Regularisierungsmasse werden. Das folgt unmittelbar 
aus der Dimensionslosigkeit von A und gilt daher in jeder Näherung. Man kann 
also bei der Betrachtung der Massenverhältnisse die Regularisierungsmas 


M — oo gehen lassen. Wäre das Thirring. 
Ano Modell (4,1) eine Gleichung für die wirklichen 
‘10 Elementarteilchen, so dürfte man tatsächlich 


nur Massenverhältnisse als Observable an. 


‘ sehen. Eine Veränderung von M ergibt nur 
eine Veränderung des gemeinsamen MaBstabs. 
[ faktors, der im System der wirklichen Ele. 
 mentarteilchen durch kein physikalische 

+N Experiment zu ermitteln ist. 
| n=3 Abb. 7 zeigt die ersten Terme des Massen- 


spektrums zum Thirring-Modell in 1. NTD. 


Pe, au Näherung. Sie nehmen bei kleinen 2’ die 
2t n=2 
Werte 
Vn’ far 27-0 (4,25) 
0 "= 
0 50 100 % an und erreichen bei großen A’ asymptotisch 
Abb. 7. Massenspektrum zum Nno __ 
Thirring-Modell für + oo. (4,26) 


5. Anwendung auf Heisenbergs Elementarteilchentheorie 


In 3. wurde ein Verfahren angegeben, aus mehrzeitigen Integralgleichungen 
auf die Differenzen höherer Eigenwerte zu schließen. Es soll jetzt versucht 
werden, diese Methode auf die Heisenbergsche Theorie*)®) mit der Be- 


wegungsgleichung der Operatoren 
+ = 0 (5,1) 
anzuwenden. 
Die t-Amplituden haben hier die Bedeutung 
Tr at! 2 | 3) ably = <n | T Ya (x) we Py | m), (5,2) 


wobei «, 6, y Spinindizes sind. Bei der Bildung der y-Funktionen können dann 
links des Strichs stehende Koordinaten (zu y gehörig) nur mit rechts stehenden 
(zu y gehörig) zu Kontraktionen verknüpft werden. Da außerdem im Mo- 
dell 4) ®) die Vertauschungsfunktion 


8) H.-P. Dürr, W. Heisenberg, H. Mitter, S.Schlieder u. K. Yamazaki, 
Z. Naturforsch. 14a, 441 (1959). 
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fir zusammenfallende Koordinaten verschwindet 
s(11) =0, (5,4) 


lautet eine spezielle, anschließend benötigte g-Funktion von 5 Variablen in 
Übertragung von (3,13) 


Yan-ı (112]1 = Tn,n-ı (11 2]12 Jagyise 
[Sae Sys t(1])e + Se (5,5) 
— [8ae 8yaT(1|)p + Spe + + 


_e-26-9 50 + 85 T(1|)al, 


_ (+1) 
6 


wobei unter s hier immer s(1 2) zu verstehen ist. 


Man kann nun durch zweimalige Anwendung von (5,1) eine Integral- 
gleichung zwischen t-Funktionen von einer und von 5 Variablen herleiten 
und anschließend (5,5) anwenden, um eine Gleichung für r(1) allein zu er- 
halten. Diese NTD-Gleichung geht aus 4) (23) bzw. ®) (109) dadurch her- 
vor, daß man die rechts stehenden Ausdrücke mit (n + 1) (m + 2)/6 multi- 
pliziert und gleichgebaute Terme mit ss-(n + 1) (n—1)/6; s~ s(n +1) 
(n — 1)/6 und s~ s~(n — 1) (n — 2)/6 hinzufügt. Zur Auswertung der Gleichung 
im Impulsraum hat man dann das in #) (29) definierte Integral X bei anderer 
Lage der Pole zu berechnen, da beim Übergang zur s--Funktion 


(x > + (x (5,6) 


übergeht. Das läßt sich neben ¢ — —e auch durch x > —x erreichen. Durch 
Anwendung der Feynmanschen Integraltricks sieht man (wenn man das an 
anderer Stelle im Integral nochmals vorkommende x vorübergehend anders 
bezeichnet), daß K bei dieser Transformation unverändert bleibt. 


Die in (5,5) auftretenden s~-Terme tragen also hier in gleicher Weise wie 
die s-Terme bei, so daß die Eigenwertgleichung für n > 1 aus derjenigen für 
n= 1 dadurch hervorgeht, daß man die rechte Seite mit 


6 (5,7) 
multipliziert. 8) (115) geht hierdurch über in 
RE. 241 - 
1=- aa X=", (5,8) 


Setzt man wie früher die in der Vertauschungsfunktion vorkommende Masse x 
der gesuchten Eigenwertdifferenz m,,n-ı gleich 


(5,9) 
10* 
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so findet man fiir die relativen Massendifferenzen 


- 


aus denen sich durch en die Massenverhältnisse 


Im = +1 6,1) 


berechnen lassen. Man erhält speziell 
n=l Ne = 1,76, Nz = 2,39. (5,12) 


Die Teilchen mit den so gefundenen Massen sind abwechselnd Fermionen und 
Bosonen, da sie sich um je ein virtuelles Teilchen unterscheiden. Diese Massen- 
eigenwerte scheinen in keiner Beziehung zu den wirklichen Elementarteilchen 
zu stehen, solange man im Anschluß an Heisenberg den niedrigsten Fermi- 
oneneigenwert mit dem Nukleon identifiziert. Wählt man dagegen das u- 
Meson als Ausgangspunkt?), so wird das Massenspektrum der Elementar- 
teilchen schon durch die 1. NTD-Näherung ziemlich gut wiedergegeben: 


Tabelle 2 


7 Mm, m/ m, 
(berechnet) (gemessen) Teilchen 
207 u 
273 
(520) 


(738) 
966 
M (1180) 
17 (1550) 
23 1836 
si , 2181 


35 2327 


39 , 2538 


Die nahe Übereinstimmung von berechneten und experimentellen Massen 
ist natürlich mit großer Vorsicht aufzunehmen, da sehr viel mehr Teilchen aus 
der Rechnung folgen als bisher beobachtet wurden. Eine endgültige Ent- 


%) K.-S. Wohlrab, Diskussionsbemerkung auf der Spezialtagung Feldtheorie, Jena 
1959. 
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scheidung über die Brauchbarkeit des Verfahrens wird man wohl erst treffen 
können, wenn es gelingt, die Lebensdauer von Teilchen zu berechnen. Die 
nicht beobachteten Masseneigenwerte müßten sich dann als sehr instabil er- 
weisen. 


Herrn Prof. Dr. W. Macke habe ich für sein Interesse an dieser Arbeit 
sehr zu danken. 


Dresden, Institut für Theoretische Physik der Technischen Hochschule. 


Bei der Redaktion eingegangen am 25. September 1959. 
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Uber die Theorie der effektiven Reichweite 
Von T. Tietz 


Inhaltsübersicht 
In dieser Arbeit leiten wir eine Formel für die Phasenverschiebung für das f Auto 
Yukawa-Potential ab. Die abgeleitete Formel für die Phasenverschiebung er- | nit 2 
laubt uns in einfacher Weise die Berechnung des Wirkungsquerschnitts der ringe 
Neutron—Proton-Streuung bei niedrigen Energien. 
Die genaue Berechnung der Phasenverschiebung fiir das Yukawa-Po- 
tential ist mühselig und muß numerisch durchgeführt werden. Eine genaue f yit! 
Übersicht über alle Arbeiten zum Thema der Berechnung der Phasenver- und 
schiebung für das Yukawa-Potential findet man in dem Artikel von Hulthén}) f yers 
und Sugawara. Wie bekannt, erlaubt uns die Phasenverschiebung die Be- 
rechnung des Wirkungsquerschnitts. In dieser Arbeit leiten wir eine Formel : 
für die Phasenverschiebung fiir das Yukawa-Potential ab. Die Schrödinger- | tg” 
Gleichung fiir das Yukawa-Potential für den Fall der Quantenzahl 7 =0 
können wir in folgender Form schreiben). Die 
mit 
z=xr, BM(ER), M=>(M,+M,), a=MBRx, 
wo E die Energie im Schwerpunktsystem, M,, die Masse des Neutrons, M, die Die 
Masse des Protons, x die Masse des Mesons und B eine Konstante bezeichnet. oe 
Die einfache analytische Näherung für y ist durch folgenden Ansatz gegeben. u 
y = sin (k x + y) — e (3) 
In der letzten Formel 7 ist die Phasenverschiebung und 5 ist eine Konstante, I 
welche wir später zu bestimmen haben. Die analytische Näherung y erfüllt ine 
folgende Randbedingungen. 
y(0)=0 und y(oco) > sin (kx + n). (4) 
Hulthen?) hat gezeigt, daß die nachstehende Schrödinger-Gleichung für Un: 
das diskrete Spektrum genau gelöst werden kann. lim 
-2 
=o, (5) 


1) L. Hulthén u. M. Sugawara, Handbuch der Physik. Band XXXIX, Berlin- 
Göttingen-Heidelberg. 1957. 
2) L. Hulthén u. B.C. H. Nagel, Phys. Rev. 90, 62 (1953). 
3) Siehe Fußnote !). 
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Für den Grundzustand bekommen wir folgenden Eigenwert a. 
a=1+ 2a. (6) 


Die entsprechende Eigenfunktion fiir den Grundzustand muB folgende Rand- 
bedingungen erfiillen. 


u(0)=0 und u(co) = 0. (7) 
Diese Eigenfunktion läßt sich wie folgt darstellen. 
=e **(1—e*). (8) 


Den Eigenwert x und die Eigenfunktion uw können wir leicht mit Hilfe der 
Methode für die Bestimmung der Eigenwerte der Schrödinger-Gleichung des 
Autors*) finden. Multiplizieren wir die D.-Gl. (1) mit w und die D.-Gl. (5) 
mit y und subtrahiert die entstehenden Gleichungen, so ergibt sich nach ge- 
ringer Umformung und Integration die folgende Beziehung 


oo 


(2? +07) | yudxr=a ( (9) 


1-e* 


Mit Hilfe der Formeln für y und u lassen sich die Integrale elementar auswerten 
und so folgt nach einigen Rechnungen folgender Ausdruck für die Phasen- 
verschiebung 7 


+2 B+ a%—a(a +b) 


+ain(2 +5+2 


(a + + 6 + 1) a+b+1) P+(a+2)? 
(1 +a)? (10) 


Die Formel (6) für den Eigenwert a = 1 + 2x vereinfacht unsere Beziehung 

und fiihrt diese auf die endgiiltige Formel fiir die Phasenverschiebung. 

+6) +6 +1) a+b+17 (a +2)? 


ctg (11) 


k . 
Die bis jetzt noch beliebige Konstante b können wir wie folgt festsetzen. Wie 


bekannt, ist die Phasenverschiebung n in der Theorie der effektiven Reichweite 
durch folgende Formel gegeben 


2 1 
ketgn =— (12) 


In dieser Formel ist d die Streulänge und r, ist die effektive Reichweite. Die 
Streulänge a läßt sich in der Gestalt schreiben 


1 
jim (ketg)=-7- (13) 
Unsere Formel (11) für n gibt uns für die Streulänge d folgende Beziehung 
1 
jim (k ctg n) = 


(@+1)(« +2) b 


‘) T. Tietz u. W. Nowak, Ann. d. Physik 1, 296 (1958). 
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Die Streulänge d, welche aus den experimentellen Messungen bekannt ist, 
bestimmt eindeutig unsere bis jetzt noch beliebige Konstante b. Die effektiv 
Reichweite r, läßt sich ebenfalls leicht aus den Formeln (11) und (12) be. 
stimmen. Die Größen n, d und r, sind nach dieser Erklärung bekannt, und 
damit erlauben sie uns die Berechnung des Wirkungsquerschnitts der Neutron— 
Proton-Streuung für niedrige Energien, welcher sich durch die Formel dar. 
stellen läßt 

4n 


Unsere Formel (11) für die Phasenverschiebung n hat den Vorteil, daß sie 


durch die richtige Bestimmung der Konstante b gute Werte für die Quer- 
schnitte für niedrige Energien der Neutron—Proton-Streuung liefert. 


in 


(15) 


sin? = 


Lödz, Polen, Institut für theoretische Physik der Universität. 


Bei der Redaktion eingegangen am 28. September 1959. 
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Eine neue elektrische Entladungsform, 
die durch Uberlagerung zweier verschiedenartiger 
Entladungen erhalten wurde 


Von @. D. Cristescu 
Mit 2 Abbildungen 


Inhaltsiibersicht 


Durch Uberlagerung eines Gleichstromes oder eines intermittierenden 
Stromes und einer Hochfrequenz-Fackelentladung gelangt man zu einer neuen 
Entladungsform, deren Eigenschaften nicht nur von den Parametern der Fackel 
allein, sondern auch vom überlagerten Strom abhängen. 

Als Spektralquelle verwendet, bietet diese Entladungsart die verschieden- 
sten Anregungsmöglichkeiten. 

Während der Untersuchung der Hochfrequenz-Fackelentladung!-®) kam mir 
der Gedanke, neues über diese Entladungsform zu erfahren, daß man der 
Fackel oder einen Teil derselben mit Hilfe eines zusätzlichen, regulierbaren 
elektrischen Feldes eine zweite Entladung überlagert und derart vielleicht sogar 
zu einer neuen Entladungsart kommt. Zu diesem Zweck wird ein Schwin- 
gungskreis (Abb. 1) im Gegentakt induktiv und darauf kapazitiv mit einem 
Resonanzstromkreis gekoppelt. Letzterer besteht aus der Spule B und dem 


Abb. 1 


1) G.D. Cristescu u. R. Grigorovici, Naturwiss. 29, 571 (1941). 

2) G.D. Cristescu u. R. Grigorovici, Bull. Soc. Roum. Phys. 42, 37 (1941). 
3) G.D. Cristescu u. R. Grigorovici, Bull. Fizica. Acad. R. P. R. 4, 419 (1952). 
4) G.D. Cristescu u. R. Grigoroviei, Com. Acad. R.P.R. 5, 515 (1955). 

5) R. Grigorovici u. G.D. Cristescu, Optika i Spektr. 6, 129 (1959). 
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Kondensator C, dessen Belegungen einerseits eine Platte oder ein leitender 
Ring, andererseits eine kleine Scheibe mit einer darauf montierten Metall. 
spitze sind. Diese Spitze dient der Entladung als Ansatzelektrode. Die Mitte 
der Spule, wo sich ein Schwingungsknoten befindet, ist durch eine Drossel mit 
dem einen Pol A der zusätzlichen Spannungsquelle in Verbindung, während 
deren zweiter Pol A’ in gleicher Weise mit einer zusätzlichen Elektrode $ in 
Verbindung steht, die zwischen die Kondensatorplatten, in verschiedenen 
Höhen, in die Fackel eingeführt werden kann. Auf diese Art wird die Unab- 
hängigkeit der beiden Stromkreise verwirklicht und die Möglichkeit ihrer 
Überlagerung nur auf das Gebiet der Entladung beschränkt. 


Der Gleichstrom der Fackel 


Mit Hilfe eines zwischen A und A’ eingeführten Meßinstrumentes kann man 
feststellen, daß die Hochfrequenz-Fackelentladung unter normalen Bedin- 
gungen im allgemeinen von einem Gleichstrom von einigen Milliampere durch- 
flossen wird, wenn der fadenförmige Kern der Entladung auch in oberem Teil 
an der Elektrode S in einem helleren Fleck endet. 

Diesen Strom scheint man einer Dissymmetrie in der Elektronenemission 
der beiden Elektroden, die die Fackel begrenzen, zuschreiben zu können. Seine 
größten Werte bei gegebener Entfernung zwischen der Spitze, auf welcher die 
Entladung ansetzt, und der Gegenelektrode S erhält man, wenn diese Elektro- 
den verschiedener Art sind, wenn ihre Oberflächen ungleich verunreinigt sind, 
oder wenn die Temperaturen der Elektrodenteile,welcheden beiden kleinen Licht- 
flecken die den Kern der Entladung begrenzen, entsprechen, verschieden sind. 


Zusätzlicher Gleichstrom in der Fackel 


Bei Anlegung einer zusätzlichen regulierbaren Spannung zwischen der 
Ansatzspitze der Entladung und der Gegenelektrode S fließt ein zusätzlicher 
Gleich- oder intermittierender Strom durch die HF-Entladung. Er führt 
zur Entstehung einer kombinierten Entladung, bei der die Möglichkeit einer 

unabhängigen DBeeinflus- 
Ina sung der verschiedenen 

Parameter besteht. 
Die Gleichstrom-Span- 
nungscharakteristiken der 
der HF-Fackel überlager- 
ten Entladung zeigen unter 
unseren Arbeitsbedingun- 
gen ein anfängliches An- 
wachsen der Spannung und 
nachher ein Abfallen der- 
selben bei wachsender 
Stromstärke. Als Beispiel 
ist in Abb. 2, Kurve a, die 
Charakteristik der Gleich- 
stromkomponenten einer 
solchen kombinierten Ent- 
ladung zwischen zwei 
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Nickelelektroden wiedergegeben. Die Entfernung zwischen den Elektroden 
betrug 35 mm, die Frequenz 4- 10° Hz und die freie Länge der Fackel in Luft 
8cm. Der negative Pol der zusätzlichen Spannung befindet sich an der An- 
satzspitze der Entladung. Eine ähnliche Kurve 6 wurde bei positiver Spitze 
erhalten. 

Bei Unterbrechung der HF-Fackelentladung bleibt oberhalb eines Mini- 
mums der Stromstärke oder Spannung eine Gleichstromglimmentladung er- 
halten, wohl infolge der durch die kombinierte Entladung bereits geschaffenen 
Ionisation der Luft. Das äußere Bild dieser Gleichstromentladung ist dem 
der HF-Fackel ähnlich ®)?). 

Die Charakteristik dieser Gleichstromentladung unter den oben ange- 
führten Bedingungen ist in Abb. 2 bei negativer Ansatzspitze mit a’ und bei 
positiver mit b’ bezeichnet. Die punktierten Linien zwischen den Kurven a 
und a’, bzw. b und b’, verbinden die Werte der Stromstärke und Spannung, 
die der kombinierten Entladung und der Gleichstromentladung allein ent- 
sprechen. 

4 Bis in die Nahe der Werte, von denen an die Gleichstromentladung selbst, 
bei Unterbrechung der HF-Spannung allein brennen kann, ist das Anwachsen 
des Stromes mit der Spannung bei der der HF-Fackel überlagerten Entladung 
beinahe linear. Bei größeren Stromstärken nimmt die Spannung der Gleich- 
stromkomponente der kombinierten Entladung beim Anwachsen des Stromes 
ab. Der Einfluß der Hochfrequenz vermindert sich mit zunehmender Strom- 
stärke. 

Eine Änderung des Elektrodenabstandes bedingt eine Veränderung des 
Kurvenverlaufs. 

Die Tatsache, daß wenigstens von einem bestimmten Wert der Stromstärke 
an die Charakteristiken der Gleichstromkomponente der kombinierten Ent- 
ladung fallend sind, führt dazu, daß beim Einschalten einer zur Entladung 
parallelen Kapazität der zur hochfrequenten Entladung überlagerte Strom 
intermittierend wird. Die Sequenz der diskontinuierlichen Entladung kann 
für eine gegebene Kapazität durch Änderung der Zusatzspannung geregelt 
werden, während die diskontinuierliche Stromstärke in der Fackel von der 
hinzugefügten Kapazität abhängt. 

Die beschriebene kombinierte Entladung kann einerseits vom Standpunkt 
ihres Mechanismus aus, sowie der erhaltenen Temperaturen, und andererseits 
als Spektralquelle und vielleicht auch dadurch, daß durch sie verschiedene 
chemische Reaktionen hervorgerufen werden können, interessant sein. 

Als Spektralquelle verwendet, bietet diese Entladungsart die Möglichkeit 
der verschiedensten Anregungen. Als sehr günstig hierzu erweist sich die Ein- 
führung eines zur Entladung parallel geschalteten Kondensators C’. 

Die Gleichstromentladung, die nach Ausschalten der Hochfrequenz bleibt, 
kann ebenfalls als Spektralquelle verwendet werden. : 

5) A. v. Engel, R. Seeliger u. M. Steenbeck, Z. Physik 85, 144 (1933). 

*) H.Thoma u. L. Heer, Z.techn. Physik 18, 465 (1933), 14, 385 (1934) ; 15, 186 (1935). 
saan’ E. Badareu, M.u. G. Giurgea u. A. Trutia, Coll. Spectr. Intern. VI, Amsterdam 
Bukarest, Physikalisches Institut der Akademie der Rumänischen Volks- 
republik. 

a Bei der Redaktion eingegangen am 14. Oktober 1959. 
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Die Energieübertragung in lumineszierenden Lösungen. II 


Lésungen mit zwei gelésten lumineszierenden Substanzen 


Von H. K. Bothe 
Mit 2 Abbildungen 


Inhaltsiibersicht 


Die in Teil I!) dieser Arbeit geschilderten Untersuchungen über den 
Energieübertragungsmechanismus in lumineszierenden Lösungen werden auf 
Lösungen mit zwei gelösten lumineszierenden Substanzen ausgedehnt. An 
Hand entsprechender Energieübertragungsschemata werden Funktionen für 
die Konzentrationsabhängigkeit der Lumineszenzintensität solcher Lösungen 
abgeleitet. Die experimentelle Prüfung dieser Funktionen an festen Lösungen 
mit Polystyrol als Lösungsmittel ergibt eine Reihe von Aussagen über den 
Mechanismus der Energieübertragung in diesen Lösungen. 


Einführung 


Wird in einer lumineszierenden Lösung eine längerwellig lumineszierende 
Substanz gelöst, so luminesziert die Lösung schon bei geringen Konzentration- 
nen der zugesetzten Substanz ausschließlich mit deren Spektrum. Dieser 
„wave-shifting‘‘ genannte Effekt läßt sich in der Weise deuten, daß die Lu- 
mineszenz der kürzerwellig lumineszierenden gelöstenSubstanz von dem länger- 
wellig lumineszierenden ‚‚wave-shifter‘‘ reabsorbiert und dieser zur Lumines- 
zenz angeregt wird. Hierdurch wird es möglich, das Emissionsspektrum 
von UV-lumineszierenden Lösungen nach längeren Wellenlängen zu ver- 
schieben und der spektralen Empfindlichkeit der gebräuchlichen Photo- 
kathoden von Photosekundärelektronenvervielfachern besser anzupassen. Diese 
Verschiebung hat jedoch nur dann praktische Bedeutung, wenn durch 
sie die Lumineszenzquantenausbeute der Lösung erhöht oder zumindest nur 
unwesentlich verringert wird. Dies hängt nicht allein von der Art des ,,wave- 
shifters“, sondern auch von den übrigen Lösungskomponenten ab. Meist 
ist das Lumineszenzvermögen von Lösungen mit zwei gelösten lumineszieren- 
den Komponenten geringer als das der Lösungen, die nur die Komponente mit 
dem besseren Lumineszenzvermögen in gleicher Konzentration gelöst ent- 
halten, in einigen Fällen jedoch auch besser (Bothe, Herforth?)). Durch 
Zusatz eines geeigneten zweiten Lösungsmittels zu einer lumineszierenden 


1) H. K. Bothe, Ann. Phys. 7 5, 339 (1960). 
2) H.K. Bothe u. L. Herforth, Referat anläßlich des internat. Kolloquiums ,,Halb- $ 
leiter und Phosphore‘‘ 1956 in Garmisch-P., ersch. Braunschweig 1958, 439. 
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Lösung kann ferner die Löschwirkung löschender Moleküle herabgesetzt 
werden (Kallmann u. Furst’)*)). Die Untersuchung dieser Effekte sowie 
die Auswahl der als wave-shifter geeigneten Substanzen und die Ermittlung 
der optimalen Konzentrationen erfolgten bisher rein empirisch (Hayes, 
Ott, Kerr°)). In der vorliegenden Arbeit wird theoretisch und experimentell 
untersucht, wie die Lumineszenzintensität von Lösungen mit zwei gelösten 
lumineszierenden Substanzen von den Lumineszenzeigenschaften der Lösungs- 
komponenten abhängt. 


Die Konzentrationsabhängigkeit der Lumineszenzintensität 


Die Konzentrationsabhängigkeit der Lumineszenzintensität J, einer 
Lösung mit einer gelösten lumineszierenden Substanz A ist bei Anregung 
mit energiereicher Strahlung gegeben durch die Funktion: 


P,{A] 


Hierbei sind P,,Q, und R, Parameter, deren Größe vor allem von der ge- 
lösten lumineszierenden Substanz A, aber auch vom Lösungsmittel abhängt 
(Kallmann u. Furst®)). Die physikalische Bedeutung dieser Parameter hängt 
davon ab, durch welche Prozesse die Übertragung der Anregungsenergie vom 
Lösungsmittel auf die gelösten Moleküle erfolgt (siehe Teil I dieser Arbeit!)). 


J (1) 


Zeigt das verwendete Lösungsmittel im reinen Zustand eine Eigenlumines- 
zenz, so ergibt sich an Stelle von (1) die Funktion 
| (1—S,) 

Ry + [A] © Qa + [A]) (Ra + 
in welcher der Parameter S, den Anteil der Lösungsmittellumineszenz an der 
Energieübertragung angibt. 


(2) 


Die Funktionen (1) bzw. (2) lassen sich ableiten auf Grund von Energie- 
übertragungsschemata, in welchen jedem an der Energieübertragung be- 
teiligten Prozeß eine gewisse Wahrscheinlichkeit zugeordnet wird (Kallmann 
u. Furst®), Birks”?)). In Abb. 1 und 2 sind derartige Schemata auf Lösungen 
mit zwei gelösten Substanzen A und B ausgedehnt. Beide Energieübertra- 
gungsschemata sind aufgestellt unter der meist erfüllten Voraussetzung, daß 
die Lösungen ausschließlich mit dem Spektrum der längerwellig lumineszie- 
renden Substanz B lumineszieren und die Konzentration der Substanzen 
A und B so gering ist, daß die gesamte anregende Strahlung vom Lösungs- 
mittel Z absorbiert wird. Neben dem direkten Energieübergang vom Lösungs- 
mittel Z zu den die Lumineszenz emittierenden Molekülen B ergibt sich ein 
weiterer Übertragungsweg von L über A nach B. Da schon bei sehr geringen 
Konzentrationen der Substanz B keine Lumineszenz von A aus der Lösung 


3) H. Kallmann u. M. Furst, Phys. Rev. 97, 583 (1955). 

4) H. Kallmann u. M. Furst, F. H. Brown, Vortrag anläßlich des internat. Kol- 
loquiums ,,Halbleiter und Phosphore“ 1956 in Garmisch-P., Braunschweig 1958, 269. 

5) F.N. Hayes, D. G. Ott u. V.N. Kerr, Nucleonics 14, 42 (1956). 

6) M. Furst u. H. Kallmann, Phys. Rev. 85, 816 (1952). 

*) J. B. Birks, Phys. Rev. 94, 1567 (1954). 
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emittiert wird, kann man annehmen, daß die Wahrscheinlichkeit des Ener. 
gieüberganges von A nach B gleich der Wahrscheinlichkeit 1/t, bzw. K,, 
der Fluoreszenzemission bei Abwesenheit der Moleküle B ist. 


Primäranregung 


Prımäranregung 


I 


Yi 
A 
Yr 8 
6 
Lumineszenzemission Lumineszenzemission 
Abb. 1. Energieübertra- Abb. 2. Energieübertragungs- 
gungsschema einer Lösung schema einer Lösung der lumi- 
der lumineszierenden Sub- neszierenden Substanzen A und 
stanzen A und Bin L. B in L. Energieübertragung 
Energieübertragung durch durch Photonen oder quanten- 
intermolekulare Stöße mechanische Resonanz 


Es mußte bei der Aufstellung der Energieübertragungsschemata berück- 
sichtigt werden, daß die Moleküle B die Wechselwirkung zwischen den Mole- 
külen A beeinflussen können. Besteht diese Wechselwirkung lediglich aus der 
Konzentrationslöschung 

A++A>A+A, 
so kann durch die Moleküle B eine zusätzliche Löschung 
At+B>A+B 
hervorgerufen werden. Dieser Löschung wird durch das additive Glied p [B] 
Rechnung getragen. 

Beruht dagegen die Wechselwirkung zwischen den Molekülen A auf einer 
Emission und Reabsorption von Photonen oder einer quantenmechanischen 
Resonanz®), so ist es zweckmäßig, das Emissionsspektrum von A in zwei Teil- 


spektren I und II zu zerlegen. Teilspektrum I wird von A nicht reabsorbiert, 
sondern mit der Wahrscheinlichkeit K, , entweder aus dem Szintillator emit- 


’) Th. Förster, Ann. d. Phys. 6, 55 (1948). 
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tiert oder bei Anwesenheit der Moleküle B von diesen reabsorbiert. Das 
kürzerwellige Teilspektrum II ist von dem Absorptionsspektrum der Substanz A 
überdeckt, so daß es außerhalb des Szintillators nicht nachgewiesen werden 
kann. Bei Anwesenheit der Komponente B tritt die Reabsorption dieses 
Teilspektrums durch die Moleküle A mit der Absorption durch die Moleküle B 
in Konkurrenz. Ein Teil der Energiequanten, welcher bei Abwesenheit der 
Moleküle B durch Reabsorption und nachfolgende innere Löschung verloren- 
geht, wird nunmehr zusätzlich von den Molekülen B absorbiert. 

Die Wahrscheinlichkeit des Energieüberganges von A nach B erhält da- 
durch die Form K,,+ K,, [B] und wird abhängig von der Konzentration [B]. 

Nach dem auf der Theorie von Kallmann und Furst beruhenden Über- 
tragungsschema (Abb. 1) ergibt sich für die Lumineszenzquantenausbeute U , , 
einer Lösung der lumineszierenden Substanzen A und B im Lösungsmittel L 
1/t, vLB] 

AB + + e[B] + «[A] + y[B] 
+ 1lt, &[4] 1/t, 
Aft, + Aft, + e[B] 1/t,+ «[A]+ y[B] 1/t; + 1/t, + B[A] + 

Legt man der Ableitung von U,, das Energieübertragungsschema Abb. 2 
zugrunde, welches für eine Energieübertragung durch Photonen oder quanten- 
mechanische Resonanz Gültigkeit hat, so erhält man 


(3) 


U,»=4 dos 4 Kız[B] 
OL Key + 2) Kııld] + + (1 — 
+ q, doz Kes Kis [A] (4) 
OL + ds) + + (1 — %ı) 


Kea + + — Ga)" 

Die relative Lumineszenzintensität J,, einer Lösung der Substanzen A 
und B im Lösungsmittel Z ist bei konstanter Anregungsenergie proportional 
der Lumineszenzquantenausbeute U,, und kann durch einen Faktor 6 zu 
dieser in Beziehung gesetzt werden 


Jap =OU 45. (5) 

In 6 geht ein die Intensität der anregenden Strahlung, die Geometrie und 
Empfindlichkeit der Meßanordnung und die willkürliche zu wählende Einheit 
der Lumineszenzintensität. Führt man die in Tabelle 1 als Funktionen der 
Wahrscheinlichkeitskoeffizienten definierten Parameter ein, so ergibt sich 
aus (3), (4) und (5) die von der Art der Energieübertragung unabhängige 


Funktion 
Is = Ry + [B] + [Bl + Qa/QılA] (6) 
P [A] 


+ + (BQ. +141 + 
oder entsprechend umgeformt 
P 

42 (Ry + (B) + + @4 FM) 

Hierbei sind die Parameter P, Q und R zugleich so definiert, daß sich für die 

Luminiszenzintensität /, bzw. J, von Lösungen der Substanzen A bzw. B 
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im Lösungsmittel L ergibt: 
P,[4] P,[B] 
A : = — 

Sie lassen sich daher aus der Konzentrationsabhängigkeit der Lumineszenz. 
intensität dieser Substanzen ermitteln. 

Der Faktor F hat bei einer Energieübertragung durch intermolekulare 
Stöße (Übertragungsschema Abb. 1) die Form 

R, 

R, + [A] + 
Bei einer Energieübertragung durch Photonen oder quantenmechanische Reso- 
nanz ergibt sich 


F=E. 


R,+@[B] 1 
R, + [41+ (10) 
wobei in beiden Fällen der Parameter E auf Grund der Definition in Tab. 1 


die Bedeutung der molekularen Lumineszenzquantenausbeute der Kompo- 
nente A hat. 


F=E. 


Tabelle 1 
Definition der Parameter in Gl. (6) und folgenden 


Energieübertragung nach Über- | Energieiibertr; nach Uber- 
Parameter tragungsschema Abb. 1 po and Abb. 2 
P,; Kea 
Qa) 
Q Ky1(1—4o1) 
to: 
Kea 
R, 1/ß(1/t; + 1/t,) 
P One on Kee 
Q 1 ; K,ı(l— 
Kes 
9») 
lt, 
+1 
K AB 
G 
a) 
Molekulare Lumineszenzquantenausbeute 
K 
der Moleküle A: 94 = Ke, 
K FB + K IB 
der Moleküle L: = 
K FL + K IL 
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Für A > 0 ergibt sich aus (6) und (8) 
P;[B] 


‘42> (Ry +(B) Qe + [B) 
und für B— 0 erhält man entsprechend 
ö P,[A] 


Durch den Quotienten 6/6’ wird die unterschiedliche Empfindlichkeit der 
Meßanordnung für die Emissionsspektren der Substanzen A und B berück- 
sichtigt. 

Sieht man zunächst von einer Beeinflussung des Lumineszenzvermögens 
der Lösung durch den wave-shifter B ab, so hat seine Verwendung nur dann 
einen Sinn, wenn 


ö 
don >1 bw 


wobei 6/0’ gegeben ist durch die relative Empfindlichkeit der Meßanordnung 
für das Emissionsspektrum der Komponente B, bezogen auf die Empfindlich- 
keit der gleichen Meßanordnung für das Emissionsspektrum der Komponente A. 
Soll Gl. (6) bzw. (7) auf Lösungen angewendet werden, deren Lösungs- 
mittel in reinem Zustand eine nicht zu vernachlässigende Lumineszenz emit- 
tiert, so muß der Einfluß dieser Lösungsmittellumineszenz auf die Energie- 
übertragung in analoger Weise, wie es im Teil I dieser Arbeit bei den Lösungen 
mit nur einer gelösten lumineszierenden Substanz geschehen ist, durch die 
Einführung der Parameter S, und S, berücksichtigt werden. S, bzw. S, 
geben nach der in Teil I dieser Arbeit gegebenen Definition den Teil der vom 
Lösungsmittel an die Moleküle A bzw. B abgegebenen Energiequanten an, 
für welche die Wahrscheinlichkeit K, , der Selbstabsorption durch das Lésungs- 
mittel Null ist, so daß er bei Abwesenheit der Moleküle A und B aus dem 
Lösungsmittel emittiert wird. 
AB Ry + + B] + + 1B] 
Rz + [B] + © Qe + + 
Diese Funktion hat für A — 0 bzw. B — 0 die gleichen Grenzwerte (11) bzw. 


(12) wie Gl. (6), wobei allerdings die Funktionen J, bzw. J, nunmehr die Form 
der Gl. (2) haben. 


(13) 


Kriterien für die Art der Energieübertragung zwischen den verschiedenartigen 
gelösten Molekülen 


Gl. (6) bzw. (13) enthält bis auf die Parameter E und @ die gleichen Para- 
meter, durch welche auch die Konzentrationsabhängigkeit der Lumineszenz- 
intensität der Lösungen gegeben ist, die nur je eine der beiden lumineszierenden 
Verbindungen im gleichen Lösungsmittel gelöst enthalten. Von Kallmann 
und Furst wurden diese Parameter bereits für zahlreiche gut lumineszierende 
organische Verbindungen in flüssigen Lösungsmitteln ermittelt®). Über die 
Bestimmung der entsprechenden Parameter für polymerisierte Lösungen wurde 
in Teil I dieser Arbeit berichtet. Die beiden einzigen, nicht aus der Konzen- 
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trationsabhängigkeit einfacher Lösungen bestimmbaren Parameter E und 4 
sind in dem Faktor F (Gl. 9) bzw. (10) der Gl. (6) und (13) enthalten. F er. 
gibt sich nach Gl. (13) aus der Lumineszenzintensität /,, einer Lösung der 
Konzentrationen A und B und den Parametern Ry, Q,, S4, Pg, Op: Rp 8, 
der beiden gelösten Substanzen A bzw. B zu: 


_ P 


— [48 
@-¢ R, + [B]’ 
Qs[A] + Q4[B] © Q4 Qs + + 
S,Q,[A] (1 — S4) Q4[B] 


Ist F = F, für eine Lösung der Konzentrationen [A,] und [B,] und F =f, 
für eine Lösung der Konzentrationen [B,] und [A,], so lassen sich die Para. 
meter E und G aus den Werten von F,, F, und R, ermitteln. 

Es ergibt sich 
G — Fi(Ra + [AD — Ba + (15) 
(F 27 F, ı) [By] 

Die Form der Funktion von E hängt von der Art der Energieübertragung ab. 
Aus Gl. (9) (Energieübertragung nach Schema Abb. 1) erhält man 


‘A 
und aus Gl. (10) (Energieübertragung nach Schema Abb. 2) 
F(R,+[A]+ @[B 
D, (17) 
Ein Kriterium dafiir, welcher der beiden Werte und damit welches der 
beiden Energieübertragungsschemata den gegebenen Verhältnissen am besten 


entspricht, ist dadurch gegeben, daß E als molekulare Lumineszenzquanten- 
ausbeute konzentrationsunabhängig, positiv und kleiner als 1 sein muß. 


(16) 


Messung der Lumineszenzintensität von Lösungen 
mit zwei lumineszierenden Komponenten 


Die experimentelle Prüfung der abgeleiteten Funktionen wurde an Poly- 
styrollösungen durchgeführt, welche die blau lumineszierenden Verbindungen 
Anthrazen oder Fluoranthen als wave-shifter und die ultraviolett lumines- 
zierenden Substanzen Stilben, Naphthalin, p-Terphenyl oder Phenyl-bi- 
phenylyl-oxadiazol (PBD) als kürzerwellig lumineszierende Komponente ent- 
hielten. Die Konzentrationen beider Komponenten wurden zwischen 0,01 
und 0,1, im Falle des PBD zwischen 0,003 und 0,03 Mol pro Liter Styrol 
variiert. In diesen Konzentrationsbereichen zeigen die betreffenden Substan- 
zen den steilsten Anstieg der Lumineszenzintensität mit der Konzentration 
(siehe Teil I dieser Arbeit). Die Herstellung der Lösungen und die Messung 
der relativen Lumineszenzintensität erfolgten in der gleichen Weise und unter 
den gleichen Bedingungen, wie in Teil I dieser Arbeit beschrieben. Hierdurch 
wurde es möglich, die in diesem Teil der Arbeit ermittelten und in Tab. 2 noch 
einmal zusammengestellten Parameter P, Q, R und S zur Berechnung der 
Parameter @ und E in Gl. (13) heranzuziehen. 
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Auswertung der MeBergebnisse 


Die in vorstehendem Kapitel beschriebenen Messungen wurden an solchen 
Lésungen durchgefiihrt, welche Substanzen enthielten, deren Parameter P, 
Q, R und S auf Grund der in Teil I dieser Arbeit beschriebenen Messungen 
und Berechnungen bekannt waren (Tab. 2). Die noch unbekannten und von 


Tabelle 2 
Die Parameter P,Q, R und S der Konzentrationsabhängigkeit der Lumineszenzintensität 
in Polystyrol bei Anregung mit £-Strahlung (siehe auch Teil I dieser Arbeit) 


Gelöste lumineszierende Parameter 
Substanz P Q R 
[1] [Mol/Lit]?) [Mol/Lit] [Mol/Lit] | 

Fluoranthen | 201 ' 0,01818 0,4675 02 
Anthrazen | 955 | RATE 2,231 0,19 
Stilben | 62 0,01632 0,1727 | 0,26 
Naphthalin 277 0,1237 0,688 3 
p-Terphenyl | 346 | 0,02482 0,282 0,136 
Phenyl-biphenylyloxadiazol | 

(PBD) | 471 ' 0,01286 0,3724 0,11 


der Art der Energieübertragung zwischen den gelösten Molekülen abhän- 
genden Parameter E und @ sind in den Faktoren F Gl. (9) bzw. (10) ent- 
halten, welche nach Gl. (14) aus den Parametern P, Q, R und S, den Konzen- 
trationen [A] und [B] der beiden gelösten lumineszierenden Substanzen A 
und B und der gemessenen relativen Lumineszenzintensität /,, berechnet 
wurden. Aus je zwei Werten von F bei gleicher Konzentration [B] und unter- 
schiedlicher Konzentration [4] wurden die Parameter @ ermittelt. Da F 
jeweils für drei verschiedene Konzentrationen [A] bestimmt wurde, ließ sich 
G für jede Konzentration [B] durch drei voneinander unabhängige Gleichungen 
berechnen. Die sich hierbei ergebenden Werte mußten bei Gültigkeit von 
Gl. (15) übereinstimmen. Dies war tatsächlich ausnahmslos der Fall, so daß 
für die weiteren Berechnungen die Mittelwerte genommen werden konnten. 
Aus G, F und R, ergeben sich nach Gl. (16) und (17) die Parameter E, und 
E,. Die gemessenen Werte der relativen Lumineszenzintensität J,, und 
die aus ihnen berechneten Parameter F,@, E, und E, sind in den Tabellen 
3—6 zusammengestellt. Die Parameter @ und E sind nach Tab. 1 so definiert, 
daß sie bei Gültigkeit von Gl. (13) unabhängig von der Art und der Konzen- 
tration [B] des wave-shifters B sein müßten. Die Parameter E müssen außer- 
dem als molekulare Lumineszenzquantenausbeute einen Wert zwischen 0 
und 1 haben. Diese Bedingungen werden bei den Lösungen mit Stilben als 
kürzerwellig lumineszierende Komponente A (Tab. 3) von den Parametern @ 
und E,„ erfüllt. Das bedeutet, daß die Konzentrationsabhängigkeit der 
Lumineszenzintensität dieser Lösungen der Funktion (12) entspricht, wobei 
der Faktor F die dem Energieübertragungsschema Abb. 2 (Energieübertra- 
gung durch Photonen oder quantenmechanische Resonanz) entsprechende 
und in Gl. (10) angegebene Form hat. 

Bei den Lösungen mit p-Terphenyl, Naphthalin bzw. PBD als kürzer- 
wellig lumineszierende Komponente A besitzen nicht die Parameter E,, son- 
dern die nach Gl. (16) berechneten Parameter E, die geforderte Größen- 
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Tabelle 3—6 

Die Konzentrationsabhängigkeit der Lumineszenzintensität von Polystyrol-Lösungen mit 

zwei gelösten lumineszierenden Substanzen und die aus ihr nach Gl. (14)—(17) berech- 
neten Parameter F,G, E, und E, 


Tabelle 3 
| Relat. 
Lumineszierende Konzentrat. | Lum.- 
Substanzen in Mol/L | inten- Parameter (berechnet) 
| sitat 
Stilben | Fluoranthen 0,01 | 0,01 211 | 0,658 14,2 | 1,24 0,68 
0,038. 231 | 0,617 
| 01 | 210 | 0515 | 
0,01 | 0,03 284 0,738 | 11,4 | 2,24 | 0,752 
0,03 | 283 0,71 
| 0,1 | 261 0,63 | | 
| 0,01 | 0,1 311 | 04 «| 12,4 | 61 | 0,76 
| 0,03 | 300 0,73 | 
| 01 | 279 0,695 | 
| Anthrazen 0,01 | 0,01. 170 0,676 | 124,12 | 0,7 
0,03 207 0,634 | 
01 | 202 0,522 
0,01 10,03 | 187 0,65 | 10,1 | 1,83 | 0,663 
0.03 | 214 0,623 | 
101 | 216 | 0,5 | | 
0,01 | 0,1 24 | 0,708 | 11,3 5,37 0,712 
0,03 258 | 0,699 
| 0,1 | 960 | 0,659 | 
Mittelwert 12 | 0,709 
Mittlerer Fehler d. Mittelwertes: absolut 1083| | 00182 
m % 2 | 2 
Mittlerer Fehler d. Einzelbestimmung: absolut | 0,56 0,035 
in % | 4,7 5 


ordnung (kleiner als 1) und Unabhängigkeit von der Art und Konzentration 
des wave-shifters B. Die Konzentrationsabhängigkeit der Lumineszenz- 
intensität dieser Lösungen wird daher ebenfalls durch Gl. (13) wiedergegeben, 
wobei jedoch der Faktor F die aus dem Energieübertragungsschema Abb. 1 
(Energieübertragung durch intermolekulare Stöße) abgeleitete Form GI. (9) 
hat. Die Parameter @ sind bei diesen Lösungen ausnahmslos negativ, und 
nicht sie selbst, sondern ihre Produkte mit der Konzentration [B] sind un- 
abhängig von der Art und Konzentration der wave-shifter-Komponente B. 
Diesem Umstand wird Rechnung getragen, wenn man statt Gl. (9) schreibt 


Ry 


F =E 1, 


(18) 


wobei @’ identisch ist mit dem absoluten Wert der Produkte @ [B] in den 
Tabellen 4—6. 
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‚Tabelle 4 
mit | Relat. 
rech- Lumineszierende Konzentrat. |Lum.- 
Substanz | in Mol/L |inten- Parameter (berechnet) 
| | sität 
— A G G[B) | #, | E, 
) Naphtha- | Fluor- ı 0,01 | 0,01 | 215 1,08 | —42,1 —0,421 (0,435 1,16 
lin anthen | 0,03 231 | | 
0,1 (245 0,817, | 
0,01 | 0,03 |290 1,335) —14,3 | —0,4295 0,523 | 1,39 
68 0,03 300 1,248 
0,1 310 1,005 
0,01 | 0,1 | 313 | 1,832) —4,586 | —0,4586 0,458 | 1,37 
a 0,03 314 | 1,22 | 
a 0,1 313 0,975 
Anthrazen | 0,01 | 0,01 |127 | 1,02 | —42,5 | —0,425 0,404 1,055 
Im 0,03 153 | 0,94 | 
0,1 197 0,765 
0,01 | 0,03 | 165 1,025 | —13,7 —0,411 |0,427 | 1,06 
7 0,03 186 0,958 | 
, 0,1 210 | 0,781 | 
0,01 | 0,1 | 2395/11 | —4,43 | —0,443 |0,408 | 1,15 
663 0,03 251 1,02 | 
2 0,1 258 0,82 | 
Mittelwert —0,431 | 0,4425 1,19 
‚712 Mittlerer Fehler d. Mittelwertes: absolut 0,007 0,018 | 
in % 1,6 | 4,1 
re Mittlerer Fehler d. Einzelbestimmung: absolut 0,017 004 | 
‚709 in % 4 1 
‚0142 
| Bei Gültigkeit von Energieübertragungsschema Abb. 1 ist@ = p/p (Tab. 1). 
035 Da @ für die betreffenden Lösungen negativ ist, kann man substituieren: 
| = —G [B] = = pP. (19) 
ration Hieraus folgt, daß die Konzentrationsabhängigkeit der Lumineszenzintensität 
szenz- der in Tab. 4—6 genannten Lösungen dann der aus dem Energieübertragungs- 
yeben, schema Abb. 1 abgeleiteten Funktion entspricht, wenn in diesem Übertra- 
\bb.1 gungsschema statt » [B] geschrieben wird —q’. Hierbei ist es zunächst gleich- 
31. (9) gültig, ob man gy’ der Wahrscheinlichkeit 6 [A] zuordnet, wie es dem Energie- 
‚ und übertragungsschema Abb. 1 entsprechen würde, oder ob man es von der 
d un- Wahrscheinlichkeit 1/t; der inneren Léschung abzieht. Im ersten Falle ergibt 
ite B. sich ; 
ıreibt = = B (20) 
(18) Die in den Tab. 4—6 enthaltenen Werte G’ =@[B] sind durchweg größer 
als die Konzentration [A] der betreffenden Lésungen (A < 0,1 Mol/1). Es 
den würden sich daher nach Gl. (20) negative Werte fiir die Wahrscheinlichkeit 
n de 


der Konzentrationsléschung der Molekiile A ergeben. Da dies physikalisch 


sun 
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Tabelle 5 
| = 
K t ‚Lum.. 
in Tinten. Parameter (berechnet) 
| sität | 
p-Ter- Fluor | 0,01 | 0,01 | 268 | 1,261| —11,6 | —0,116 0,786 | 1,34 
phenyl anthen | 0,03 318 | 1,13 
01 305 | 0,873 
0,01 | 0,03 1300 | 1,06 | —4,48 | —0,1334 0,59 | 1,18 
0,03 '309 | 0,935 | | 
0,1 273 | 0,68 
0,01 0,1 | 312 | 0,934| —1,345 | —0,1345 0,52 | 0,996 
0,03 | 0,81 | 
| 01 | 271 | 0,596 | 
Anthrazen | 0,01 | 0,01 274 | 1,63 | —10,26 | —0,1026 1,09 | 1,72 
0,03 1384 | 1,475 | 
0,1 388 11 | | 
0,01 01 | 0. 0,03 246 | 1,305) —3,69 | —0,1108 0,885 | 1,38 
0,03 | '313 | 1,17 | 
01 ‚310 | 0,871) 
0,01 1268 | 1,11 -1,02 | —0,102 0,746 | 1,17 
| 0,03 298 1,0 | | | 
01 | 1279 | 0,758) | 
Mittelwert — 0,1165 | 0,77 1,29 
Mittlerer Fehler d. Mittelwertes: absolut 0,006 ‚0,085 
in % 6 11 
Mittlerer Fehler d. Einzelbestimmung: absolut 0,015 |0,2 
in % 15 26 


keinen Sinn hat, kann ’ nur der Wahrscheinlichkeit 1/t, der inneren Löschung 
zugeordnet werden. Das bedeutet, daß die Wahrscheinlichkeit der inneren 
Löschung der Moleküle A durch die Anwesenheit der Moleküle B herab- 
gesetzt wird. Durch was für einen Prozeß dies geschehen könnte, läßt sich 
an Hand der bisherigen Vorstellungen von der inneren Löschung nicht sagen, 
zumal der betreffende Prozeß nicht mit einer Energieübertragung auf die 
Moleküle B verbunden sein kann, sondern lediglich zu einer Erhöhung der 
molekularen Lumineszenzquantenausbeute E = i, + BER... +h der Molekiile A 
fiihrt. Der in Frage kommende ProzeB ist in dem untersuchten Konzentrations- 
bereich (B = 0,01 — 0,1 Mol/l bzw. bei PBD: B = 0,003 — 0,03) unabhängig 
von der Konzentration [B]; er muß aber bei geringeren Konzentrationen [B] 
konzentrationsabhängig sein, denn bei einer reinen Lösung der Substanz A 
tritt der Parameter G’ = p’/ß nicht auf; 9’ muß also für B = 0 verschwinden. 


Interessant ist ein Vergleich der sich aus der Konzentrationsabhängigkeit 
der Lumineszenzintensität ergebenden Parameter E mit den Werten der mole- 
kularen Lumineszenzquantenausbeute, welche Birks auf Grund der an kri- 
stallinen Substanzen gemessenen Lumineszenzabklingzeiten ermittelte (Tab. 7). 
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Tabelle 6 
Lumineszierende Konzentrat Lum.- 
| in Mol/L * | inten- Parameter (berechnet) 
sität 
[4] | [BJ | lun | F G GEB) | B, | F, 
PBD Fluor- 
(Phenyl-| anthen | 0,003) 0,01 | 442 |4,7 | —122,6 | —0,3675|0,94 | 71,5 
biphe- | 0,01 452 | 2.96 | 
nylyl- | 0,03 308 | 0,94 | 
di | 
| 0,033, 0,03 |377_ | 3,8 | —36,6 | —0,366 |1,02 | 59,4 
| 0,01 446 | 2,54 | 
| 0,03 327 | 0,98 | 
0,003 0,1 | 327 | 2,64 | —120 | —0,36 |1,17 35 
0,01 346 | 1,92 
| 0,03 318 | 
Mittelwert —0,3645 1,04 
Mittlerer Fehler d. Mittelwertes: absolut 0,0028 0,073 
in % 0,8 7 
Mittlerer Fehler d. Einzelbestimmung: absolut 0,0045 0,115 
in % | 


Da der Parameter E ebenfalls die Bedetung der molekularen Lumineszenz- 
quantenausbeute hat, und diese als molekulare Eigenschaft nicht wesentlich 
von den äußeren Bedingungen abhängen darf, müßte er mit dem von Birks’) 
für die gleiche Substanz auf völlig andere Weise ermittelten Wert überein- 
stimmen. Derartige vergleichbare Werte existieren für die Substanzen Stilben, 
Terphenyl und Naphthalin. 
Tabelle 7 
Vergleich der von Birks’) ermittelten molekularen Lumineszenzquantenausbeuten mit 
den Parametern 


Molekulare Lumineszenzquantenausbeute 
Lumineszierend E berechnet auf Grund der | re 
Lumineszenzintensität von | von 
Lösungen | a 

Stilben | 0,71 | 0,6 
p-Terphenyl 0,77 | 0,6 
Naphthalin | 0,44 | 0,9 
PBD 1,0 


Wie aus Tab. 7 hervorgeht, stimmen die Werte fiir die auf verschiedene 
Weise ermittelten molekularen Lumineszenzquantenausbeuten von Stilben 
und p-Terphenyl verhältnismäßig gut überein. Hieraus geht hervor, daß die 
Definition des Parameters E richtig ist, und daß die molekulare Lumineszenz- 
quantenausbeuten von Stilben und p-Terphenyl nicht wesentlich vom Aggre- 
gatzustand (kristallin oder gelöst) abhängen. Naphthalin hat dagegen in 
Lösung offenbar nur eine halb so große molekulare Lumineszenzquantenaus- 
beute wie in kristallinem Zustand. Dies steht in guter Übereinstimmung mit 
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der Beobachtung von Herforth®), nach welcher das Lumineszenzvermögen 
eines Naphthalin-Kristalls beim Übergang in den amorphen Zustand auf die 
Hälfte abnimmt. Für PBD, eine der Substanzen mit dem besten Lumines. 
zenzvermögen in Lösung, ergibt sich eine molekulare Lumineszenzquanten. 
ausbeute von eins. 

Zusammenfassung der Ergebnisse 


Die Konzentrationsabhängigkeit der Lumineszenzintensität von Lösun- 
gen mit zwei gelösten lumineszierenden Substanzen läßt sich in gleicher Weise 
an Hand der entsprechenden Energieübertragungsschemata ableiten, wie die 
durch Kallmann und Furst bzw. Birks bereits für Lösungen mit nur einer 
gelösten lumineszierenden Substanz geschehen ist. Die sich hierbei ergebende 
Funktion (13) enthält die Parameter P, Q, R und S, durch welche auch die 
Konzentrationsabhängigkeit der Lumineszenzintensität der Lösungen charak- 
terisiert ist, die nur je eine der beiden betreffenden lumineszierenden Substan- 
zen enthalten. Hinzu kommen noch die beiden für die kürzerwellig lumines. 
zierende gelöste Komponente spezifischen Parameter E und @, von denen E 
die Bedeutung der molekularen Lumineszenzquantenausbeute dieser Kompo- 
nente hat. Der Parameter @ tritt entweder nur im Nenner oder im Zähler 
und Nenner eines Faktors F der betreffenden Funktion auf, je nachdem ob die 
Energieübertragung durch intermolekulare Stöße oder durch Photonen bzw. 
quantenmechanische Resonanz erfolgt. 

Die experimentelle Prüfung der zunächst theoretisch abgeleiteten Funktion 
ergab, daß ihr die Konzentrationsabhängigkeit von Polystyrol-Lösungen mit 
zwei gelösten lumineszierenden Substanzen in dem betrachteten Konzentra- 
tionsbereich innerhalb der Genauigkeit der Messungen und Berechnungen 
völlig entspricht. Bei den Lösungen, die Stilben als kürzerwellig lumineszie- 
rende Komponente enthalten, hat die Konzentrationsabhängigkeit der Lu- 
mineszenzintensität die Form, welche bei einer Energieübertragung durch Pho- 
tonen oder quantenmechanische Resonanz zu erwarten ist. Bei den übrigen 
suchten Lösungen stimmt die Konzentrationsabhängigkeit der Lumineszenz- 
intensität mit der für den Fall einer Energieübertragung durch intermolekulare 
Stöße abgeleiteten Funktion überein, wenn diese bezüglich des Vorzeichens 
und der Konzentrationsabhängigkeit des Parameters G modifiziert wird. 

Die sich aus den Berechnungen ergebenden Werte für die molekularen 
Lumineszenzquantenausbeuten einiger Substanzen stimmen mit den von 
Birks’) auf andere Weise ermittelten Werten gut überein. 


9) L.Herforth, Dissertation TH Berlin 1948; L.Herforth u. H. Kallmanı, 
Ann. Phys. 4, 231 (1949). 
Leipzig, Institut für angewandte Radioaktivität. 


Bei der Redaktion eingegangen am 16. Oktober 1959. 
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Zur Frage des thermischen Gleichgewichts im Eisenbogen 
Von K. H. Krysmanski 


Inhaltsiibersicht 


Entgegen der weitverbreiteten Ansicht, daß in Lichtbögen, die in freier 
Atmosphäre brennen, das thermische Gleichgewicht praktisch stets erfüllt ist, 
wird gezeigt, daß in Eisenbögen hinsichtlich der Besetzungszahlen der einzelnen 
Niveaus beträchtliche Abweichungen von der Boltzmann-Verteilung auf- 
treten. Das hat zur Folge, daß Temperaturmessungen unter Zuhilfenahme von 
Übergangswahrscheinlichkeiten zu fehlerhaften Werten führten. Andererseits 
sind daher alle Bestimmungen der Übergangswahrscheinlichkeiten aus Linien- 
intensitäten mit einem systematischen Fehler behaftet. 


Einleitung 


Die Bestimmung von Bogentemperaturen ist gegenwärtig noch immer als 
ein unbefriedigend gelöstes Probelm anzusehen, besonders bei Metallbogen 
von niedriger Stromstärke (1—10 A), die in freier Atmosphäre brennen. Die 
Schwierigkeiten liegen nicht nur auf experimentellem Gebiet, obwohl es hier 
Meßverfahren gibt, die unter Umständen Temperaturunterschiede von 
weniger als 100° zu messen gestatten, sondern sind vor allem in der Theorie 
zu suchen. Es geht dabei um die Frage, ob sich für ein heißes Bogengas eine 
sinnvolle Temperaturdefinition geben läßt oder nicht. Strenggenommen 
ist der Temperaturbegriff nur für den schwarzen Körper anwendbar, da in 
diesem Falle alle temperaturabhängigen physikalischen Vorgänge zum selben 
Zahlenwert der Temperatur führen. Unter dieser Voraussetzung spricht man 
von einem thermischen Gleichgewichtszustand. Die Bedingungen hierfür 
lauten in allgemeinster Form: Für eine gegebene Temperatur hat jedes an 
der Wechselwirkung beteiligte Teilchen einen bestimmten Mittelwert der 
kinetischen bzw. potentiellen Energie, die Gesamtheit der Energiewerte unter- 
liegt einem bestimmten Verteilungsgesetz, die Strahlung ist eine eindeutige 
Funktion der Wellenlänge. 

Der Frage des thermischen Gleichgewichts in Bogensäulen ist in den letzten 
30 Jahren viel Aufmerksamkeit geschenkt worden, was in der Literatur ent- 
sprechenden Niederschlag gefunden hat. Wegen der bei Gasen im Vergleich 
zu festen Körpern wesentlich geringeren Wechselwirkung kann in einer Bogen- 
entladung das thermische Gleichgewicht niemals streng erfüllt sein. Da die 
Energieaufnahme hauptsächlich durch die Elektronen erfolgt, die Energie- 
abgabe dagegen im wesentlichen durch die Gasatome (Strahlung, Konvek- 
tion) muß die Elektronentemperatur stets größer als die Gastemperatur sein. 
Fernerhin ist die Abweichung der Strahlung von der eines schwarzen Körpers, 
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wenn man von einzelnen sehr starken Linien mit Selbstabsorption absieht, 
offensichtlich. Es ist daher lediglich eine Frage der Zweckmäßigkeit, ob man 
den Zustand des Bogengases trotzdem noch als thermisches Gleichgewicht 
bezeichnen darf oder nicht. Es ist also zu untersuchen, ob die Abweichung 
vom thermischen Gleichgewicht so groß ist, daß sie über die experimentelle 
Meßgenauigkeit hinausgeht oder bei theoretischen Ansätzen berücksichtigt 
werden muß. 


Für ein thermisches Gleichgewicht muß nicht nur die Gleichheit von Gas- 
und Elektronentemperatur gefordert werden, sondern es müssen auch Anre- 
gungstemperatur und Gleichgewichtstemperatur des lIonisationszustandes 
mit der Gastemperatur übereinstimmen, ein Punkt, der häufig übersehen 
wurde. Bei thermischem Gleichgewicht erfolgt die Besetzung der angeregten 
Zustände nach der Boltzmann-Formel. Dabei führen Stöße von Elektronen, 
Ionen und Atomen und Strahlungsabsorption zu höheren, Stöße zweiter Art 
und Strahlungsemission zu tieferen Zuständen. Für ein Plasma, das von 
einer Hohlraumstrahlung gleicher Temperatur durchsetzt wird, sich also im 
vollen thermischen Gleichgewicht befindet, sind Emissions- und Absorptions- 
prozesse, wie auch die Stoßprozesse unter sich gleich groß. Bei einem Fehlen 
der Hohlraumstrahlung sind die Emissionsprozesse zahlreicher als die Ab- 
sorptionsprozesse, was dazu führt, daß die Zahl der angeregten Teilchen 
kleiner ist, als es der Boltzmann -Formel entspricht, oder was gleichbedeutend 
ist, daß die Anregungstemperatur kleiner als die Gas- bzw. Elektronentem- 
peratur ist. Für sehr starke Linien mit nicht zu hohen oberen Termen, ins- 
besondere Resonanzlinien, kann das Plasma eine ausreichende Tiefe besitzen, 
so daß in seinem Innern die Emissions- und Absorptionsprozesse im Gleich- 
gewicht sind. Bei schwachen Linien dagegen erfolgt die Emission aus optisch 
dünner Schicht, d. h. ohne Selbstabsorption. Wie stark die Abweichung von 
der Boltzmann-Verteilung ist, läßt sich bei den lückenhaften Kenntnissen 
über die Eigenschaften des Plasmas, besonders hinsichtlich der Wirkungs- 
querschnitte kaum berechnen. Ansätze hierzu finden sich bei Maecker!), 
wo allerdings bei Stoßvorgängen nur Elektronen berücksichtigt werden. 
Danach ist die Abweichung von der Boltzmann -Verteilung nur klein, wenn 
bei einer großen Zahi von Anregungsstößen der größte Teil durch Stöße 
zweiter Art rückgängig gemacht wird. Das ist um so mehr der Fall, je größer 
die Elektronenkonzentration und die Wirkungsquerschnitte für Anregung 
und Stöße zweiter Art sind und je kleiner die Summe der Übergangswahr- 
scheinlichkeiten von Linien des betreffenden Terms ist. Im gleichen Sinne 
wirkt sich eine stufenweise erfolgende Anregung aus. Nach Ornstein’) 
spielen allerdings die Elektronenstoßprozesse bei niedrigeren Temperaturen, 
und dazu zählen auch noch die Temperaturen im Niederstrombogen, eine 
untergeordnete Rolle im Vergleich zu Stößen von neutralen Atomen und 
Molekülen. 


Ähnliche Überlegungen treffen für die Gleichgewichtstemperatur des 
Ionisationszustandes zu. Die Saha-Formel gilt ebenfalls nur für ein thermi- 
sches Gleichgewicht. Eine quantitative Überprüfung der Anwendbarkeit 


1) H. Maecker, Handbuch der Physik, herausgeg. v. S. Flügge, Bd. XXII, S. 305 
(1956). 
2) L.S. Ornstein u. H. Brinkman, Physica 1, 797 (1934). 
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dieser Formel für die Ionisierung im Bogen läßt sich kaum durchführen, da 
die Bestimmung von Partialdrucken und effektiver Ionisierungsspannung 
bei zusammengesetzten Bogengasen auf größte Schwierigkeiten stößt. 

In fast allen Arbeiten, die sich mit der Temperaturbestimmung von 
Bogenentladungen befassen, wird angenommen, daß die Bedingung des ther- 
mischen Gleichgewichts praktisch erfüllt ist. Diese immer wieder übernommene 
Behauptung geht auf die Untersuchungen von Ornstein?), Mannkopf?°), 
Witte’) u.a. zurück. Während Ornstein den thermischen Mechanismus 
im Bogen theoretisch betrachtet, erbringen Mannkopf und Witte für den 
frei brennenden Kohlebogen den Nachweis der thermischen Anregung und 
schließen daraus auf die Gleichheit von Elektronentemperatur und Gastem- 
peratur. Diese Bedingung ist zwar notwendig, aber noch nicht hinreichend 
für ein thermisches Gleichgewicht im obigen Sinne. 

Für den Eisenbogen hat Burhorn®) den Beweis des thermischen Gleich- 
gewichts zu erbringen versucht, indem er die Gastemperatur aus der Doppler- 
breite bestimmter Linien und die Elektronentemperatur aus der Kontinuums- 
intensität ermittelt und beide Temperaturen gleich groß findet. Die Tem- 
peraturbestimmung aus der Dopplerbreite ist nicht sehr genau, und zwar 
nicht nur meßtechnich, was die Streuung der Ergebnisse unterstreicht (nach 
Bergstedt®) + 20%), sondern auch, weil man nie ganz sicher ist, ob der 
Einfluß von Strahlungsdämpfung und Stoßverbreiterung ganz eliminiert 
werden kann. Entsprechendes gilt für die Elektronentemperatur. Das Elek- 
tronenbremskontinuum ist nicht die einzige Ursache für die Kontinuums- 
strahlung, da neben Molekül- und Rekombinationskontinuum auch in die 
Bogensäule eintretende Fliissigkeitstrépfchen beträchtliche Störungen hervor- 
rufen. Außerdem wird zur Bestimmung der Elektronentemperatur die Saha- 
Formel, deren Anwendung für Bogengase ohnehin Schwierigkeiten bereitet, 
herangezogen und damit thermisches Gleichgewicht voraussetzt, was erst 
bewiesen werden soll. 


Anzeichen für Nichterfüllung des thermischen Gleichgewichts 


Es gibt eine Reihe von Beispielen, die zeigen, daß der Zustand zumindest 
in der Niederstrombogensäule weiter vom thermischen Gleichgewichtszustand 
entfernt ist, als man bisher annahm. In einer Arbeit von Egorow, Koles- 
nikow und Sobolew’) zeigen die Verfasser, daß ein Bogen, der zwischen 
Kohleelektroden in Edelgasatmosphäre brennt, erheblich vom thermischen 
Gleichgewicht abweicht. Während sich für die Anregungstemperatur nach der 
Boltzmann-Formel ein Wert von 11000 °K ergibt, soll die Elektronen- 
temperatur nur 3500 °K betragen. Wie die Verfasser betonen, dürften hierfür 
die geringen Wirkungsquerschnitte für Stöße zweiter Art und die Existenz 
von langlebigen metastabilen Zuständen bei Edelgasen verantwortlich sein, 
= eine stufenförmige Anregung mit relativ langsamen Elektronen ermög- 
ichen. 


®) R. Mannkopf, Z. Physik 76, 396 (1932). 

‘) H. Witte, Z. Physik 88, 415 (1934). 

5) F.Burhorn, Z. Physik 140, 440 (1955). 

*) K. Bergstedt, Z. Physik 155, 23 (1959). 

) W.N. Egorow, W. N. Kolesnikow u. N. N. Sobolew, Dokl. Akad. Nauk SSSR 
Moskva 121, 440 (1958). 
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Da bei Metallen ebenfalls metastabile Zustände auftreten und ihre Wir. 
kungsquerschnitte im Vergleich zu den im Kohlebogen vorherrschende 
M olekülgasen geringer sind, ist hier ähnliches zu erwarten. 

Bestimmt man daher relative Übergangswahrscheinlichkeiten aus den 
Intensitätsverhältnis von Linien, die die Bogensäule aussendet, so darf ma 
sich nicht wundern, daß die erhaltenen Werte erheblich abweichen von solchen, 
die, soweit möglich, quantenmechanisch berechnet oder durch grundsätzlic 
andere Meßverfahren ermittelt wurden. Vergleicht man z. B. die relative 
Ubergangswahrscheinlichkeiten für 33 Fel-Linien, die sowohl von King 
und King®) als auch von Sobolew’) bestimmt wurden, so zeigen sich beträcht. 
liche Unterschiede. King und King errechnen Oszillatorenstarken aus de 
Gesamtabsorption durch Eisendampf, der sich in einem geheizten K.ohlerohr 
befindet, während Sobolew Übergangswahrscheinlichkeiten aus Linie 
intensitäten eines Bogens bekannter Temperatur ableitet. Im Kingscha 
Kohleofen ist trotz geringerer Stoßzahl der Zustand sicher näher dem ther. 
mischen Gleichgewicht als im Lichtbogen, da hier nur die Besetzungszahle 
der unteren Terme eines Überganges der Boltzmann-Verteilung gehorchen 
müssen und die nahen geheizten Wände das Strahlungsgleichgewicht begin. 
stigen. Bei Sobolew wird eine Boltzmann-Verteilung für die oberen 
Terme eines Überganges vorausgesetzt. Die Diskrepanzen sind so groß, dab 
sie nicht mehr durch Meßfehler bedingt sein können. Setzt man nämlich die 
relativen gf-Werte nach King und Sobolew für A = 2983 A gleich (gfggg; =36), 
so ergibt sich g fg5g5 = 1800 nach King und gfasg; = 66 nach Sobolew. Die 
Werte unterscheiden sich um den Faktor 27,5. Da die Linienintensitat pro- 
portional dem Produkt aus der Besetzungszahl N, und den gf-Werten ist, 
läßt sich diese Erscheinung auch so deuten, daß der obere Term 2° G9 der 
Linie A = 3565 A im Vergleich zum oberen Term y® D$ der Linie 2 = 29834 
zu schwach besetzt ist (kleinere Lebensdauer). Selbstabsorption sollte ausge- 
schlossen sein, da der Bogen zwischen Cu-Elektroden mit einem Zusatz von 
0,3%, Fe brannte. 

Auffälligerweise ist die Übereinstimmung innerhalb eines jeden Multipletts 
recht gut. Es ist auch möglich, darauf weist Ornstein?) bereits hin, daß die 
Ubergangswahrscheinlichkeiten durch die Felder der Ionen und Elektronen 
im Bogen beeinflußt werden, ihre Werte daher von der Temperatur abhängg 
werden. Das trifft besonders für Interkombinationslinien zu. 

Auch ohne direkte Benutzung von möglicherweise fehlerhaften Werte 
der Übergangswahrscheinlichkeiten lassen sich Temperaturmessungen nach 
einer Vergleichsmethode vornehmen. Dazu ist notwendig, das Intensitätsver- 
hältnis zweier Linien mit bekannten oberen Termwerten zu messen. Man er- 
hält eine Gleichung von der Form 


Ist die Bogentemperatur bekannt, so läßt sich der Wert von K angeben, de 
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für das betreffende Linienpaar eine Konstante ist. Durch Bestimmung def —— 


Intensitätsverhältnisses in einem beliebigen Bogen läßt sich dessen Tem} 


8) R.B. King u. A. S. King, Astr. Journ. 87, 24 (1938). 
9) N.N.Sobolew, Sh. exp. teor. Fis. 18, 131 (1943). 
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peratur dann berechnen. Ausgehend von einer angenommenen Temperatur 
von 5020 °K für einen 1 A-Bogen erhielten Gordon, Cotter und Parker?®) 
für Stromstärken von 2,5 A 


und 4 A folgende Temperaturen Tabelle 1 
(°K) (Tabelle 1). Linienpaar | 1A | 2,54 4A 
Die Differenzen liegen nach 
Angaben der Verfasser weit 3402 er 
der Meßgenauigkeit, Fe 3415 
auch die Selbstabsorption ist in 3425 
den meisten Fallen vernach- Fe 2415 5020 5190 5310 
lässigbar klein. Für die Er- 
scheinung selbst wurde keine Fe 3028,7 5020 5300 5460 
Erklärung angegeben. 3998 
3292,0 
Für Serienmessungen von Fe 016,2 4230 | 4068 39101) 
Bogentemperaturen ist wegen 3236,2 
seines verhältnismäßig geringen 3962 
experimentellen Aufwandes und Sn 5034 5020 6580 94501) 


seiner guten Reproduzierbar- 

keit das Verfahren am besten geeignet, das von der Messung des Inten- 
sitätsverhältnisses von Linien mit bekannten Ubergangswahrscheinlichkeiten 
ausgeht, bzw. die im vorigen Abschnitt beschriebene Vergleichsmethode. Dabei 
ist die Frage zu klären, ob dieses Verfahren für den Eisenbogen überhaupt 
angewendet werden darf. Ausgehend von der Gleichung 


9o 


erhält man nach Umformung die Zahlenwertgleichung 


0,625 - AR 
7 


T= 


Hierbei ist 7’ die Temperatur in °K, AE die Energiedifferenz der oberen Term- 
werte E, je zweier Linien in em!, P = 4,,9;v = fx; I die gemessene 
Intensität und lg der dekadische Logarithmus. Trägt man für verschiedene 


Linien die Werte le als Funktion von E, auf, so liegen die Punkte (im 


Idealfall) auf einer Geraden mit der Steigung — 7: 


Für die Vergleichsmethode wird das Intensitätsverhältnis a derselben 


0 
Linie im Bogen mit der unbekannten Temperatur 7, und im Vergleichs- 


®) G.M. Gordon, O.A.Cotter u. D.R.Parker, The welding Journal (Res. 
suppl.) 35, 109 (1956). 

book Diese Temperaturen sind unmittelbar aus den gf-Werten von King und King 
abgeleitet. 

2) Sn war als Verunreinigung im Cu-Mantel der Eisenelektrode. 
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bogen mit der Temperatur 7, ermittelt. Für verschiedene Linien erhält man 

mit Ig?2 =/(E,), wobei die Differenzen der graphischen Darstellung zu ent. 
0 

nehmen sind, 


AE 
T= 
AE 
Die Steigung der resultierenden Geraden hat die Größe En _ ER 


0 
Für die Versuche wurden Elektroden verschiedener Abmessungen au 
technischem Eisen benutzt, die Stromstärken von 1,5 bis 200 A variiert. Die 
Intensitätsmessung erfolgte unter Verwendung von Quarz- und Glasspektre 
graphen hoher Auflösung auf photographischem Wege. Die Eichung der spek. 
tralen Empfindlichkeit erfolgte unter Verwendung der Strahlung des positiven 
Kohlebogenkraters. 


Um den Einfluß der Selbstabsorption wenigstens größenordnungsmäßig 
bestimmen zu können, wurde der Bogen mit seinem eigenen Licht durch 
strahlt. Der Bogen befand sich im Krümmungsmittelpunkt eines Hohl. 
spiegels, so daß die zentralen Teile in sich abgebildet wurden. Durch eine 
schmale Blende unmittelbar hinter dem Bogen, d.h. zwischen Bogen und 
Hohlspiegel, konnte erreicht werden, daß der Spektrographenspalt über eine 
Abbildungslinse teils durch den Bogen, teils durch das Licht von Bogen und 
seinem Hohlspiegelbild beleuchtet wurde. Die Spektrallinien waren daher 
quer zur Dispersionsrichtung in ihrer Helligkeit abgestuft. Aus diesem Inten- 
sitätsverhältnis läßt sich die Größe der Selbstabsorption abschätzen. 


Von den etwa 250 Fel-Linien mit bekannten relativen Übergangswahr- 
scheinlichkeiten (aus den gf-Werten von King und King®), Carter!) und 
Aarts") abgeleitet) wurden zunächst die Linien ausgesondert, die durch eng 
benachbarte Linien auch solche von Legierungsbestandteilen des Eisens, 
gestört wurden, ferner diejenigen, die zur Selbstumkehr neigen. Mit den ver- 
bleibenden Linien wurde versucht, die Temperatur zu bestimmen, indem durch 
graphische Mittelung die Temperaturgerade gezogen wurde, da die Einzelwerte 

sehr stark streuen. Für den mittleren Teil der 


Tabelle 2 Bogensäule ergaben sich in Abhängigkeit von det 

A °K Stromstärke folgende Temperaturen (Tab. 2). 
Die Mittelwertbildung unterliegt einer gewissen 
2 4100 Willkür, so daß die Temperaturen um 300° diffe. 
5 4700 rieren können. Die Achstemperatur, die sich durch 
2 == Umrechnung von seitlicher Intensität auf radialen 
100 5600 Emissionskoeffizienten ermitteln läßt, ist je nach 
200, 5700 Breite der Bogensäule bei 2 A etwa 500° bei 1004 


etwa 200° höher als die „‚mittlere‘‘ Temperatur. 
Die so erhaltenen Temperaturen erscheinen jedoch im Vergleich zu Werten, 
wie sie in der Literatur zu finden sind, zu niedrig. Die Abweichungen sind um 
so größer, je niedriger die Stromstärke ist. Auch die von Sobolew®) angege- 
benen Übergangswahrscheinlichkeiten führen nicht zu besseren Ergebnissen. 


13) W.W.Carter, Physic. Rev. 76, 962 (1949). 
14) J. Aarts, D. Hartung u. €. J. Bakker, Physica 20, 1250 (1954). 
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Da keine leicht ionisierbaren Bestandteile im Bogengas nachzuweisen waren, 
sind diese Werte auch hinsichtlich des thermischen Mechanismus bedenklich. 

Um zu ,,richtigen‘‘ Temperaturen zu kommen, kann man sich auf geeignete 
Linien beschränken. Dabei zeigt sich, daß Linien, deren obere Terme zum 
Triplettsystem gehören, im allgemeinen zu niedrige Intensitäten aufwiesen. 
Auch solche Linien sind zu schwach, deren oberer Term durch andere Über- 
ginge mit großer Übergangswahrscheinlichkeit nur geringe Lebensdauer hat. 
Es ist jedoch sehr schwer, allgemein gültige Regeln für diese Abweichungen 
zu finden. Linien desselben Multipletts zeigen dagegen, abgesehen von even- 
tueller Selbstabsorption, immer ein gleichartiges Verhalten. 


Man könnte zunächst annehmen, daß die Schwierigkeiten sich verringern, 
wenn man aus dem Spektrum eines 2 A-Bogens solche relativen Übergangs- 
wahrscheinlichkeiten festsetzt, daß alle Punkte eine Gerade von gewünschter 
Neigung ergeben, doch zeigt sich, daß für höhere Stromstärken wieder starke 
Streuungen auftreten, die mit wachsender Stromstärkendifferenz zunehmen. 

Bei Stromstärken über 100 A nehmen die Streuungen der Meßpunkte eine 
andere Form an, sie dürften hier weniger durch die Abweichung vom thermi- 
schen Gleichgewicht verursacht werden, sondern durch die zunehmende Selbst- 
absorption. Auch die Intensitätsmessung wird durch Linienverbreiterung 
und kontinuierlichen Untergrund erschwert. 

Unter Verwendung der gf-Werte aus (8) und (13) lassen sich mit Hilfe der 
Linien 4227, 4247, 4375, 4427 Ä einigermaßen brauchbare Temperaturwerte 
herleiten, die bei Stromstärken über 30 A dem 


wahren Wert näher kommen dürften (Tab. 3). Tabelle 3 
Die Achstemperaturen sind um 500—200° höher. A oK 
Die reine MeBgenauigkeit läßt eine Temperatur- 
angabe auf weniger als + 150° zu. Daher ist diese 2,5 | 4800 
Methode zur Bestimmung von Temperaturdifferenzen 5 5050 
besser als jede andere’ geeignet. Echte Temperatur- = 
schwankungen (+ 200°) können bei kurzer Be- 100 6000 
lichtungszeit nachgewiesen werden. Vor allem lassen 200 6100 


sich bequem Untersuchungen über die Abhängigkeit 
der Temperatur von Bogenlänge, Elektrodendurchmesser usw. durchführen. 
Meßtechnisch wäre die bereits genannte Vergleichsmethode vorzuziehen, 
da dann die spektrale Empfindlichkeit nicht berücksichtigt zu werden braucht. 
Die Genauigkeit verringert sich jedoch dadurch, daß die Temperatur des mit 
nicht zu kleiner Stromstärke betriebenen Vergleichslichtbogens selbst Schwan- 
kungen unterliegt. 
Da zur Vergleichsmethode Übergangswahrscheinlichkeiten nicht benutzt 
werden, sind die Beweise für die Nichteinhaltung des thermischen Gleich- 
gewichts noch eindrucksvoller. Es sind immer die gleichen Linien, deren MeB- 


punkte Igge außerhalb der MeBgenauigkeit aus der Reihe tanzen. Anderer- 


seits fügen sich Linien mit beträchtlicher Selbstabsorption noch verhältnis- 
mäßig gut ein. 

Gleiche Erscheinungen traten auf, wenn an Stelle der Fe I-Linien Fe 
II-Linien benutzt werden. Mit den von Sobolew®) angegebenen Ubergangs- 
wahrscheinlichkeiten ließ sich überhaupt keine Temperatur bestimmen, da 
die Punkte im Diagramm derartig verstreut waren, daß sich keine Tem- 
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peraturgerade ziehen ließ. Hinsichtlich der Vergleichsmethode zeigten sich 
ähnliche Abweichungen wie bei Fe I-Linien, obwohl wegen der geringen Kon. 
zentration von Fe*+-Ionen die Selbstabsorption sicherlich kleiner ist. 

Möglicherweise sind die metastabilen Zustände a?F und a®F stärker be. 
setzt, als es nach der Boltzmann-Formel der Fall sein sollte, zumindest in 
der Randzone. Alle stärkeren Linien, die auf diesen Niveaus enden, zeigen 
eine Selbstabsorption bzw. Selbstumkehr von solcher Größe, wie sie bei starken 
Linien, die mit dem Grundniveau a®D kombinieren, auftritt. 


Zusammenfassung 


Da die Bedingung des thermischen Gleichgewichts im Eisenbogen nicht aus- 
reichend erfüllt ist, führt die Temperaturbestimmung mit Hilfe von relativen 
Übergangswahrscheinlichkeiten, die nicht aus Bogenentladungen abgeleitet 
sind, im allgemeinen zu Temperaturen, die im Vergleich zu anderen Meß. 
methoden zu tief liegen. Dabei ist die Abweichung um so größer, je kleiner die 
Stromstärke ist, woraus zu entnehmen ist, daß Annäherung an thermisches 
Gleichgewicht erst bei großen Stromstärken zu erreichen ist. 

Da der Grad der Annäherung außerdem von der jeweiligen Zusammen- 
setzung der Bogengase abhängig ist, zeigen aus Bogenentladungen ermittelte 
Übergangswahrscheinlichkeiten unterschiedliche Werte, je nachdem, ob das 
zu untersuchende Element selbst das Elektrodenmaterial bildet oder als 
Zusatz zu anderen Materialien (Kohle, Kupfer usw.) dient. Auch die um- 
gebende Atmosphäre (Luft, Edelgas) dürfte von Einfluß sein. In dieser 
Richtung sind weitere Untersuchungen im Gange. 

Diese Arbeit ist zum Teil in enger Zusammenarbeit mit dem Zentral- 
institut für Schweißtechnik in Halle entstanden. Herrn Dipl.-Ing. Wein- 
schenk und Herrn Ing. Hirschfeld sei an dieser Stelle für ihre Unterstützung 
bei der Durchführung von Messungen, wie auch für ihre wertvollen Diskussio- 
nen gedankt. Desgleichen habe ich Herrn Prof. Dr. Schallreuter für sein 
förderndes Interesse zu danken. 


Greifswald, Physikalisches Institut der Ernst-Moritz-Arndt-Universität. 


Bei der Redaktion eingegangen am 13. November 1959. 
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The Role of Initial Stresses in Lattice Dynamics 
By E. 8. Rajagopal 


Summary 


Considerations of the initial stress are introduced in the general lattice 
dynamics in discussing the symmetry of the elastic constants. It is proved 
that if one describes the infinitesimal elastic behaviour from the state S, 
the symmetry of the behaviour follows that of the state S but is indepen- 
dent of the magnitude of the initial stresses frozen in S. The bearing of the 
result on the current formulations of lattice elasticity is indicated. 


$ 1. Introduction 


In an infinite crystalline lattice, the condition that each particle is in equili- 
brium, leaves the initial (uniform) stresses indeterminate. Under the presump- 
tion that the classical theory is not applicable when “frozen” stresses are 
present and that it requires more than 21 constants to describe the elastic 
behaviour of such bodies, Huang?) has discussed the conditions for the vanis- 
hing of these initial stresses by a comparison of his equations for the wave 
propagation in the (initial) stressed state and in the (natural) stressfree state. 
Since these conditions have been reproduced later also*)%) it is necessary to 
point out that the considerations of Huang are incomplete‘) and that in the 
classical infinitesimal theory the elastic behaviour of an infinite lattice is 
independent of the amount of the initial stresses, provided the behaviour is 
described consistently in a coordinate system imbedded in the initial (stressed) 
state. 


§ 2. Initial stresses in infinitesimal Deformations 


It is first necessary to point out that almost all the work done in lattice 
elasticity is valid only for infinitesimal deformations. Murnagan®) proved 
in 1937 that the indiscriminate use of the derivative of the Helmholtz free 
energy to obtain the stress components is valid only for infinitesimal defor- 
mations and Birch®) pointed out this defect in the work of Born and others. 


1) K. Huang, Proc. Roy. Soc. London 208 A, 178 (1950). 
*) M. Born and K. Huang, Dynamical Theory of Crystal Lattices, Oxford 1954) 
p. 237. 
; 3) G.Leibfried, Handbuch der Physik, Ed. Flugge, Springer, Berlin 1955, VII. 
198. 
*) Regarding the incorrectness of Huang’s remarks on the role of initial forces in 
the validity of Cauchy’s relations, see Stakgold, I. Quart. Appl. Math. 8, 176 (1950). 
5) See his Finite Deformations of an Elastic Solid — Wiley, New York 1951, p. 56. 
*) F. Birch, Physic. Rev. 71, 809 (1947). 
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This fact was again emphasized by Truesdell’). In spite of this, one finds 
Born and Huang’) using the derivatives of the free energy and still claiming 
the results to be true for finite deformations. The same charge appears to 
hold against the work of Lowdin?) also. 

It is not out of place to emphasize the approximations involved in the 
infinitesimal deformations!). The usual Lagrangian description (material des. 
cription) of the strain is = 4 (us,g+ ug,a+ Uy, « Uy,s) and the Eulerian 


description (spatial description) is &,5 = + — N Uy, « Uy,a). If 
Y 


one neglects the powers of the displacement gradients greater than 1, and 
writes &,5 =% (Us, + Ug,.) one gets the classical infinitesimal elasti- 
city. Here no distinction need to be made between the material and spatial 
co-ordinates. If one uses the full strain tensor without this approximation, 
one gets the geometric non-linearity. In both these cases the stresses are 
lienar functions of the strain (Hooke-Cauchy law). If one considers the si- 
tuation where the stresses contain higher powers of the strains, one gets 
physical non-linearity, as for example in the third order elastieity of 
Murnagan®). In the last two cases, the deformations are finite, the mathema- 
tical equations are non-linear and the treatment is quite complex!!). In all 
the cases the strains are numerically small if plasticity is not to be reached. 


The subject of initial stresses is important in technical applications and in 
geophysics. A very clear discussion of the problem of initial stress, its relevance 
to finite deformations but not to infinitesimal elasticity, is given in the review 
of Truesdell!2). It is very unfortunate that in spite of this work, uncritical 
remarks on the role of initial stresses have been made in several works on 
lattice dynamics. The present situation is a simplifications of the general 
theory. (1) In an infinite atomic lattice, the initial stresses are uniform every- 
where and the three dimensional periodicity is maintained. The lattice remains 
homogeneous and Hookean. (2) The second main difference is that in the 
general theory of initial stresses the situation is such that the frozen stresses 
take the body from its “natural” (stressfree) states to the initial (stressed) 
state, the applied strains takes the body to the deformed state and one is re- 
quired to connect the deformed state to the natural (stressfree) state. In the 
present case of infinite crystall lattices, the hypothetical natural state is 
unknown while every property of the initial state can be computed. One is 
required to connect the deformed state to the initial state, which alone has 
observational reality. Such a situation may arise, for example, from the zero 
point and thermal energies. A full discussion of this point is contained in 
the report “The validity of the recent theories of crystal elasticity” by R.S. 
Krishnan and E.S. Rajagopal (Proc. Ist Symp. Solid State Physics, 
Bangalore 1960) (referred to as Part I). A review of the subject, to be pub- 
lished soon (Part V), also contains a full discussion. 


7) C. Truesdell, J. Rational Mech. Math. 2, 602 (1953). 

8) M. Born and K. Huang (I. c.) — Chap. VI. 

®) P.O. Léwdin, Ark. Mat. Astron. Fys., 35, No. 30, p. 11 (1948). 

1) V. V. Novozilov, Foundations of the Non-linear Theory of Elasticity, Graylock, 
Rochester 1953, p. 128, § 4. 

11) C, Truesdell, Z. Angew. Math. Mech. 36, 97 (1956). 

12) C. Truesdell, J. Rational Mech. Math. 1, 173, 217 (1952); 2, 605 (1953). 
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First consider a solid of finite extent, which is in equilibrium in its initial 
state, but still has frozen stresses 12g = 1, 2,3) and body forces 7° 
(per unit mass). The balance of the forces and the moments acting on a volume 
element in this state leads to the usual relations 


at, 
Y 
at the points in the interior of the body and to the relations 
(2, cos (By) = #8, (2) 


at the surfacees of the body. Let now infinitesimal deformations be applied 
to this body. The usual convention? (it is only a convention) is to refer the 
strains to the state of the body before the deformation (material description) 
and the stresses to the state after the deformation (spatial description). But 
if the displacement gradients are infinitesimal, the two descriptions are iden- 
tical. Consequently the equilibrium conditions for the body will be 


oT, + > 5 =F,; tas—tsa=9; 


where 

and the coordinate system x, is still referred to the initial state containing 
the frozen strains. In view of eqs. (1) and (2), these become 


at, 
oT a+ =Fai ten =0 (3) 


at the interior points, and 
lap cos = lay (4) 


at the boundaries. Thus the infinitesimal elastic behaviour can be discussed as 
though the initial stresses were absent!#). The result is actually a trivial case 
of the statement for finite deformations. In that case, one must consider the 
changes in 9, v, 7%, 2 etc., caused by the applied deformations and the problem 
is necessarily complicated. Biot!5) has proved the independence from the 
initial stresses if one consistently uses the Lagrangian description. 
The equations of vibrations are obtained by leting F, denote the d’Alem- 
bertian force, and in the absence of body forces!*), one gets 
eu, 
= (5) 
These equations must be supplemented by prescribing the boundary conditions. 


u? I. S. Sikolnikoff, Mathematical Theory of Elasticity, McGraw Hill, New York 
, p. 29, 64. 

M4) P. Appel, Traite de Mecanique Rationelle, Gauthier-Villars, Paris 1928, III, p. 
ot S. Timoshenko and J. N. Goodier, Theory of Elasticity, McGraw Hill, New York 

1, p. 238. 

35) M. A. Biot, Philos. Mag. 27, 486 (1939). 

16) In the atomic case, the initial state of the body has no force acting on the atoms. 
In the genral case, the presence of body forces leads to forced vibrations and does not alter 
the situation. See A. E. H. Love, Mathematical Theory of Elasticity-Dover, New York 
1944, p. 304. Truesdell) has pointed out that Love’s treatment of the initial stress 
problem leaves much to be desired. 
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All these statements hold good for infinite media also, except that there 
are no boundary conditions. Since the impressed stresses and the strains in 
the Hookean body are connected by 


tag = Capyö&yd (6) 
yö 


the equations of wave propagation are 


eu = Capy dlp ls Uy (l; = direction cosines). (7) 


Thus in the general case the effective elastic constants for wave propagation 
in the infinite media have symmetries independent of the initial stresses. The 
only assumptions involved are that the body is Hookean and homogeneous 
and that the additional deformations are infinitesimal, which, as discussed 
earlier, are satisfied in the infinite atomic lattice under consideration. 

Though c,gy9 have the symmetry in «ßyö independent of the magnitude 
of the initial stresses, the number of independent constants and their numerical 
values depend on the symmetry and nature of the initial state of the body. 
There can be at most 21 constants in the general case. (In the present 
atomic case, the infinite lattice can at most become triclinic. But it can never 
become heterogeneous since the three dimensional periodicity is maintained). 
The constants will be the same as those of the hypothetical natural state only 
if the initial stresses are infinitesimal; for then, one deals with the sum of two 
infinitesimal deformations which is still an infinitesimal affine transformation. 


§ 3. Some comments on the Work of Huang 


The argument of Huang!) questioning the statement of Brillouin”) 
that the symmetry of the elastic constants is not affected by the initial stresses 
is based on the conclusion that if one subjects the medium (with the frozen 
stresses) to a homogeneous deformation and applies the rigid rotational in- 
variance of the strain energy, one gets 


+ Caypa —Cyapa = 0. (8) 


Hence ¢,yg4 4 Cyaga in geneal, if the initial state itself is strained. But 
this argument, though perhaps necessary, is not sufficient by itself (It is just 
like building up a system of mechanics with the only condition of rigid rotatio- 
nal invariance). Its compatibility with the fundamental laws of mechanics is 
not at all considered and as the theory stands, it is not even applicable to 
ideal crystals at 0 °K, since the stresses due to the zero point energy still 
persist. Infact, the very idea of a solid state presupposes the forces acting 
across surface elements from one part of the body into the other. This gives 
the background level of stresses present in any concept of a solid (cohesion). 
(See also Part I and V.) 

The real source of error in Huang’s work lies in his eq. (6.4) which in the 
present notation is 


= Exp + Wap + = Ony (9) 


be L. Brillouin, Tenseurs en Mecanique et en Elasticite, Masson, Paris 1949, 
p. 232. 
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The summation term $, is of a higher order infinitesimal than the rest of the 
terms. Huang has used infinitesimal deformations at other places, for example 
his eq. (6.12) (see also $ 1), and should have neglected the $, terms. Instead, 
he retains them to get eq. (8). In fact if one wants to retain the higher order 
terms, one must make a clear distinction between the Lagrangian and Eu- 
lerian descriptions. The resulting wave equation is no longer of his form 
(eq. 6.12) viz., 

o Vu, = 2, (10) 


The exact wave equation in the case of initially stressed bodies has been given 
in Lagrangian description by Biot!®) and in Eulerian description by 
Sakadi and Takizawa’). It involves the initial stresses also and even 
in the absence of initial stresses the equation contains terms of higher orders 
in addition to the terms of eq. (7). If a wave equation of the type eqs. 7. or 10 
is to be obtained at all, one must necessarily use infinitesimal deformations 
and in such a case the Superposition Theorem is valid, as indicated in § 2. 

Hence one must conclude that Huang’s discussion of the initial stresses in 
a crystal lattice is wrong. An important consequence of this conclusion is 
that the method used by Huang to introduce the necessary symmetry into 
the expressions of the elastic constants deduced by Begbie and Born?) 
is not valid®!). The real difficulty in the expressions for the elastic constants 
lies elsewhere, as will be discussed in Parts III and IV. 


Encouraging help from Prof. R. 8. Krishnan and Dr. 8. K. Laxmana 
Rao and the award of a Senior Research Fellowship by the C. S.I. R. are 
gratefully acknowledged. 


18) M. A. Biot, J. Appl. Phys. 11, 522 (1940). 

1) Z. Sakadiand E. Takizawa, Mem. Fac. Engg. Nagoya 1, 95 (1949); 2, 51 (1950). 

2) G.H. Begbie and M. Born, Proc. Roay. Soc. London 188 A, 179 (1947). 

21) Dr. G. Leibfried (private communication) has suggested that the necessary 
symmetry can be introduced into the expressions of Begbie and Born, if one proves, 
by using the rotational invariance condition, the relations given by him (Leibfried, 
l.e., p. 148). We shall exhibit the difficulty with this procedure also in Part. III. 


Bangalore-12 (India), Department of Physics, Indian Institute of Science. 


Bei der Redaktion eingegangen am 27. November 1959. 
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The Rotational Invariance Condition in Lattice Dynamics 
By E.8. Rajagopal 
With 1 Figure 


Summary 


The note discusses the validity of the rotational invariance condition 
used in lattice dynamics. It is pointed out that rotations of non-central 
fields cannot be represented by displacements alone. Formulating the equa- 
tions of lattice dynamics referred to Lagrangean co-ordinates imbedded in 
the crystal it is shown that the usual formulae do not strictly represent 
the physical requirement of rotational invariance. The result has consider- 
able effect on several applications of the lattice theory. 


§ 1. Introduction 


In order to develop the basic equations of the general lattice dynamics, 
where no assumptions regarding the interatomic potential functions are made, 
it is necessary to employ the various invariance relations obeyed by the po- 
tential energy of the lattice. It is well-known that in lattice dynamics, the 
crystal can be considered to be made up of atomic nuclei whose co-ordinates 
can be represented in terms of the cell indices m, p,q, ... and the basis indices 
j, k,l... The rotational invariance condition!)?) can then be represented as 


The aim of the present note is to show that eq. (1) has no basis in the general 
lattice dynamics, since the mathematical formulation does not strictly represent 
the pyhsical condition of the invariance of the potential energy of the lattice 
under a rigid body rotation. 

Consider, as an example, a calculation of the properties of an arbitrary 
(triclinic) lattice restricting oneself to the first neighbour interactions. To 
an atom (m,), let the first neighbour, which is one in number, be (px). Re- 


membering that D,, wr ? ) is the force in the « direction on the (m,) atom 
due to the (p;) atom for unit relative displacement of the atoms in the ß diree- 
tion, the force constant scheme ®, , fas P ) (in terms of x, ß for fixed ye P )) 


jk 
is symmetric (see eq. 9), but otherwise completely arbitrary, since there is 


1) G. H. Begbie and M. Born, Proc. Roy. Soc. London 188A, 179 (1947). 


2) M. Born and K. Huang, Dynamical Theory of Crystal Lattices, Oxford 1954, 
p. 222. 
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no reduction by the crystal symmetry. But using eq. (1), one gets ®,,&, = 
@,,&3. Writing the relations in full and using eq. (9), one finds that eq. (1) 
gives only one independent force constant and that the force field becomes 


a central (radial) force scheme. 


The genesis of the strange results is easily found. The relation (1) is ob- 
tained by subjecting the nuclei to the displacements 


=X 


(2) 


and equating the resultant change in potential energy to zero. But the displa- 
cements (2) do not strictly constitute a rigid body rotation for an arbitrary 


lattice. To see this, first consider a system of 
point masses with the initial configuration 
shown by full circles in Fig. 1(a) and let it be 
subjected to displacements of eq. (2) with 


000 
(brocken circles) in Fig. 1(a) is seen to be a rigid 
rotation of the system. The same is true for a 
continuous distribution of matter’). Now con- 
sider a system of oriented particles (spin, 
magnetic moment, quadrupole moment, etc., 
are simple quantum mechanical examples of 
directive properties of the nuclei) with the 
initial configuration (full circles) as in Fig. 1(b). 


0650 
Wap -(~ 0 ) The resultant configuration 


The result of applying the displacements u, 6) 


(eq. 2) is also shown and the resulting con- 
figuration has certainly a difference from the 
rigid body rotation Fig. 1(c). It is to be remem- 
bered that this is a striet adherence to the ab 
initio principles and that there is no question 


3) In fact, eq. (2) is derived from the continuum 
concept of affine transformation by replacing thevolume 
elements by their mass centres to form a system of dis- 
crete elements. Now the anisotropy of matter, which 
is visualised in the atomic theory by the atomic lattice 
configurations, can be obtained un the continuum 
theory only if “every particle of the continuum 
possesses distinct characteristic properties 
along distinct directions emanating from it 
and therefore the specification of their orien- 
tation isessentially included in the description 
of the state of the continuum”’. (Italics those of 
E.Hellinger, Encyk. Math. Wiss. IV—4, 609 (1914). 
It is obvious that in such a case the replacement of 
the volume element by their mass centers is not per- 
missible, since a displacement always preserves the 
orientation of the mass point whereas the affine trans- 
formation may involve a rotation of the volume 
element. 


\ 


Figure 1. Effect of applyin, 


the Displacements u, = 


a) a system of 

point masses, b) a system of 

dipoles (oriented particles). 

The rigid body rotation is given 

in c) and is to be compared 
with b) 
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of any external fields involved. In fact, Corbon and Stehle®) have given a 
transformation which does represent rigid rotation for non-central fields. It 
will be seen that it cannot be represented in terms of displacements alone 
and this is the motivation for the artifice introduced in § 2. For a calculation 
of the energy change caused by eq. (2), see Part V. 


At this juncture, it is necessary to digress and to point out the connection 
with the quantum mechanical adiabatic approximation. In the original 
treatment of Born and Oppenheimer®), the nuclear spin and magnetic 
moments were not considered. But with the advent of h. f. s*) and more re. 
cently NMR’) studies, it was recognised that the ground state Hamiltonian 
must include the potential energy, not only due to the nuclear configurations 


Ex (*) but also due to the nuclear parameters like spin and magnetic moment. 


(Indeed, the nuclear dipole and quadrupole moments are being invoked in 
lattice dynamics by several workers’). The representation of the rigid body 
rotation must then include besides the variations of the nuclear positions the 
changes in the nuclear parameters (internal co-ordinates), as for example 
given by Corbon and Stehle‘*). The present problem is intimately con- 
nected to the “slip effect” in molecular spectra, i. e., the fact that the inner 
spherical shells of atoms point to fixed spatial orientation even when the mole- 
cule is rotating about itself®). The energy changes involved in such nuclear- 
nuclear interactions are in the microwave region and small compared to 
electronic transitions. But in a theory of perfect lattices where the electronic 
interactions are completely removed and only nuclear interactions are con- 
sidered, such types of energies cannot be neglected, especially when the theory 
is to be applied to acoustic waves of practically zero frequency. The matter is 
quite serious since the energy density of the deformation of eq. (2) can be 
comparable to the strain energy density (see Part V). 


Before giving a rigorous proof of the invalidity of eq. (1) directly, it is 
pertinent to comment that Laval?) has explicitly rejected the condition (1) 
as being inadmissible in the general theory. But he has not offered a proof of 
this, nor has he noticed the error in eq. (2). Earlier, Harrison (quoted by 
Launay!!) has noticed that the general atomic force constants do not satisfy 
eq. (1). But instead of investigating the anomalous situation why any force 
field should vary under the displacements (2) when all physical considerations 
demand that there be no variation for a rigid rotation Launay has commented 
that it is the general force field which must be carefully studied. 


‘) H.C. Corbon and P. Stehle, Classical Mechanics, Wiley, New York 1950, p. 324. 

5) M. Born and J. R. Oppenheimer, Ann. Physik 84, 457 (1927); A. Haug and 
G. Sauerman, Z. Physik 158, 269 (1958). 

6) C. H. Townes and B. P. Dailey, J. chem. Physics 20, 36 (1952). 

7) J.S. Waugh, Annals N. Y. Acad. Sci. 70, 900 (1958). 

8) J. M. Luttinger and L. Tisza, Physic Rev. 70, 954 (1946); 72, 257 (1947). V.S. 
Mashkevich, and K. B. Tolpygo, Soviet Physics JETP, 5, 435, 707 (1957).; M. Lax, 
Physics Rev. Letters 1, 133 (1958); H. Kornfeld, Z. Physik 22, 21 (1924). 

®) C.H. Townes and A. L. Schawlow, Microwave Spectroscopy, MrGraw Hill, 
New York 1955, p. 17. 

10) J. Laval, J. Physique Radium 18, 247 (1957). 

11) J. Launay, Solid State Physics 2, 248 (1956). 
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Both Launay and Laval have realised that eq. (1) can be geometrically 
satisfied, when there are special symmetry elements at each atom. But in 
the general case (class C,) which will be considered here, no such symmetry 
elements can be assumed. 


§ 2. Imbedded co-ordinates in lattice Dynamics 


A consistent formulation of lattice dynamics can be made by introducing 
a set of imbedded co-ordinates. One refers all the mechanical operations per- 
formed on the lattice to a system of cartesian co-ordinates €, unalterably 
fixed in the body of the crystal. (It is surprising why this idea, so natural in 
hydrodynamics and elasticity has not been used in lattice dynamics.) When a 
rigid translation or rotation is applied to the crystal the co-ordinate frame £, 
moves along with the crystal so that the quantities referred to these axes will 
not be altered. In other words, the physical condition of translational and 
rotational invariance is automatically incorporated in this system. Since the 
object is to test the validity of eq. (1), the plan of the proof will be to for- 
mulate the equations governing some phenomenon, say the propagation of 
long acoustic waves, and then to compare these exact equantions with the 
equations given by the usual formulation with and without using eq. (1). 
This comparison will show whether the condition (1) is equivalent to the 
invariance under rigid body rotation or not. 

Apart from the usual lattice hypothesis, the only condition that need be 
assumed is that the lattice is in equilibrium initially and after an infinitesimal 
deformation. For the infinite lattice, this condition leaves the initial uniform 
stresses indeterminate (e. g., due to zero point and thermal energies). But it 
is well-known that the infinitesimal elastic properties do not depend upon 
the magnitude of these frozen stresses (Rajagopal, Part. II). The present 
formulation will reveal itself to be independent of the initial stresses, quite in 
accordance with the continuum theory. Regarding the interatomic forces, 
many body forces will not be considered and so the general notation of Born 
and Huang?) (except for minor obvious changes) can be preserved. 

The potential energy ® of the lattice considered as made up of two body 
interations can be represented as (see also Part IV) 


j 


where the derivatives evaluated at 7 = 0 can be represented by ®, ve 


®,5 ry ), «++, Since the body is in equilibrium with no net force on each 


®,("*) == @, (4) 


The atom (p,) having a displacement 9 ( 2) exerts on the atom (m,) with a 


displacement 7 ( ry) the cohesive force f = given by 


& 


| 
| 

©) 
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The total force exerted on (m,) by all the other atoms is 


The summation extends over all the atoms except (m,) and to remove this 
restriction, it is sufficient to take 


so a 

(6) 
an 


a] 


The relation (6) can be introduced in other ways also, but they are all 
equivalent. 
In the initial equilibrium state, two atoms (m,) and (p,) exert on each 


other cohesive forces Sr Ge and Sr G 7 which from the well-known law 
of mechanics are equal and opposite. After the infinitesimal deformation, 
Sr’ r = 1) = — Ir 4 u where Sr’ (5 2 the new cohesive force is equal to 


Sr +f So = 4). Using the expression (5), 


one has 
But since 
one has 
Viswanathan?) has deduced eq. (9) by considering the potential between 


ay 
m p 
and hence can be made symmetric in x and ß. Laval?®), besides deducing 
eq. (9), has discussed how this symmetry is demanded by several other factors 


also. Though the fact that (7%) a. 167) has been mentioned by 


the two particles and observing that 


depends only on 


k 
Born), it is surprising that Begbie-Born!) have not used this relation to 
deduce eq. (9). 
Another relation can be easily obtained as 


= £, bet (10a) 


2) K.S. Viswanathan, Proc. Ind. Acad. Sci. 389A, 199 (1954). 
13) M. Born, Atom Theorie des festen Zustandes, Encyk. Math. Wiss., V-3, p. 536 
(1923). 
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This results from the fact that 


Now, in the sum over a, the pairs of quantities with lattice translations + a 
cancel each other. Moreover, using (6), one can rewrite (10a) as 


1) = 0. (10b) 


As mentioned earlier, no relations can be obtained by rigid body move- 
ments since the reference frame is imbedded in the body itself. 


§ 3. Long Acoustic wave Propagation 
The equations of motion of the (m,) atom, using eq. (7) becomes 
(m\ _ m—p p 
(11) 


The plane progressive wave solutions, taking into account the lattice periodi- 
city are 


(7) "Ca (12) 
which yield the 3 g equations (g = No. of particles in the cell) 
ot (i) = 4) (b) (13) 
where 


Because of (8), (9), 
(15) 


The equations for long acoustic waves can be solved by the Born approxi- 
mation procedure?). It is sufficient to indicate briefly the operations. 

Writing ¢ y for y and expanding the various quantities in powers of e and 
finally letting e > 1, one has 


= 08 


and 


(n) =eo, (n) (16b) 
Ge (j|n) = Coa + € lias (j|n) + Co, + (16e) 


ween 
) (16a) 
where 
uc 
(17) 
d by 
on to 
0 
(10a) 
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It will be seen that I>, IT, --- are real and symmetric in ja,kf 
while J, --- are imaginary and antisymmetric in the pair (j «), (k ß). Their 
symmetry is more than that of the corresponding quantities in Born and 
Huang. Since eq. (13) have 3 g values of w and the associated £, (7), these 
solutions are specified by w (n), Cx (j|n), with n = 1, 2,... 3g. 

The perturbation equations are obtained as 


0 = ST, (24) (k|m) (20) 
ik - j - 

0 = ET + en) 


To 
From (6) and (17), one has 


(23) 
and from (10) and (18) 


(22) 


>» jk = 
(Hs Hx) ( 3) 0. (24) 
The zero-order equations (20) always have the solutions 


Coa (k|n) = ur" (n) (25) 
which represent the translations. 
With the zero-order solutions, the first order equations (21) become 


where £1, (k|n) are to be evaluated. The associated homogeneous part is 

identical with (20) and has solutions of the form (25). So by the algebraic 

theorem, the solvability condition of (26) is obtained by multiplying the in- 

homogeneous part by u", summing over j and equating the result to zero, 
that is, 

1 jk: - 

which in virtue of (24) is satisfied. Hence eq. (26) are solvable. One also 


notices that when eqs. (26) are multiplied by u" and summed over j, both 
sides vanish. This means that among the g equations for a fixed value of 4, 
only g — 1 can be independent. Without loss of generality set 


fia (o|n)=0; «=1,2,3 (28) 
and consider the remaining 3g — 3 equations in (39 — 3) unknowns. In 
general, the (3g — 3) x (39 — 3) matrix r,(2)) (j,k =1,2,..., 9-3 


is non-singular and let its inverse be @’ the elements of which satisfy 
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Multiplying eq. (26) by a(i ) and summing over j, x one gets 


Introduce the 3g-3g matrix @, by bordering G’ with zeros as follows: 
jk jk 
if jk+0 
== otherwise 


One also has = 


which itself is a consequence of the symmetry of a2 er With G, eq. (30) 


) which follows from a symmetry of @’ 


becomes 


kl Bu 
where the indices k/ go over 0,1,....g—1. (The summation dummy index 
u should not be confused with the mass of the !th particle u). 
The second order eq. (22) becomes, on using (25) and (31), 


Ik\ 4, 
+ 
where (>, (k|n) are the inknowns. The solvability condition yields 


It is now convenient to consider the terms the right side. 

(u, Tg = 472 &, (3) &5 Yy Ys (34) 


where v is the volume of the cell. From (19), it is seen that 


[~y,B 6] = [By, 6] = [x 6, By] = [B 4,07] (35) 


but no interchange between x and y or between ß and 64 is permissible. 


= (ay, B vy ys- (36) 


(33) 


It is seen that 

(~y, B 6) = (B46, a) (37) 
but an interchange of x and y or of ß and 6 is not permissible (since the rota- 
tional invariance is already included). 
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Introducing 


[a y, 6] + (ay Bd) = {x y, Bd} (38) 


0 = tat {a y, BO} Yy Yons- (39) 


(eq. 33) becomes 


It is appropriate now to test eq. (1). A comparison of the exact equations 
(34—39) with the equations deduced from the usual formulism (Born and 
Huang?) of referring to an external co-ordinate system shows that 

(i) the exact equations agree with such of those results of the 
usual theory which do not make use of eq. (1) or (2), and 

(ii) the use of eq. (1) or (2) introduces additional symmetry 
relations which are not warranted by the exact procedure of trea- 
ting the physical conditions of translational and rotational in 
variance of the potential energy. 

Hence in the general case, eq. (1) is not strictly equivalent to the invariance 
condition under rigid body rotations. The reason for this fact is the use of 
the wrong representation, eq. (2), of the general rotational movements. Any 
other relation which alters the nuclear co-ordinates alone without introdu- 
cing the corresponding changes in the nuclear internal parameters or vice 
versa, will not lead to an invariance relation. Evidently the translational 
invariance has a counterpart in the present theory. 


§ 4. Discussion 


The so-called rotational invariance condition (1) has been widely used in 
many applications of the general lattice dynamics, some of which are affected 
by the present note. In particular, one may mention the photoelastic) 
and the piezoelectric!®) behaviours of crystals. The calculations on diamond 
by Smith!) are not seriously affected because of the symmetry of the lattice 
as pointed out by Launay!!) also. 

Of special importance is the work of Herring!”). According to the current 
ideas, the method of homogeneous deformation can be applied only to finite 
lattices!8) and in such a case the expressions for the elastic constants differ 
from those obtained by the long wave methods. Using the translational and 
rotational “invariance” relations eqs. (1, 6), Herring showed that the two 
expressions can be reconciled if one considers the contributions from the surface 
terms. The present note invalidates the proof of the equivalence and brings 


4) O. Theimer, Proc. physic. Soc. 65, 38 (1952). M. Born and K. Huang, l.e., 
p. 374, The reason why Theimer’s expressions give zero value for the photoelastic con- 
stants for NaCl type crystals is that Theimer considers only the first order effects. 

5) M. Born and K. Huang (I. c.), p. 262. 

6) H.M. J. Smith, Philos. Trans. Roy. Soc. London 241A, 105 (1948). Mr. _ 
Rajagopal has kindly pointed to the author that Smith’s coupling coefficients Ca 
are in error by the omission of exp [i (9, + q, + 4.) 2/2]. This affects the eigenvectors 
and the Raman effect calculations. Z. Physik 158, 473 (1960). 

') C. Herring (unpublished, 1937). The author is indebted to Dr. Conyers Herring 
for giving him an outline of this work. k 

13) There appears to be some cloudy thinking about this matter also. The issue will 


be set at rest in Part IV, where a homogeneous deformation method will be developed 
which is convergent even for an infinite lattice. 
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into focus the difference between the two expressions. Such a difference cannot 
even be conceived of in the classical theory of elasticity which the current 
lattice dynamics alleges to establish on a sound atomistic basis. 

Returning to the wave propagation in the lattice, one notices that the 
square bracket symbols [x y, 8 6] (34, 35) take 36 independent values for a 
triclinic crystal, while the round brackets (« y, 6 6) (36, 37) and hence the 
curly brackets {a y, 6 6} (38) take 45 independent values. If each atom is 
a centre of inversion, the round bracket terms vanish and {a y, B 6} = {ß y, a 6} 
with only 36 independent values. (Strictly it is sufficient if a center of in- 
version at each atom together with some translations brings the crystal 
jnto self coincidence.) If one puts eq. (39) in the form 

then it is easily seen that this equation, corresponding to the classical Chri- 
stoffel equation involves 36 independent constants [{x y, 6 6} + {a 6, B y}}/2. 
From such facts, Laval!) and Viswanathan!) have concluded that a 
crystal must have 45 elastic constants and Krishnan et. al!®) have suggested 
numerous experiments for distinguishing between the exact equation (40) 
and the classical Christoffel equation. But in Part IV it will be shown that 
the matter is not so simple as was believed and that some remarks need be 
made on the works of Laval and Viswanathan. 

In Part IV of this work, a rigorous theory of elasticity will be developed 
which will give a rational explanation of the facts (eg., the difference between 
the constants in the homogeneous deformation and the wave propagation 
methods, the absence of additional symmetry in eq. (39) ete.) which plague 
the current theories of lattice dynamics. 


The author thanks Prof. R. 8. Krishnan for his encouragement and 
the €. S. I. R. for the award of a Senior Research Fellowship. 


2) R. S. Krishnan, V. Chandrasekharan and E. S. Rajagopal, Golden 
Jubilee Research Vol. Ind. Inst. Sci, p. 150 (1959). 


Bangalore-12 (India), Department of Physics, Indian Institute of Science. 
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The Existence of Interfacial Couples in Infinitesimal 
Elasticity 


By E.S. Rajagopal 


Summary 


The existence of interfacial couples in general crystalline media is estab. 
lished and an extension of the Cauchy stress principle is formulated in 
this paper. A strict method of treating the homogeneous deformation of an 
infinite lattice is also presented. The relevance of interfacial couples in the 
theory of elasticity is discussed together with a critical survey of the recent 
work in this field. 


§ 1. Introduction 


The theory of elasticity has as one of its fundamental corner stones, the 
Cauchy stress principle?) that in an elasic body the surface reactions acting 
across an elementary area can be represented by a force denoting the stress 
vector. The balance of the moments of the surface and body forces acting on 
an elementary volume leads then to the symmetry of the stress tensor and 
consequently the rotational components of the deformation of the volume 
elements do not play any part in the strain energy. Alternatively, one may 
start with the Boltzmann’s axiom?), that in an elastic body the surface 
reactions can be represented by a symmetric (stress) tensor. Though the 
implication of Boltzmann’s axiom is not evident, the axiom can be proved 
to be equivalent to Cauchy‘s principle and hence can form the starting point 
in the strict axiomatic development of the theory of elasticity. Indeed, it 
will be proved (see § 4) that even when body couples are present, a modified 
symmetric tensor can be found satisfying Boltzmann’s axiom, though some 
authors have fondly expressed their belief to be otherwise. 

The thermodynamic approach can be made by proving the existence of 
the strain energy function. But thereafter one has to make specific assump- 
tions about the form of the strain energy function®), which are equivalent to 
the Cauchy-Boltzmann axioms, though they do not have the same direct 
appeal. There are other reasons which force one to accept Cauchy’s principle 
as the most satisfactory basis for the study of elasticity in rational mechanics‘). 

In the present note, it will be proved that Cauchy’s postulate needs an 
extension by including interfacial stress couples. The consequences of this 


1) ©. Truesdell, J. Rational Mech. Math. 1, 157 (1952); A.E. Green and 
W. Zerna, Theoretical Elasticity, Oxford 1954, p. 63. 

2) G. Hammel, Handbuch der Physik Eds. Geiger and Scheel, V. p.9 (1927). 

», A.E. Green and W. Zerna, (l. ce.) p. 74, 155. 

4) C. Truesdell, (l. ec.) 1, 171 (1952); 2, 599 (1953). 
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jogical extension will be briefly discussed with a view to refute the various 
uncritical suggestions which have been recently made about the foundations 
ot fhe theory of elasticity. 

It is obvious that the justification or otherwise of the Cauchy-Boltz- 
mann postulates cannot be carried out within the framework of the continuum 
mechanics, but requires a detailed consideration of the atomic basis of matter. 
On the basis of the simplified model of central force interactions among the 
atoms constituting the crystall lattice, Born®) proved the Boltzmann’s 
axiom. Thereafter nobody seems to have paid any attention to this problem, 
in spite of the increased knowledge of the atomic interactions, especially the 
necessity of noncentral forces. 

Surprisingly, several current workers in lattice dynamics®) have expressed 
the view that the strain energy of a crystal cannot be defined in the general 
lattice theory. The naive arguments are as follows: Consider the homogeneous 
deformation of a crystal of unbounded dimensions, which one must consider 
if one wants to avoid the delicate conditions at the boundary atoms. The dis- 
placements increase with distance indefinitely as one moves from some origin 
and the energy change associated with this deformation is unbounded. Since 
this energy cannot be normalised to unit volume in the general lattice theory, 
it has been concluded that the homogeneous deformation method cannot be 
applied in the general theory. In addition, it is alleged that conditions must 
be imposed for the vanishing of the uniform ‘frozen’ stresses which can be 
present in the initial state. In Part II’), their contentions about the question 
of initial stresses were refuted and in the present note, their considerations on 
the applicability of the homogeneous deformations will be discussed. 

Strictly, the representation of the total potential energy of the infinite 
crystalline lattice in the form (II. 3) 


k 
ned as the derivative of ® which is unbounded. None of the workers, who have 
discussed infinities of deformation energy have paid attention to this. This 
doubtful procedure is completely avoided by introducing ®,, e ee as the 
force constants between the (m j) and (p k) atoms as has been done at the very 
outset of Part III. It is surprising that though the force constant aspect is 
known for a long time®), it has not been exploited at all. Amusingly, it is this 
force constant consideration which leads to the complete solution of the va- 


pm M. Born, Atomtheorie des festen Zustandes, Encyk. Math. Wiss. V-3, p. 546 
). 

*) K. Huang, Proc. Roy. Soc. London 208 A, 178 (1950); M. Born and K. Huang, 
Dynamical Theory of Crystal lattices, Oxford 1954, p. 237. 

?) The note “The Role of Initial Stresses in Lattice Dynamics” will be referred to 
as Part IT, the note “The Rotational Invariance condition in Lattice Dynamics” as Part IlI, 
and the present note as Part IV. For example, eq. (II. 3) refers to eq. 3 of Part II. The 
report “The validity of the recent theories of crystal elasticity (Proc. Ist Symp. Solid 
~~" u Bangalore. 1960) will be referred to as Part I and a forthcoming review 
as Part V. 

®) G.H. Begbie and M. Born, Proc. Roy. Soc. London 188 A, 179 (1947). 
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rious problems. For, one can apply a deformation to the crystal and calculate 
directly the stresses acting across an interface from the displacements and the 
force constants. Alternatively, since the energy associated with the relative 
displacements of two atoms is known, one can calculate the energy associated 
with the atoms of a single unit cell due to the displacements of all the atoms 
of the lattice. This energy change divided by the volume of the cell is directly 
the strain energy per unit volume. There seems to be an impression among 
the recent workers in lattice dynamics that all the quantities of interest must 
be derived from the (divergent) potential energy of the infinite lattice and it is 
very difficult to explain why the straightforward, almost trivial, procedures 
indicated above have not been used. After all in rational mechanics force 
isas fundamental a concept as energy. 


§ 2. Direct Evaluation of the Stress in an Atomistic Model 


It is well known that a homogeneous deformation of a lattice must be re- 
presented (with usual notation of Part IIT) as 


Born and others®)*) have used u,,, themselves to denote the strain 
components, but it will be necessary to use here the (infinitesimal) strain 
tensor Eup = (wx,6 + Ug,.)/2. The force on the (mj) atom due to the (pk) 
atom is given (IIT. 5) by 


+ a =m — p. (2) 


In the deformed state, the net force on (mj) is zero; so summing eq. (2) over 
all (p k) and using (III. 6) that 


= (3) 


5 Pap Pap = 0 (4) 
Remembering (IIT. 17) that = (u, ® =T, and that the 


one gets 


matrix J, E 5) is singular, it is obvious that all the 3 g quantities 0, (k) cannot 


be determined. But actually only the relative interlattice displacements are 
of physical interest and hence one can set 95 (k) = 0. Also, using (III. 29a), 
one has 


Now equations (4—5) give 


No confusion need arise between the dummy index u and the mass u,. 


E 
tion: 
the 
thes 
stud 
lat 
to 1 
are é 
oth 
sinc 
incl 
erys 
wit! 
forn 
| The 
bef 
incl 
| 
ot, 
The 
| car 

Sin 


E. 8. Rajagopal: The Existence of Interfacial Couples in Infinitesimal Elasticity 195 


This relation enables one to eliminate the 6 terms in the subsequent equa- 
tions. As is well known, physically the 6 term represents the readjustment of 
the partial lattices so as to lead to the equilibrium of the system. Though 
these interlattice displacements are known for a long time, they have not been 
studied in detail®) probably because they are macroscopically inobservable. 

The calculation of the stress in an atomic basis dates back to Cauchy”). 
later Poincare and others!) proved in a general way that the stress is equal 
to the suitable component of the force exerted on a suitably oriented unit 
area by one half of the infinite crystal on the infinite cylinder of crystal in the 
other half crystal with unit area as the base. The edge effects do not come in 
since the stress contribution will increase as L? while the fringing effect will 
increase L. The stress method has been used by Born?) and others. 

Consider the infinite crystalline plane with the normal ox, dividing the 
crystal into two infinite half crystals, say A and B. Take a primitive cell 
with the edges r,, t, in the x, plane and consider the semi-infinite parallelepiped 
formed by stacking similar cells in the r, direction into say the A semicrystal. 
The x component of the total force per unit area exerted on this parallelepiped 
by the semi crystal B gives t,,. Strictly this unit area refers to the state 
before applying the infinitesimal deformation and the forces refer to the 
incremental stress caused by applying the deformation. This point connects 
up the present discussion with Part II. The stress is thus 


a 


t A jk 


The summation over n the number of the cell in the parallalepiped can be 
carried out and yields the well known result 


1 
= A jk mo dy fo (5). 
Now £, (4) = &, (a) + &,(j) —&, (k), where &,(a)=a,2,. So replacing 
|ty A by v, one has 


Since f, =—f, (see discussion before III.9), one allows the summation 
over all a, and using in addition e q. (3) 


The force f, (Gi is caused by the homogeneous deformation and with eq. (2) 


Now (III. 34) gives 


culate 
1d. the 
lative 
ciated 
atoms 
rectly 
mong 
must 
d it is 
dures 
force 
be re- 
(1) 
strain 
strain 
| 
(2) 
| over 
(3) 
wf 
are 
29a), 
(5) 
(6) 
13* 


196 Annalen der Physik. 7. Folge. Band 6. 1960 


Moreover, 


which with eq. (6) and (III. 36) yields 
1 b b i hl a 
Consequently, using (III. 38) 
Stay = {La BO] + (ay; BO} una = BO} uno. (2) 


It is well to remember that {x y, ß 6} = {ß6, « y} and that there are 45 inde. 
pendent curly bracket constants in the general case. No interchange of « and 
y or of ß and ‘6 other than the above one is permissible. 

In terms of the infinitesimal strain tensor 


= BO} + {x y, B}] 0/2 = Cay, BS ep (8) 
<a y, Bd> = <a y, Öß) but <a y, B d> <y x, Od. (9) 


Consequently in the general case 


where 


tay < bys 
that is, Boltzmann’s axiom is not valid in the case of a triclinic crystal with 
general noncentral interatomic forces. 


§ 3. The strain Energy in lattice Dynamics 
Before discussing (see § 4) the implications of eq. (8—9), it is instructive 
to see how the strain energy density can be obtained in a straightforward way. 
Incidentally, this will complete the refutation of the statements prevalent in 
recent literature on lattice dynamics that the strain energy density cannot be 
formulated in the general theory. 


The force exerted by the displaced atom (p,) on (m,) is fa er z “a = 


p—m 
= Dis oe ak »( kj ) and if the atom (m,) undergoes a relative displace 


ment ua(7 5) the work done is shared between the two atoms. So the 


energy change associated with all the atoms in the mth cell is 

ly m—p p—m m—p 

This is independ ent of m i. e. , the energy change is the same for all cells, as 
it should be. The energy Are is 


eq. (1), one 
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Using eqs. (4) and 6), 
=> (i) & = 


agh day 


b b 
Dus (; Un,y Up, 6 
Now eqs. (III. 34). (III. 36) and (III. 38) give 
6U = 2,0 y, Una,y Up, 6- (10) 


This is entirely consistent with eq. (7) and this proves the contention that 

a rigorous and straightforward development for the energy density can be 
ven. 

F It is also to be noted that only the incremental energy density has been 
defined. Indeed, it is a well known result in the classical theory that one can 
pass from 6U = [ff tap dV to if one knows the 
state of zero €,g in the reference frame used. Indeed, it was discussed in Part IT 
that so long as one remains in infinitesimal deformations, the frozen stresses 
in the initial state do not at all matter and one will reckon all the quantities 
form this initial state. A similar situation exists in the new theory developed 
in the present papers (viz. ,the inclusion of surface couples), but since the method 
of homogeneous deformations cannot lead one very far in the elasticity 
theory (see next section) except to be of convenience in the lattice dynamic 
calculations, one must now investigate the general theory of elasticity. 


'$4. Some effects of interfacial couples in Infinitesimal Elasticity 


The axiomatic development of the theory of elasticity can be made on the 
basis of the following stress principle: The reactions across a surface 
elements dS in a strained elastic body consist of a force t dS (vector part of the 
stress) and a couple y ds (pseudovector stress). In addition there may be body 
forces T and body couples & (per unit mass) acting on the volume elements. 

Let t,„ denote the component of the force along the x, direction when the 
surface element has a normal n. Then, a classic line of argument, due to 
Cauchy, regarding the equilibrium of a tetrahedron with one face normal to 
n and the other three faces normal to the x, xg x, axes, gives the well known 
result 


= au) 
A similar consideration yields 

Yan= Yap (12) 
Here t,; and y,g denote respectively the tensors representing the vector and 


the pseudovector stresses. The balance of forces acting on a volume element 
(€ is the impressed external force) 
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gives, on using eq. (11) and the familiar Gauss’s theorem 


2 6x, 0 (T, ug F,) = 0. (13) 


The balance of the moments acting on the volume elements 
SS (Xa ton — Xp tan) ds + ff fo [xs (Ts — Fp) — xp (T.—F,)] dv 
+ SS Yynds + JJ feGapdv =0 
becomes, on using eqns. (11—13) and Gauss’s theorem 


The afy cyclic order is to be preserved in eq. (14) except in the dummy index. 

It is now convenient to point out some features of the results established 
by eqns. (8, 9). They state that when the reactions across an element of 
surface are calculated there exists a finite difference tz, — t,g which is due to 
the surface dependent terms. On the otherhand, eq. (1 4) considers the reac- 
tion acting over the bonding surface as well as the reactions over the enclosed 
volume and establishes a difference t;,—t, 3 contributed by surface dependent 
as well as volume dependent factors. In fact it is a very well known character 
of the calculations on the lattice models, that all interactions among the 
various atoms lead only to surface stresses while interactions with external 
sources alone lead to volume forces or moments. For example, taking the case 
of the gravitational attraction, the interaction between the various individual 
atoms lead only to surface forces while body forces are obtained from the 
interactions of the atoms with an external gravitational field (which in special 
cases, may be the smoothned potential of the aggregate of atoms). Hence 
the calculation of the difference t;, —t,g on the lattice model con- 
stitutes by itself the proof of theexistence of the surface couples. 
(See also Part V). 


But the fact that the interactions among the atoms lead only to a surface 
stress system has been overlooked by some recent workers who equate the 
difference t;,—txg (the lattice dynamic calculations which they have used 
to deduce this result are not free from defects, as will be shown in the sequel) 
to the volume couples and claim that the classical elasticity theory needs an 
extension. It will now be proved that the introduction of volume couples does 
not alter the structure of the classical theory and that the introduction of 
surface couples is an innovation. This is easily deduced from the consideration 
of strain energy. 

The work done by the volume and surface stresses when the displacement 
u is increased by 6u resides as dW the strain energy increment in the medium, 
W being reckoned per unit volume. 

Now 


[ff dwdv = ff du-t, dS + ff de-y, dS + [ff 
+ fff (I — F) dv. 
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Using eqns. (12, 14) and the commutativity of 6 and d/éx, ete. 
yudS + fff + Go de = Jf yon ds + [ff E deny G, dv 


ö 


dw, 

Ina similar way 
ou 
SS du-t, ds + fff du o dv = ae 
Consequently, 
du, dw, 

(=z) FEL (=) buy (tap (15) 

This can be put in a more convenient form by introducing the stress, strain 


components. 
Let tx be the symmetric part of t,; and r,; the antisymmetric part so 


that 


Where Tag = (tag+tpa)/2; Tap = (tag — tpa)/2. (16) 
The components of the strain tensor are, of course 
Ep = (Ux, 8 + Up, «)/2. (16a) 
With « By in cyclic order, 
Wy = Wap = (Up,a — Ua,p)/2. (16b) 
Then x 
~ lap Ux,p = ae ab (Tap + Tap) — Wag) 
ae &ap — (tag —tpa)/2 (17) 


=< Exp — (tap — tea). 


Ow 
ow = +E (18) 


It is well to remember that the expansion of the first term on the right side!) 

Equation (18) requires careful discussion. The first point to be noticed is 
that the volume couples do not enter into the strain energy. This provides a 
complete refutation of the recent suggestions that the classical theory 
neglects the volume couples and that their introduction leads to more than 21 
elastic constants for the triclinic crystal. This point has already been laid bare 


®) R.S. Krishnan, V. Chandrasekharan and E. 8S. Rajagopal, Proc. National 
Inst. Sci. India (under publication). 
") A. E. H. Love, Mathematical Theory of Elasticity, Dover, 1946. 
1) See, for example, I. S. Stakgold, Quart. Appl. Math. 8, 185 (1950). 
my Compare for example, H. Jeffreys, Cartesian Tensors Cambridge 1952 Chap. 
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by Tiffen and Stevenson’). In the absence of interfacial couples the 
elastic behariour is adequately described by t,, and e,g and in the rela. 
tionship between them there will be only 21 constants. Starting from eq. (18) 
one can obtain all the equations of the elasticity theory by the well know 
variational procedures and they will be identical with those of the usual 
classical theory except for slightly redefined shear stresses as proved by 
Tiffen and Stevenson. Secondly even in a static deformation departur 
from the classical theory can occur due to the presence of interfacial couples 
Since the rotation derivatives can be expressed in terms of strain derivatives, 
it is seen that homogeneous deformations require only 21 constants for their 
adequate description (as also seen from eq. 14), and only non uniform strains 
require more constants. This explains the notations of irrotational elastic 
constants introduced earlier by the present author’) in preference to the 
term static constants used by some other workers. The terminology irrota- 
tional strains for ¢,g and rotational strains for w.,s seems quite appro. 
priate. @,, is not necessarily a higher order infinitesimal as can be checked 
for a wave propagation. 


After the rigorous foundation of the elasticity theory by Cauchy, Green 
and others, the possibility of volume couples was considered by Maxwell, 
Duhem and others!®). A more recent discussion of the elastic dielectric is 
due to Toupin!®) who also shows the symmetry of the elastic constants in the 
presence of volume couples. The possibility of surface moments was first 
considered by E. Cosserat and F. Cosserat}®) who deduced the equivalent 
of eq. (14). But since rigorous atomic discussions and careful experimental 
evidences were not available, the classical theory came into universal acceptance. 


The modern phase of the foundation of elasticity is due to Born and his 
collaborators®)®) who first considered radial atomic forces and later general 
interatomic forces. While the radial force treatment is correct as far as it goes 
their general theory is plagued with inconsistencies as we have exposed in this 
series of papers. In the last decade. Laval!’), Viswanathan!8), Le Corre), 
Raman”), Wooster”) and others have suggested modifications of the 
classical theory, the first two making lattice dynamic studies also. 


13) R. Tiffen and A.C. Stevenson, Quart, J. Math. Appl. Mech. 9, 306 (1956). 
The present note corrects an error in this paper. Apart from some notational slips in 
their eq. (5), Tiffen and Stevenson assume in their eq. (12) y„a = — Ya. (1 —2 4,4) 
for which there is no basis. In a private letter, Dr. Tiffen has agreed to this point. 

14) R.S. Krishnan, V. Chandrasekharan and E. S. Rajagopal, Golden Jubilee 
Res. Vol. Ind. Inst. Sci., p. 150 (1959); Proc. Nat. Inst. Sci. India (in press). These 
papers are formally based on the Laval-Le Corre — Viswanathan-Raman theory 
and even though that basis is proved to be wrong in the present note; it will be seen that 
the whole set of calculations can be reformulated on the basis of the present note. 

15) The references may be found in C. Truesdell, J. Rational Mech. Math. 1, 158 
(1952). 

16) R. A. Toupin, J. Rational Mech. Math. 5, 909 (1956). 

17) J. Laval, L’etat Solide, Solvay Congress, Brussels 1951, p. 273. 

18) K.S, Viswanathan, Proc. Indian Acad. Sci. 89 A, 196 (1954); 41 A, 98 (1955). 

1%) Y. Le Corre, Elasticite et piezoelectricite cristallines, Thesis, Paris 1953, J. 
Physique Radium 19, 541 (1958). 

2) C.V. Raman and K.S. Viswanathan, Proc. Indian Acad. Sci. 42 A, 1, 5l 
(1955). 

21) W. A. Wooster and N. Joel, Nature 180, 430 (1957); 182, 1078 (1958). 
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As far as the lattice dynamic studies are concerned, both Laval and Vis- 
wanathan have not met the question of initial stresses (Part II). In fact, 
Viswanathan assumes a finite crystal, when he must have considered 
Herring’s work (Part III § 4). Secondly while Laval in 1957, has recognised 
the failure of the so called rotational invariance condition (Part III) he has not 
given any proof for this. Viswanathan has completely overlooked this failure. 
He applies this condition at one place, but fails to notice that the condition 
will leave only 36 independent constants and not 45 as correctly pointed out 
by Laval. Both Laval and Viswanathan have not noticed the fallacy 
in the description of lattice rotations as considered in part III § 1. 

As far as the continuum theories are concerned all these authors attribute 
the departure from the classical theory to the volume couples. This is cer- 
tainly completely wrong, as conclusively established by Tiffen and Ste- 
venson and in the present note. Their theories lead to a strain energy in- 
volving the components of rotation w,, (as in eq. 17) so that a rigid body ro- 
tation will alter the energy content, a situation which has little semblence to 
reality. The correct procedure is that only differential rotations, just as diffe- 
rential displacements, are involved in the strain energy. 

It may be added that Le Corre*) has once mentioned surface couples, 
but he has not developed the matter correctly. McClintock, and his co- 
workers**) have considered the effect of surface couples (eq. 14), without 
realizing the possibility that such couples may be present in perfect lattices 
also. 

Reissner) had discussed a case where the assymmetry of the stress 
tensor contradicts the postulate of the continuity of the strain. The present 
author has not been able to prove (for an assymmetric stress tensor) the 
existence of the Airy’s function as assumed by Reissner. 

Huntington*) has recently reviewed the theory of crystal elasticity 
without noticing the points raised in the present papers. In particular, the 
anomaly noted by him on p. 281 is due to a notational slipin Kellermann’s 
paper as pointed out by Blackman”), 
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Uber die Polarisierbarkeit und Suszeptibilitat 


freier Atome und lonen in der Theorie 
von Fermi und Amaldi 


Von T. Tietz 


Inhaltsiibersicht 


In dieser Arbeit berechnen wir die Polarisierbarkeiten wie auch die Sus. 
zeptibilitäten von freien Atomen und positiven Ionen. Mit Hilfe eines glück- 
lichen Ansatzes für die Elektronendichte des Atoms oder Ions ist es in dieser 
Arbeit dem Verfasser gelungen. gute Ergebnisse für die Polarisierbarkeiten 
wie auch für die Suszeptibilitäten von freien Atomen und Ionen im Vergleich 
mit der Erfahrung zu bekommen. In Tab. 2 und 3 werden unsere Ergebnisse 
mit anderen theoretischen, wie auch experimentellen Werten verglichen. 
Zum Vorteil unserer Methode erhält man für das atomare Mittel von r?, eine 
einfache geschlossene Formel. 


Die exakte Berechnung der Polarisierbarkeit von Atomen und Ionen mit 
mehr als einem Elektron ist zur Zeit wegen mathematischer Schwierigkeiten 
unmöglich. Bei der quantenmechanischen Berechnung der Polarisierbarkeit 
muß man die Kenntnis der Eigenfunktionen und der Eigenwerte aller Elektro- 
nenzustände voraussetzen, deren exakte Berechnung zur Zeit noch aussteht. 
In einer kürzlich erschienenen Arbeit hat Gombäs!) die Polarisierbarkeit 
freier Atome in einem erweiterten statistischen Modell mit Hilfe einer Störungs- 
rechnung berechnet. Mit Hilfe dieser Störungsrechnung gelangt Gombäs 
für die Polarisierbarkeit der freien Atome zu einer Formel, die keinerlei 
empirische oder halbempirische Parameter enthält. Die aus dieser Formel 
berechneten Polarisierbarkeiten der freien Edelgasatome Kr und Xe stimmen 
mit den empirischen Werten befriedigend überein, für die leichteren freien 
Edelgasatome Ne und Ar erhält man jedoch viel zu große Werte. In dieser 
Arbeit leiten wir einen geschlossenen Näherungsausdruck für die Polarisier- 
barkeit von freien neutralen Atomen und von Ionen mit Hilfe der Kirkwood- 
schen?) Formel ab. Wie bekannt gilt nach Kirkwood für die Polarisierbarkeit 
x von freien neutralen Atomen und Ionen 


4 — 
(r2)2, (1) 
wo N die Zahl der Elektronen des Atoms oder Ions und a, den kleinsten 
Borschen Wasserstoffradius und r? das Mittel über r? bezeichnen. Für r 


1) P. Gombas, Acta Physica Hungarica 9, 101 (1959). 
2) J.G. Kirkwood, Physik. Z. 33, 57 (1932). 
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gilt in der statistischen Theorie von Thomas und Fermi?) 
oo 
r= [0 (r)dv = 4a r8Q(r) dr. (2) 


Hierin bedeutet o(r) die Elektronendichte. Es ist zweckmäßig, solche Modelle 
für o(r) zugrunde zu legen, bei denen die Austauschkorrektionen bzw. Korre- 
lationskorrektionen für die äußeren Gebiete des Atoms durch die Fermi- 
Amaldische Korrektion ORDER werden. Dann ist: 

N y" (x) 


o(r )= (3) 
Z ist die Ordnungszahl des Atoms oder Ions und es ist x =, , wo 
Die Th Ba ch Funktion y(x) genügt folgenden Randbedingungen: 
Z—N+1 
y(0) =1, y(%) =0 und %:% (%) = (5) 


Diese Randbedingungen lassen sich für die Approximation des Verfassers*) 
leicht befriedigen®). Für y(x) haben H in diesem Falle folgende Näherung 


1 (; +a 
(6) 
wo a eine Konstante bezeichnet. Die Br Y(%) = 0 ist exakt 
erfüllt. Bei 2 = 0 muß noch die Randbedingung y(0) = 1 erfüllt sein. Unsere 
Approximation y(x) gibt dagegen 


1 
= (7) 


da aber ax, erheblich größer als 1 wird, unterscheidet sich die rechte Seite 
nur ganz belanglos von 1. Als dritte Randbedingung erhalten wir 
bax —N+1 
Diese Beziehung erlaubt uns in einfacher Weise die Berechnung von x, 
und damit wegen 2 —=r/u auch die der Radien ry. Für kugelsymmetrische 


Elektronenverteilung erhalten wir für r? mit Hilfe von (2) und (3) durch 
partielle Integration folgenden Ausruck: 


% 
| ry" (x) da = A [(23 y'(a)|2— (3.2 +6 f xy(x)dx] (9) 
0 


mit 


(10) 


Hieraus ergibt sich durch Einsetzen unserer Approximation bei Beachtung 
der entsprechenden 


Z—N+1)2 


6 3 (11) 


3) P.Gombäs, Die statistische Theorie des Atoms und ihre Anwendungen. Wien, 
Springer- Verlag (1949). 

‘) T. Tietz, Nuovo Cimento 1, 968 (1955). 

Tietz, J. chem. Physics 26, 430 (1957). 
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| ‘2%@8 Die Polarisierbarkeit folgt dann aus der der Kirkwoodschen Formel (1). 
“af „ist aus Gl. (8) zu bestimmen. Für die bisher noch freie Konstante a werden 
folgende vom Verfasser bereits in®) gründlich diskutierte Zahlenwerte ange- 
Sma | nommen: a =($) = 0,536257 oder. 1/a = 1,7566. Tab. 1 gibt einen Ver- 
gleich unserer x,-Werte für verschiedene Ordnungszahlen Z mit den numeri- 
schen Werten von Fermi’) und Kobayash 
— Ersichtlich ist die Übereinstimmung ziemlich gut. Tab. 2 enthält die 
| 22% § Polarisierbarkeit « von freien Edelgasatomen und Ionen pro Grammatom in 
Einheiten 10-4 cm? und gibt einen Vergleich unserer «-Werte mit theoreti- 
schen und experimentellen Ergebnissen von anderer Seite. 


gue Ersichtlich sind unsere Ergebnisse nicht schlechter als die von anderen 
#5 mit verschiedenen statistischen Modellen numerisch berechneten a-Werte. 
man Für die diamagnetische Suszeptibilität von Atomen und Ionen pro Gramm- 
Naga fp atom gibt bekanntlich folgender Ausdruck 

L e2 — 
258 1 (12) 


wo L die Loschmidtsche Zahl, c die Lichtgeschwindigkeit, m und e Masse 


~3%% § und Ladung des Elektrons bezeichnen. r? ist auch in diesem Falle durch (2) 
gegeben. Die berechneten Suszeptibilitäten von freien Edelgasatomen und 
freien Ionen sind in Tab. 3 angegeben. Hier sind auch unsere y-Werte mit den 
auf andere Weise numerisch berechneten verglichen. Zum Vergleich haben 
wir einige von Angus°) mit Hilfe halbempirischer Eigenfunktionen berechnete 
ot Werte angegeben, an denen Eucken noch kleinere Korrektionen angebracht 
hat. 
Auch in diesem Falle stimmen unsere x-Werte für die freien Atome und 

=a! Ionen gut mit dem empirischen Werten überein. 
Hits Aus dieser befriedigenden Übereinstimmung unserer x- und y-Werte mit 

«+. 4 der Erfahrung kann man den Schluß ziehen, daß unsere Näherungsformel fiir 
|S! die Elektronendichte des Atoms oder Ions zufolge eines glücklichen Ansatzes 
für o(r) eine gute Annäherung an die wahre Elektronendichte ist. Die theore- 
tische Begründung der korrigierten Elektronendichte im Thomas-Fermi- 
schen Modell ist ein schwieriges Problem. Einen Überblick über die Korrek- 
tionen und Erweiterungen des statistischen Modells findet man in dem Buche 
von Gombäs. 


Der Verfasser dieser Arbeit ist Herrn Professor D. P. Gombäs für einen 
Besuch in seinem Institut sehr zum Dank verpflichtet. 


5) Vergleiche hierzu (4) 
?) Vergleiche hierzu (3) 
8) S.Kobayashi, J. Phys. Japan 14, 1039 (1959). 
*) Vergleiche hierzu (3). 


Lödz (Polen), Institut für theoretische Physik der Universität. 


Bei der Redaktion eingegangen am 7. Dezember 1959. 


Nach Angus halbempirisch (Korrigiert) . 


Experimentelle Werte. . . . 
Unsere Werte (1/a = 1,7566). 


Theoretische Untersuchungen 
über das Verhalten geladener Teilchen 
in Sattelpunkten elektrischer Wechselfelder 


Von Adolf Müller 
Mit 17 Abbildungen 


Inhaltsübersicht 


In dieser Arbeit wird die Bewegung eines geladenen Teilchens in elektri- 
schen Gleich- und Wechselfeldern unter Berücksichtigung der Reibung und 
des Erdfeldes untersucht; insbesondere interessieren periodische Lösungen 
der (nicht-linearen) Bewegungsgleichung auf den Symmetrieachsen von eben- 
und axialsymmetrischen Potentialfeldern in Sattelpunktnähe. Es werden 
unter anderen die Bedingungen angegeben, unter denen sich das geladene 
Teilchen im Sattelpunkt eines Wechselfeldes ,,frei aufhängen‘ läßt. Auch das 
Verhalten von zwei gleichen Teilchen in der Symmetrieebene axialsymme- 
trischer Potentiale wird kurz untersucht. Zur Bestimmung von Ladung und 
Masse des Teilchens ergeben sich neue Meßmethoden. 


I. Einleitung 


Nach Experimenten von H. Straubel!)?2) können elektrisch geladene 
Teilchen mit Massen zwischen 10-8 und 10-5 gr (Durchmesser etwa 3 - 10° bis 
10? cm) und Ladungen um 10-12 As in einem Drahtring oder einem System 
von zwei parallelen Drähten frei „aufgehängt‘‘ werden, wenn man an die 
Drähte eine niederfrequente Wechselspannung (bis etwa 20 kV) anlegt. Man 
stellt verschiedene Typen von periodischen Bewegungen mit der ganzen und 
auch der halben Frequenz des Wechselfeldes fest, welche von der Zahl der 
Teilchen, der gewählten Anordnung und den Werten von Ladung, Masse, 
Reibung, Frequenz und Amplitude abhängen. Die Teilchen schwingen zum 
größten Teil auf Symmetrieachsen oder in Symmetrieebenen der Anordnungen 
in der Nähe eines Sattelpunktes. Das Erdfeld bewirkt ein ,,Durchhangen“ der 
Teilchen, bei mehreren Teilchen kommt die gegenseitige Abstoßung hinzu. 
Als außerordentlich günstig für Messungen empfiehlt sich die Anordnung eines 
Kreislochblendenkondensators. Es soll im folgenden versucht werden, einige 
dieser Phänomene mathematisch zu erfassen. Die Bewegungsgleichungen dieser 
Teilchen sind im allgemeinen nichtlineare Differentialgleichungen mit periodi- 
schen Koeffizienten. Ihre Behandlung führt in den relativ jungen Zweig der 


1) H. Straubel, Naturwissenschaften 1955, 42. 
2) H. Straubel, Z. Elektrochemie 60, H. 9/10 (1956). 
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theoretischen Mechanik, die nichtlineare Mechanik*%-14). Leider existiert 
gegenwärtig noch keine Theorie zur Lösung der uns hier interessierenden Dif- 
ferentialgleichungen. Unser Vorgehen wird daher so sein, daß wir den von 
H. Straubel durchgeführten Experimenten Lösungsansätze entnehmen und 
dann die Parameter bestimmen, wobei wir uns auch auf Lösungstypen der 
(linearen) Mathieuschen Differentialgleichung stützen können. 


II. Die Bewegungsgleichung eines geladenen Teilchens 


Es wird die Gültigkeit des linearen Reibungsgesetzes —x - it angenommen 
(x > 0: Reibungskonstante, r: Ortsrektor, +: Geschwindigkeit des Teilchens). 
Das elektrische Feld € (r, ¢) wird aus einem vorgegebenen Potential der Art 
g (t, t) = p= (t) + 9 (t) - cos w t berechnet (quasistationär, w: Kreisfrequenz, 
t: Zeit). Ist g die (gerichtete) Erdbeschleunigung und €, ein zusätzliches homo- 
genes Gleichfeld, so lautet also die Bewegungsgleichung eines geladenen Teil- 
chens 


mt + xt —e (r,t) + G&}—mg =0 (1) 


(m: Masse, e: Ladung des Teilchers). Mit der Transformation w -¢ = 2 t und 
den Abkiirzungen 


ergibt sich aus (1) 
+ 2yr’ + 20 {grad (t) + cos2r- grad (t)} 2) 


ein System von 3 im allgemeinen nichtlinearen, gewöhnlichen und inhomogenen 
Differentialgleichungen 2. Ordnung mit periodischen Koeffizienten. 


II. Lösungen bei linearer Annäherung des Feldes in Sattelpunktnähe 


Wir beschränken unsere Betrachtungen auf Schwingungen in Symmetrie- 
Achsen oder -Ebenen. Wird mit q die entsprechende transformierte Koordinate 


y a: u .Stoker, Non-linear vibrations in mechanical and electrical systems, New 
or 

*) N. Minorsky, Introduction to Non-linear Mechanics, Ann Arbor (USA) 1947. 
ie’ N. Krylov u. N. Bogoljubov, Introduction to Non-linear Mechanics, Princeton 


int S. Lefschetz, Contributions to the Theory of Non-linear Oscillations, Princeton 


”) N. W. McLachlan, Ordinary non-linear differential equations in engineering and 
physical sciences, Oxford 1950. 

8) A.A. Andronow u. C. E. Chaikin, Theory of Oscillations, Princeton 1949. 

®) Proceedings of the Symposium on Non-linear Circuit Analysis, Polyt. Inst. of 
Brooklyn, New York 1953 (April). 

”) I.G. Malkin, Metoda Liapunowa i Poincare w Teorii Nelinejnych Kolebanij, 
Gostechisdat 1949. 

1) K.F. Teodortschik, Awtokolebatelnye Sistemy, Gostechisdat 1952. 

”) K. Klotter, Neuere Methoden und Ergebnisse auf dem Gebiet der nichtlinearen 
Schwingungen, VDI-Berichte Bd. 4 (1955). 

8) G.J. Cohn u. B. Saltzberg, J. Appl. Phys. 24, 2 (1953); 25, 2 (1954). 

4) L.W. Kantorowitsch u. W.J.Krylov, Näherungsmethoden der höheren 
Analysis, Berlin 1956. 
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bezeichnet, ist die Bewegungsgleichung des Teilchens in diesem Falle vom Typ 
q(t)” + 2y q(t)’ + 28 - cos 2t)-F (q)+k=0, (3) 
F (q) charakterisiert den Feldverlauf, « ist der angelegten Gleichspannung 


U_ und ß der Wechselspannung U. proportional. In Sattelpunktnähe ist 
bei sehr kleinen Auslenkungen F (q) = g, so daß mit 


q’ +2yq' + (a—2B-cos2t)g+k=0 (4) 
eine inhomogene Mathieusche Differentialgleichung mit Reibungsglied 
vorliegt. Mit y=k=0 wurde (4) bereits weitgehend untersucht) 
deshalb führen wir nur die für weitere Betrachtungen wichtigen Ergebnisse an. 
Die beiden linear unabhängigen Lösungen sind fiir y = k = 0 


+ co 
41.2 =. Con (x, B) eitanz»)r, (5) 


falls der ,,charakteristische Exponent‘ »—=A-—iu reell) keine ganze 
Zahl ist. Ist u = 0 (d. h. y reell), sind beide Lösungen beschränkt, man spricht 
von „stabilen“ Bereichen in der «, ß-Ebene (siehe Abb. 1, schraffiert). In den 

és nichtschraffierten Be- 


reichen ist v komplex, beide 
Lösungen also in —oo<r 
WUT <oo unbeschrankt (u + 0, 
| mG YY 4 » A ganze Zahl). Diese Be- 
a” 2 DE reiche werden durch Grenz- 
| Ga kurven &.n (ß) bzw. Xn (8) 
getrennt, auf denen =4 
N I eine ganze Zahl ist; q, und 

0 s N. . 
v. EN 9 sind dort linear abhän- 
Gey | gig; die linear unabhangige 
s / Lösung von g,, ist nicht 
beschränkt (~ r). Auf &,2n 
0 und &s2en+2 gibt es Lésun- 
Poel gen der Periode x, auf 
Abb.1. Stabilitätskarte der Mathieuschen Differen- %s2n+ı Und Oegn+1 der 
tialgleichung Periode 2 r(n = 0,1,2...). 


Falls y > 0, geht (4) mit 
k=0 durch die Transformation q = e-’*-@ in eine Mathieusche 
Differentialgleichung mit 4 = « — y?, 8 =A über. Für g ergeben sich damit 
„Dämpfungsgrenzkurven“, gekennzeichnet durch || = y = u (x — )]. 
Auf ihnen gibt es eine stabile periodische Lösung, die andere klingt ab. Ist 
|| >y, sind die Lösungen unbeschränkt (instabile Bereiche); ist dagegen 
|#| <y, klingen beide Lösungen ab (q — 0, „dämpfungsstabilisierte‘‘ Bereiche). 
Die Abb. 2 zeigt die Grenzkurven (mit «7 bezeichnet) fürn =1beiy = y 0,004; 


0,2; y0,2 und y0,4 (gestrichelte Kurven hinreichend genaue Werte aus”), 
15) J. Meixner u. F.W.Schäfke, Mathieusche Funktionen, Berlin-Göttigen- 


Heidelberg 1954. 
1) N.W.McLachlan, Theory and Application of Mathieu Functions, Oxford 


47. 
17), M.J.O.Strutt, Lamésche, Mathieusche und verwandte Funktionen, Berlin 
18) G. Kotowski, Z. angew. Math. Mech. 1948, 23. 
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ausgezogene Kurven nach der 1. Näherung f* — (af — 1)? = 4 y?). Der Fall 
k + 0 gibt Anlaß zu einer partikulären Lösung g, der Periode z (siehe auch }§)), 
die im allgemeinen beschränkt ist; 


auf den Grenzkurven, wo Lösungen 23 

der homogenen Differentialgleichung £A 

mit der Periode x vorkommen, tritt ag r 

jedoch ,,Resonanz“ auf, dort geht q, 20 

beiy = 0 mit t? und bei y>0 mit SG KL / 

1 gegen Unendlich. 7 
Zusammenfassend kann man sagen: 15 zZ Z / 

für |äl<y klingen alle Lösungen [ 

von (4) ab und es bleibt nur die / [ /' 

partikuläre Lösung g, der Periode x; 10” l l 

fir |ä| > y wächst immer eine Lösung ,, u ae | a 

unbegrenzt; für || =y ist eine Lö- ! 

sung periodisch, die andere klingt ab; \ \ 

auf a3, zeigt q, Resonanzverhalten N\ \ 

(n =0,1,2...). Im Falle einer eben- N \ 

symmetrischen Anordnung (x, y-Ebene, SN q » \ 

Feld soll nicht von z abhängen, z.B. ? AW \ 

Paralleldraht-System) sind die Be- N \\ \ 

wegungsgleichungen eines Teilchens N NN 


im Sattelpunkt in linearer Näherung -05 ’ 

(z soll nicht weiter interessieren) 05 
"“19y a2’ —(x—2B-cos2t)x 

0, (6) Abb. 2. Grenzkurven a” für = 1 und 

” , y = 0,004; 0,2; /0,4 in 1. Nähe- 

y y +(a—2B-cos2t)y (7) rung (gestrichelt) 


Liegt mit « auch — x in stabilen Bereichen, d. h. ist |u (a — y?, B)| << y und 
zugleich |u (— a — y®, B)|<y, dann klingen die Lösungen der homogenen 
parametergekoppelten Gleichungen ab, es bleiben nur die partikulären Lö- 
sungen der Periode x. Fallsk, = k, = 0, sind dann x = 0 und y = 0 asympto- 
tisch stabile Lösungen!®). Dies bedeutet physikalisch, daß man das geladene 
Teilchen im Sattelpunkt frei „aufhängen“ kann (siehe auch!)2)20). Nur mit 
Gleichspannung (8 = 0) ist ein „Aufhängen‘ nicht möglich, da bei ß = 0 
mit « nie zugleich —« in einem stabilen Bereich liegt (Abb. 1), nur mit Wechsel- 
spannung ist dies jedoch für bestimmte -Bereiche möglich, z. B. 0 <|ß|<ß% 
(Bn: Schnittpunkt der Dämpfungskurve a} mit der ß-Achse). Für 6 erhält 


man in 1. Näherung By = V 1+ 4y®, in 2. Näherung 


(Bx)? = (99 + 12 y%) — VE (99 + 12 — (1 + 474) (81 + 369%) 
0,8245 + 472 (8) 


19) Es wird eine Lösung g, dann als stabil bezeichnet, wenn eine „Störung“ öq fiir alle 
T> 1» beschränkt bleibt; geht außerdem dg — 0 für t > oo, ist gg asymptotisch stabil. 
2) R. F. Wuerker u.a., J. Appl. Phys. 80, 3 (1959). 
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(bis py ee 2 ist der Fehler kleiner als 1%, bei 8} ~ 1 nur 0,03%, BY = 0,90805). 
Wird bei a = 0 |8| > ß%, so werden die Schwingungen in x- und y-Richtung 
gleichzeitig instabil. 

Im Falle einer axial-symmetrischen Anordnung (r,y, 2: Zylinderkoordi- 
naten) sind die Bewegungsgleichungen in linearer Näherung im Sattelpunkt 
(mit y’ = 0 und nur k, + 0) 


2’ + 2y2’'+ («—2B-cos2t)z+ k, =0, (9) 
r’ + (10) 


Hier müssen mit x und f auch — 3 « und — } in stabilen Bereichen liegen, 
damit das Teilchen „aufgehängt“ werden kann. Ein Blick auf Abb. 1 zeigt, 
daß man das Hauptaugenmerk auf den 1. stabilen Bereich zu legen hat, 
d.h., es muß wieder 0 < < sein bei « = 0 und y <1, damit das 
Teilchen ,,schweben“ kann. Es ist zu bemerken, daß die Lösung z= 0 (bei 
k, = 0) eher als r=( instabil wird, wenn man die Wechselspannung (ß) 
erhöht, da sich die Parameter um den Faktor — # unterscheiden (r = 0 ist 
also noch für 0 < |ß| < 2% stabil). Dieser Sachverhalt wurde experimentell 
bestätigt. 

Sind die Auslenkungen aber nicht mehr sehr klein, müssen die nicht- 
linearen Glieder der Bewegungsgleichungen berücksichtigt werden, was im 
folgenden geschehen soll. 


IV. Periodische Lösungen auf den Symmetrieachsen der Anordnungen in 
nichtlinearer Annäherung des Feldes 

Das Experiment zeigt, daß auf den Symmetrieachsen der Anordnungen 
unter bestimmten Bedingungen periodische Schwingungen eines Teilchens 
mit größerer Amplitude möglich sind, obwohl nach der linearen Näherung die 
Lösungen dann unbeschränkt wachsen sollten. Dies ist eine Folge der Nicht- 
linearität des Feldes, wodurch die Parameter « und ß effektiv verändert 
werden. Wie man nun leicht ableiten kann, läßt sich der Feldverlauf auf den 
Achsen in der Nähe des Sattelpunktes q = 0 (|g|< 1) durch 


F_(q-) oder Fy (1 + (11) 
darstellen (q.. entsprechend transformierte Koordinate). F_ entspricht einer 
„weichen“, F, einer ‚harten‘ Rückstellkraft. Somit ist bei nichtlinearer An- 
näherung des Feldes die Differentialgleichung 

q’ + 2yq' + (x— + k=0 (12) 
zu untersuchen, da (12) mit  — —q? und k? + —k? auch q, beschreibt. 

Wir gehen nun so vor, daB wir die Schwingung mit bekannter Frequenz 
ansetzen und die Amplituden in Abhängigkeit von den Parametern in Näherung 
ermitteln. Ohne Berücksichtigung der Reibung ist dies eine gewisse Willkür, 
da im Experiment die erforderlichen Anfangsbedingungen wohl kaum realisiert 
werden können. Bei Reibung jedoch stellt sich die periodische Lösung 
(Periode x oder 2) ein. Die äußerst umständlichen Rechnungen werden hier 
ganz weggelassen und nur einige Ergebnisse angeführt. Es wird angesetzt 


© Ay + A, cost + B,-sint + 4A,:cos2r-+ B,-sin2r, (13) 
q- © a + cost + b,-sint + a,- cos 27+ by-sin 2r. (14) 
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1. Fall: Lösungen der Periode 2 in 1. Näherung bei y = k = 0 (Ay = 
Man erhält q_ B,-sint nahe «,, mita+ß 21; 


q- A,- cost nahe a,, mit a—ß 21; (15) 


6,-sint nahe a, mit «+P <1; 
g+ a,-cost nahe «,, mit <1; 

hierbei ist 


a? = — A?; b? = — B?. Die Abb. 3 zeigt B, über der angedeuteten «a, P- 


1 
Ebene. Man erhält die cos- aus der sin-Schwingung, wenn man ß durch —ß 
ersetzt. Die Stabilität der Lösungen wurde untersucht, das Ergebnis zeigt 


Abb. 3. Amplitude B, über «, $-Ebene 


Abb. 4. Ganz allgemein kann man sagen, daß die Schwingungen stabil sind, 
wenn &, f in ,,instabilen‘‘ Bereichen der linearen Differentialgleichung (Abb. 1) 
liegen; anderenfalls sind sie instabil (was experimentell bestätigt wurde). 
2. Fall: Lösungen der Periode 2x in 1. Näherung bei k = 0, y > 0. Mit 
A, = A, = B, = Oergibt sich für A,/B, eine Bestimmungsgleichung 4. Grades. 
Speziell für x = 0 erhält man (8 > 0; y!< 1) die Näherung 
q ~C,-sin(t +4) mit t 
1 x £%X i+ B’ 
B(1 + 


Die Abb. 5 zeigt C, für y? = 0 und y? = 0,025. Für q_ existiert diese Lö- 
sung nur für £2 — («x — 1)? > 4° (im ,,instabilen“ Bereich der Abb. 1) und 
ist dort asymptotisch stabil. Für q, müssen «, ß im „stabilen“ Bereich liegen, 
jedoch ist diese Schwingung instabil und nur g, = 0 dort asymptotisch stabil. 
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Solange also (bei 1 =0) 0<|B|<P% gilt, ist q_=0 die sich für 
Tt + co einstellende asymptotisch stabile Lösung. Diese wird für |ß| > 
instabil und q_ beginnt zu wachsen, bis der nichtlineare Faktor (1 — q?) den 


ans! 
70,025 
i 
Abb. 4. Bereiche der möglichen 0 
stabilen und instabilen Schwin- 05 10 15 B 20 
gungen in der «, ß-Ebene (1. Nahe- Abb.5. Amplitude C, mit y? = 0 und 0,025 in Ab- 
rung) hängigkeit von ß 


Wert von ß effektiv so verkleinert hat, daß wieder eine periodische Lösung 

möglich ist, die wegen der Reibung für große t auch angenommen wird. Die 

Abb. 5 zeigt, daß die Amplitude C, bei fy sprunghaft einsetzt (x = 0), es gilt 

aC, 

pay B 

Hiermit ist eine sehr empfindliche Meßmethode gegeben (siehe Abschnitt VIIIa). 

Zur Art der Lösung q, ist zu sagen, daß (z. B. beix = 0)g; = 0 für0 <|B|<pr 

asymptotisch stabil ist. Im Intervall 0 < |p| < By ist auch noch eine Qa 

periodische Lösung möglich, sie ist aber instabil, wie im Experiment bestätigt 
wurde. 

3. Fall: Es wurden ebenfalls die z-periodischen Lösungen bei y = k = 0 
in 2. Näherung untersucht in der Nähe der Grenzkurve &,, ~ — }/?. Es 
wird jedoch nicht näher darauf eingegangen, da sie in diesem Zusammenhang 
geringeres physikalisches Interesse haben. 

4. Fall: Die ,,partikulare“ Lösung q, unterscheidet sich nur geringfügig 
(bei gf <1) von der partikulären Lösung der linearen Differentialgleichung 
(z. B. fiir « = 0,0 < |p| < py: 


I cos (2 T +n}: 


=oo. (17) 


und zwar ist gp, etwas kleiner und g,_ etwas größer, z. B. (x = 0) 
= A,, + Ay, 008 27 + B,-sin2t; —y-A,,; 


Man beachte a +k? >0 und # + —k2 < 0! Beide Lösungen g,, und 
q,. sind asymptotisch stabil, sie wurden im Experiment beobachtet, die Art 
der Abweichung durch Messungen bestätigt. Das Erdfeld bewirkt also das 
„Durchhängen‘“ der Teilchen. 
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Dies bewirkt, daß die Schwingung der Periode 27 nicht genau auf der 
Dämpfungsgrenzkurve einsetzt, da das nichtlineare Glied die Parameter 
effektiv ändert. Für |k?| < 1 (für Erdfeld sehr gut erfüllt) wird (beix = y = 0) 

= 


wobei also jetzt die Schwingung bei || = - einsetzt. 


V. Spezielle symmetrische Potentialfelder | 


Den bisherigen Untersuchungen lagen ganz allgemein symmetrische 
Potentiale zugrunde, jedoch beschränkten wir uns bei der Diskussion der Be- 
wegungsgleichungen auf die unmittelbare Umgebung eines Sattelpunktes. 
Wir wollen nun größere Amplituden betrachten und auf spezielle Potential- 
felder eingehen, die im Experiment Verwendung finden (Paralleldrahtsystem, 
Kreisring, Schlitzblenden-Kondensator und Kreislochblenden-Kondensator). 
Es sollen kurz die Koordinatentransformationen (auf g), die Parameter «, 


und die den Feldverlauf auf den Achsen beschreibenden Funktionen F (q) : 


angegeben werden. 
a) Paralleldrahtsystem (siehe Abb. 6) 
Die Drähte sollen das Potential U_ + U coswt haben und sich in 


einem großen, geerdeten Hohlzylinder vom Durchmesser 2 D befinden. Dann 
ist 


88eUl, 4SeU,, 
Der Feldverlauf auf den Achsen ist F = ET. fiir g =y/l und F= I für 


q=2/l. 


Abb. 6. Paralleldrahtsystem Abb. 7. Kreisring 


b) Kreisring (Abb. 7) 
Beim Kreisring ist 


auf der z-Achse gilt für Werte g = \> . = < 0,5 ebenfalls mit ausreichender 
Genauigkeit F (q) ~q/(1 (Fehler < 1,7% 
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e) Schlitzblenden- Kondensator (Abb. 8) 
Es wird D>h>R vorausgesetzt. Man erhält 
4eU_ 
wobei 9 =qg- =0 und =U_+ U.-coswmt angenommen wurde, 


& = 


Den Feldverlauf auf der y-Achse kann man wieder annähern durch F = i 


mit q = und auf der x-Achse durch asi mit a7 der 
Fehler ist für |qg{ < 0,5 kleiner als 1%. 


h | 
|7 ° | © 
x r 
h | h | 
t 
Abb. 8. Schlitzblenden-Kondensator Abb. 9. Kreislochblenden-Kondensator 
d) Kreislochblenden- Kondensator (Abb. 9) 
Wie unter c) wird hier 


Eine gute Näherung für den Feldverlauf auf der z-Achse ist F x a fiir 


m 
la| < 0,5 mit q -\5 + für sehr große Werte ist F =<. signum q eine 
gute Näherung. Auf der r-Achse gilt hier genähert F x = an — mit q = ia 
— q . 

Zum Schluß soll darauf hingewiesen werden, daß bei der Verwendung von 
einem symmetrischen Vierpolfeld mit hyperbolischen Elektroden ?1)22)2) bei 
exakter Ausführung keine nichtlinearen Glieder auftreten, was sich in der 
Praxis nur annähernd verwirklichen läßt. 


VI. Periodische Lösungen auf den Symmetrieachsen der speziellen Anordnungen 
Im vorigen Abschnitt wurde gezeigt, daß sich die Feldverläufe auf den 


Achsen der speziellen Anordnungen für |q| S 0,5 durch F, “rie und 
durch F, = i £ Fr annähern lassen. Wir untersuchen daher jezt die periodischen 
4 


2!) W. Paul u. Steinwedel, Z. Naturforschg. 1958, 8a. 
22) W. Paul u. M. Raether, Z. Physik 140, 262 (1955). 
2) W. Paul u.a., Z. Physik 152, 143 (1958). 
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Lösungen der Differentialgleichung 
177 q 
q + (a — - cos 2t)- (18) 


(für @ <1 ergibt sich (12)). 

Die Amplituden der Schwingung B,-sint-+ A,-cost bei k=a =0 
sind in 1. Naherung durch 

2-8 
gegeben. In der Abb. 10 ist die Amplitude A aus q = A-sin (r+yx,) mit 
d, =tgyı füry = 0; 0,1; 0,3 und 0,5 eingezeichnet. Die allgemeinen Aussagen 


des Abschnittes 3 treffen 20 


hier natürlich auch zu (z. B. 


Stabilität). 
= 
4 


Die 2. Näherung wurde 
nr fry=a=k=0 
durchgeführt. Außer der Er- Pr 
höhung der Genauigkeit für 10 7 
B(> 0) nahe 0,908 erklärt 
diese 2. Näherung zusätzlich as LO Ta 
eine experimentell beobach- 


44% 


tete Schwingung mit großer / 

Amplitude zwischen den 0 

parallelen Drähten (auf 2- B * 
Achse, d.h. q,). Für |B]< Abb. 10. Amplitude A für y = 0; 0,1; 0,3 und 0,5 in 
0,908 ist g, = stabil; er- Abhängigkeit von ß bei a = 0 


reicht |3| den Wert 0,9098 

(der exakte Wert ist 0,90805!), wird g; = 0 instabil, q; beginnt sich auf- 
zuschwingen. Dadurch wird || effektiv vergrößert, bis etwa |B|/(1 — |g+ max|*) 
s 7,51 wird (Wert der Grenzkurve), wo wieder eine stabile Schwingung mög- 
lich ist. 


pr = 
04 
ty Naherung > 
Näherung 02; 
Qs . 
01 
/ 
0 l 0t- — — 
09 m 13 15 17 19 27 09 095 0 RB 
Re Abb. 12. Vergleich zwischen gemessener 
Abb. 11. Maximale Amplitude in 1. und 2. Nähe- und berechneter Amplitude für drei 
rung bia=y=k=0 Teilchen 


In Abb. 11 ist die maximale Amplitude von q_ = B,-sint + B,-sin3r 
für 6 > 0 eingetragen, man erkennt die Abweichung der ersten Näherung 
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von der zweiten. Der Wert B=|B,| + |B,| läßt sich sehr gut durch 


(B — 0,908) (19) 
annähern (siehe kleine Kreise in Abb. 11). 

Das Experiment bestätigt innerhalb der Meßgenauigkeit und der gemachten 
Näherungen den Verlauf der Amplitude. In Abb. 12 sind Meßpunkte von drei 
verschiedenen Teilchen und die Kurve nach (19) eingetragen. Die Reibung 
wurde dabei vernachlässigt, was bei diesen Teilchen einen Fehler kleiner als 
3%, bedeutet (k = 0 läßt sich durch geeignetes homogenes elektrisches Gleich- 
feld E, immer erreichen). 

Die Schwingung der Periode x, verursacht durch das Erdfeld, wurde mit 
x=y=0 und || <1 untersucht, Bezeichnung q_—= Ap + Ag - c082r 
bzw. 9, =@,-+ a,-cos2t. Die Abb. 13 und 14 geben A, und a, über ß auf- 

getragen für k = 0,01 und zum Teil k = 0,04 


Naherung an, die Stabilitat der méglichen Schwingung 
fir» / ist gleichfalls angegeben (k > 0). Experimentell 
40.01 / wurde bestätigt, daß für A, < 0 das Teilchen 
3 + / bei Verkleinerung der Wechselspannung immer 
A. stärker „durchhängt‘, bis es bei einem be- 
ii stimmten A, (nach Rechnung » —0,362) aus 
2 7 
/ ~~ lnstail~> 102 
/ -04 u 
/ . 
0 -06 > 
0362 %-004 | stabil frlineare instabil 
pe 0 02 04 06 08 
Abb. 13. A, als Funktion von ß ag 
beiy=a= 0 und k= 0,01 so- Abb. 14. A, bzw. a, als Funktion von ß bei y = « =0 
wie k = 0,04 und k = 0,01 


der Anordnung fällt (Messungen ergaben den Wert —0,35 + 0,02). Die Schwin- 
gungen mit A, > 1,529 wurden ebenfalls beobachtet. Für g; (x- oder 
r-Achse) sind z-periodische Lösungen möglich, falls 0 < |a,| < 1. Das Teil- 
chen fällt nicht aus dem System, sondern auf den Rand der Anordnung, 
falls 8 > 0 geht (da x = 0). Diese Lösung ist stabil. 


VII. Zwei Teilchen gleicher Ladung und Masse in axial-symmetrischen 
Potentialfeldern 
Bisher wurde nur die Bewegung eines einzigen Teilchens verfolgt. Es soll 
nun kurz darauf eingegangen werden, wie sich mehrere (insbesondere 2) 
Teilchen, die zur Vereinfachung der Rechnung alle gleiche Masse und Ladung 
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haben mögen, in der Umgebung von Sattelpunkten axialsymmetrischer Po- 
tentiale verhalten. Das Experiment zeigt, daß die Teilchen für |8| < AY 
in der Ebene z = 0 mit der aufgeprägten Frequenz w gegeneinander schwingen 
(Bezeichnungen wie oben). Zwei Teilchen ordnen sich so an, daB 1, = —t% 
gilt (siehe Abb. 15). Bei n Teilchen beträgt der gegenseitige Winkelabstand 
360°/n. Es können sich bei vielen Teilchen auch einige ,,konzentrische Ringe“ 


yy 03 \ ix 
full Y, \ 10" 
Yj \ —— 


~ 


— 
— — 


IA 


N 
N 


0 | 

0 07 02 Q3 04 05 06 07 0809 10 
Abb. 15. Anordnung zweier Teil- ad 
chen in der Symmetrieebene Abb. 16. a,(ß) bei « = k=0 mit p= 10; 10-4; 
des axialsymmetrischen Systems 10 bei y = 0 und p = 10> bei y? = 0,5 


aus mehreren Teilchen bilden?0). Die Bewegung von zwei gleichen Teilchen 
können wir bei 3=%3=(0 und y; =y3 =0 durch eine eindimensionale 
o 


Differentialgleichung beschreiben. Mit den Bezeichnungen q, = 2 r 


( = V 4 für Kreislochblenden-Kondensator, o = \+ fiir Kreisring) und der 
Feldannäherung q,/(1 — 9?) erhält man die Gleichung (Bezeichnung wie oben) 
p 
mit 


R 


Wir machen den durch das Experiment begründeten Näherungsansatz 


q+ = 4) (1 + d,- cos 2t + dy - sin 27). 


In den Versuchen von H. Straubel ist 10% > p> 10%. Die Abb. 16 zeigt 
einige Kurven a, (8) mit p und y als Parameter. Die Pfeile geben die Amplitude 
an. Mit wachsender Wechselspannung nähern sich die Teilchen dem Sattel- 
punkt r—=z= 0, wobei die Amplitude wächst. Für |f| > Ay treten auch 
sprunghaft 2 -periodische Schwingungen in z-Richtung auf, wie das Ex- 
periment bestätigte. 
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VIII. Meßmöglichkeiten zur Bestimmung von Ladung e und Masse m 
des Teilchens 


In die Bewegungsgleichung eines Teilchens gehen e/m und über die Rei- 
bungskonstante x die Masse m als gewöhnlich zu messende Größen ein, die 
anderen Parameter sind bekannt. 


Für x ~ 0 kann man aus den vorangegangenen Untersuchungen nur Meß. 
möglichkeiten für e/m ableiten. Nur wenn die Reibung berücksichtigt wird, 
welche als Funktion der Masse und der Geschwindigkeit des Teilchens bekannt 
sein muß, ist mit Hilfe dieser Untersuchungen eine Bestimmung von Ladung 
und Masse selbst möglich. Wir können hier x als Funktion von m nur dann 
explizit angeben, wenn wir Stokessche Reibung mit allen Einschränkungen 
der Gültigkeit annehmen“). Bei den von H. Straubel durchgeführten Ex. 


perimenten waren die Parameter so gewählt, daß |B] <2 und k2? <1 erfüllt 
war. Der Parameter y = = wurde unter der Annahme x = 627,1% 


(Stokessche Reibung) für Bewegung in Luft bei Zimmertemperatur und 

w = 100.571 abgeschätzt, 

qn siehe hierzu Abb. 17 (d; = 

0 10 10 10 10 10 10 Dichte des Teilchens, n, = 

Zähigkeit der Luft ~ 1,82 

- 10-4 Poise; ry: Radius des 


tia Teilchens zwischen 3 - 10% 
dr ST und 102 cm). Man sieht, 

(5) 4 daß fast für alle Teilchen 
10 y<1 gilt, meist sogar 

+ 0,8... 1,0 gr} 
0" ie em’). Somit sind die in den 
vorangegangenen Abschnit- 

tenabgeleiteten Näherungs- 


brauchbar. 
Abb. 17. Reibungskonstante y in u vom Es sollen nun kurz die 
Teilchenradius r, bei 50 H sich aus den obigen Unter- 


suchungen ergebenden MeB- 
möglichkeiten skizziert werden. Diese Meßmethoden können natürlich noch 
aus anderen Zusammenhängen abgeleitete hinzugefügt werden (z. B. direkte 
Wägung oder Ausmessung des Durchmessers bei festen Teilchen, Streu- 
lichtmessungen, langsamer Fall der Teilchen im Erdfeld u. a.). 


a) „Sprungpunktmessung“ 


Das Einsetzen der subharmonischen Schwingung (Periode 22) auf der 
z-Achse axialsymmetrischer Anordnungen liefert eine sehr empfindliche und 
genaue Meßmöglichkeit (siehe Abschnitte 3 und 5). Zweckmäßig hebt man durch 
ein homogenes elektrisches Gleichfeld die Wirkung des Erdfeldes auf (k, = 0)), 
so daß sich das Teilchen für |8| < AY im Sattelpunkt in Ruhe befindet. Die 


*4) W. Gerlach, Handbuch d. Physik XXII/1, 1933. 
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Annahme Stokesscher Reibung ist dann gerechtfertigt. In brauchbarer 
Näherung gilt 
By 0,8245 + 47%. 


Hieraus ergibt sich wegen ß ~ e UL!m w* und y ~ 1/o - m*/s der Zusammen- 
hang 


g e- U2 


worin S, und S, bekannte Konstanten sind (ohne Reibung S, ~ 0). Diese 
Messung liefert also einen Zusammenhang zwischen m und e/m. Für relativ 
große Teilchen ist y?< 1 und man erhält mit S, ~ 0 eine gute Annäherung 
für e/m. Durch eine weitere unabhängige Messung kann man dann e und m 
bestimmen. 


b) Amplituden- und Phasenmessungen 


Die Messung der Amplitude der subharmonischen Schwingung liefert einen 
weiteren Zusammenhang zwischen e/m und m, der den entsprechenden Ab- 
schnitten zu entnehmen ist. Diese Methode ist aber weniger empfindlich und 
ungenauer durch die Näherungsansätze. Natürlich können auch die anderen 
abgeleiteten Amplitudenbeziehungen zu Messungen benutzt werden. Des- 
gleichen können Messungen der Phasenverschiebung x (siehe oben) durch- 
geführt werden; in Näherung gilt oft (y?< 1, |ß| # 6%) tg y = y, so daß hier 
y (und damit die Masse) direkt abgeschätzt werden kann. 


e) Statische Messung von e/m 


Man hebt unter der Bedingung |B| < BY die Wirkung des Erdfeldes (g, 
etwa in z-Richtung) durch ein homogenes elektr. Gleichfeld auf (E,). Dann 
gilt wegen k, = 0 die Beziehung e/m —=g/E,. Für |B| nahe py ist diese Methode 
sehr unempfindlich, deshalb muß man soweit als möglich fp (d.h. UL) ver- 
ringern. Diese Methode hat den Vorteil, daß man mit dem Wechselfeld das 
Teilchen noch ,,festhalten kann, solange nicht exakt k, = 0 ist. Die Durch- 
führung von Messungen bestätigte insbesondere die Vorteile der vom Ver- 
fasser vorgeschlagenen Kreislochblenden-Anordnung für diese Messungen. 
Zur Messung von e und m empfiehlt sich also eine Kombination der ,,Sprung- 
punkt-“ und der statischen Methode, die am empfindlichsten sind und kurz 
nacheinander angewendet werden können. 


d) Relativmessungen 


Oft interessiert man sich für langsame zeitliche Veränderungen von e, m 
oder e/m (z. B. Verdampfen u.a.). Durch Sprungpunktmessungen zu den 
Zeiten t, und t, ergibt sich die Relation 


le} ___ |U| \ —flel wi 
V0,8245 + 472), \m /0,8245 + 472 


fir y? < 1 gilt 


(lel), 


173 


SE 
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e) Bei annähernd gleichen Teilchen kann man e?/m bestimmen durch 
Messung des gegenseitigen Abstandes 2 A und der Amplitude C (siehe hierzu 
Abschnitt VII), es gilt 


= 
m 


Ist e/m bereits bekannt, genügt bei y? << 1 die Messung von A; man berechnet 
dann e?/m (und damit e bzw. m) aus 


Zum Schluß sei noch darauf hingewiesen, daß sich die Ladung des Teilchens 
auch bei konstanter Masse im Wechselfeld zeitlich merklich ändert. 


Die vorliegende Arbeit wurde im Theoretisch-Physikalischen Institut der 
Universität Jena angefertigt. 

Herrn Prof. Dr. K. Schuster möchte ich an dieser Stelle meinen Dank 
für die Förderung dieser Arbeit aussprechen. Für wertvolle Anregungen und 
Unterstützung von der experimentellen Seite her bin ich Herrn Prof. Dr. 
H. Straubel zu Dank verpflichtet. 


Berlin-Pankow, Wissenschaftlich-Technisches Büro für Reaktorbau. 


Bei der Redaktion eingegangen am 18. Dezember 1959. 
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Untersuchungen zum verallgemeinerten 
Albedo-Problem’) 


Von M. Richter 
Mit 5 Abbildungen 


Inhaltsiibersicht 


Mit Hilfe einer 500 mC Ra—Be-Neutronenquelle in einer Paraffinanord- 
nung wurden Untersuchungen durchgeführt, die Aufschluß über die Genauig- 
keit theoretischer Betrachtungen zum verallgemeinerten Albedo-Problem 
bringen sollten. Die Ausmessung des Neutronenfeldes erfolgte mit Dysprosium- 
Sonden. Aus den Meßergebnissen ergibt sich die Brauchbarkeit der Theorie, 
deren Anwendung bei anderen Diffusionsproblemen demnach durchaus 
gerechtfertigt ist. 


Problemstellung 


Der Begriff der Albedo wurde schon von Fermi in die Neutronenphysik 
eingeführt. Man stelle sich einen ausgedehnten Paraffinblock vor, in dem 
Quellen thermischer Neutronen isotrop über den gesamten Raum verteilt sind. 
Die Neutronendichte an einem Punkte mitten im Paraffin möge mit o, bezeich- 
net werden. Für den Fall, daß sich unmittelbar daneben ein dickes Kadmium- 
blech befindet, sei die Dichte 04. Das Kadmium soll so dick sein, daß alle 
darauf fallenden Neutronen absorbiert werden. Seine Kantenlänge soll 
mehrere Diffusionslängen betragen. Als Albedo bezeichnet man das Ver- 


hältnis . 
Die Dicke y des Kadmium-Absorbers kann nun zwischen 0 < y < co 
variiert und die Neutronendichte o, in Abhängigkeit davon gemessen werden. 
In diesem so verallgemeinerten Albedo-Problem interessiert das Verhältnis > > 
Für das Albedo-Problem ergibt sich aus der linearen Fermischen Diffu- 
sionstheorie?) 


=14+YN, (1) 
wobei N das Verhältnis des Streuquerschnittes o, zum Absorbtionsquerschnitt 
6, ist. Bei Berücksichtigung der Anisotropiekorrektur aus der Transport- 
theorie wird aus (1) 


— 14 (1-c0sB). (2) 
oo 
Hierbei ist cos # der Mittelwert des Kosinus aller möglichen Streuwinkel. 


!) Aus der Dissertation, eingereicht in Rostock 1959. 
*) E. Fermi, Ric. Sci., VII—II, 13 (1936). 


asi 


222 Annalen der Physik. 7. Folge. Band 6. 1960 


Auch für das verallgemeinerte Albedo-Problem läßt sich eine Lösung 
finden, die allgemein 


— 1 +VN (1 — cos 8) - f(y) 
lautet. Aus der linearen Theorie ergibt sich 
3 

f(y) =Tg (4) 


Es ist leicht zu sehen, daß fiir yoo f(y)=1 wird, Gl. (3) also in (2) über- 
geht. 

Aus prinzipiellen Erwägungen heraus nahm nun Lyons?) an, daß (4) im 
Gebiete 0 < y < oo ungenau wird, da bei der Summierung im linearen Modell 
nicht alle tatsächlich vorhandenen Neutronengruppen erfaßt werden. Er 
brachte (4) in eine andere Form und erhielt schließlich 


Darin ist [A,] eine Abkürzung für den Ausdruck 
[Ay] = 0,706 [(hy + 2 hy + hg) — 0,217 (hg + hy — hy — hg) (6) 


— 0,115 (hg — 2 hy + hy) — 0,065 (hg — hy — 2 hig + 2 + — hg) 
Die Abkürzungen h, bedeuten 


1 oo 
d 
h(y) = f que *=ev—y 
0 y 


fir n > 0 


1 


h, (y) = n+1 — (y)] = (Y)- (7b) 


Bekanntlich ist das n-te Exponential-Integral fiir positive Werte des Argu- 
mentes y definiert. In Abb. 5 sind die Ergebnisse sowohl aus der linearen 


V3 


Fermischen Theorie (te? v) als auch aus der exakten Theorie graphisch 


dargestellt. 
Die Größe f (y), die auf der Ordinate aufgetragen ist, ergibt sich aus (2) 
und (3) aber auch zu 


f(y) = (8) 
fo _4 
P00 


und kann also aus experimentellen Werten ermittelt werden. Merkliche Ab- 
weichungen der Ergebnisse aus beiden Theorien lieBen eine experimentelle 
Nachpriifung als sinnvoll erscheinen. 

Um die Frage zu beantworten, bei welchen Absorberdicken y die Messungen 
am ehesten eine Entscheidung zulassen, wurden die Differenzen der aus det 


3) D. Lyons, Ann. Physik (6), 4, 379 (1949). 


4) Dieser Ausdruck gilt fiir o = - 1: >0,. 


0, +64 


sie al 
In d 
renze 
gehö 
einge 
Inteı 
E 
Expe 
sung 
L stant 
von 
einig 
etwa 
I 
eine 
der 
zuer: 
mit 
fiillt 
von 
daß 
| 
kom 
werd 
Stöß 
ordn 
I 
such 
war. 
dadu 
Krei 
läßt. 
über 
Der 
alen 
s 


Osung 


1e Ab- 
entelle 


ungen 
us der 


M. Richter: Untersuchungen zum verallgemeinerten Albedo- Problem 223 


linearen und der exakten Theorie errechneten Werte gebildet. In Abb. 1 sind 
sie aufgetragen. Die Kurve 1 enthält 


f (Wex 
f (YWex 


In der Kurve 2 sind die Diffe- 
renzen der zu jedem f(y)-Wert 
gehörenden Neutronendichten 
eingetragen. Danach ist das 
Intervall 0,5 <y<2 für die 
Messungen am geeignetsten. 


- 100%. 


Say 


Eine Voraussetzung bei den ? J 
Experimenten ist, daß die Mes- Abb. 1. Relative Differenzen, die sich aus der 


Experimente | 
¥ 
linearen und der exakten Theorie ergeben 


sungen in einem Gebiet kon- 
stanter Dichte der thermischen Neutronen durchgeführt werden. An Hand 
von Abschätzungen läßt sich zeigen, daß das Feld eine Ausdehnung von 
einigen Diffusionslängen haben muß. In Wasser oder Paraffin sind das 
etwa 90—100 mm. 

Die Herstellung eines räumlich konstanten Neutronenfeldes läßt sich durch 
eine bestimmte Form der Neutronenquelle erreichen. Man kann das Prinzip 
der Helmholtzschen Stromspulen zugrunde legen. Diese Möglichkeit griff 
zuerst Volz®) auf, der zwei koaxiale Glasrohre 
mit Berylliumflittern und Radium-Emanation 
füllte. Die Form und die Maße dieser Quelle 
sind der Abb. 2 zu entnehmen. Bei dieser Art 
von Quellenanordnung sind die Radien der 
Ringe und deren Abstände so zu bestimmen, 
daß zunächst eine konstante Dichte von Neu- 
tronen mit der Energie von 1 MeV zustande 
kommt. In einem geeigneten Bremsmedium 
werden die Neutronen dann nach wenigen 
Stößen thermisch. Die Dimensionen der An- 
ordnung hängen vom Bremsmedium ab. 

Im vorliegenden Falle stand für die Unter- 
suchungen eine 500 mC Ra—Be- Quelle zur Ver- 
fügung, die in einem Platinzylinder eingelötet er, 
= Nun läßt sich ein Ring der Volzspule uch Volk wr Herstellung 
adurch schaffen, daß man die Quelle auf einer eines räumlich konstanten Neu. 
Kreisbahn mit gleichem Durchmesser umlaufen tronenfeldes 
läßt. Eine Summierung der Neutronendichte 
über einen Umlauf ergibt eine konstante Quellstärke auf dem gesamten Ringe. 
Der andere Ring wird erhalten, wenn die Quelle auf einer zur ersten koaxi- 
alen Kreisbahn in einer zweiten Ebene rotiert, wobei der Abstand der Ebenen 
dem der beiden Volzringe entspricht. Summiert man die Neutronendichten 


5) H. Volz, Z. Physik 121, 201 (1943). 
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jedes Punktes im Innern der Spule über je einen Umlauf auf beiden Ringen, 
so zeigt sich, daß eine annähernd konstante Dichte erhalten wird, wenn die 


t, < T), (Dy) 


erfüllt ist, wobei t, die Zeit für einen Umlauf und 7), (Dy) die Halbwertszeit 
des Dysprosium sind. Die Meßergebnisse dieser Pseudo-Spule können denen 


Bedingung 


Paraffinglock orenb. Paraffinzytinder 


| Absorber 


Abb.3. Versuchsapparatur mit rotierendem 
Paraffinzylinder 


der echten Volzspule gleichgesetzt 
werden. 

Die Abb. 3 zeigt eine schematische 
Darstellungder benutzten Anordnung, 
Die Sonde und der Absorber sind 
zwischen den beiden Ringebenen an- 
geordnet. Es ist auf Grund des Auf. 
baues notwendig, die Aktivierungen 
und Ausmessungen zu teilen und für 
jeden Ring einzeln durchzuführen. 
Das Ergebnis, das zur weiteren Aus- 
wertung herangezogen wird, besteht 
aus der Summe dieser beiden Teil- 
ergebnisse. 

Zu beachten ist dabei, daß 

1. Der Absorber und die Sonde 
sich stets an genau dem gleichen Ort 
befinden. 

2. Die beiden Ringebenen ein- 
wandfrei fixiert und stets reproduzier- 
bar sind. 

3. Die Umlaufgeschwindigkeit kon- 
stant ist. 


Im vorliegenden Falle sah die technische Ausführung noch etwas anders 
aus, wobei das Prinzip aber nicht verändert wurde. Die Quelle blieb fest und 


fanden, an ihr vorbei. 


Ro-Be-Quelle rotierte nicht, vielmehr 
drehte sich der Zylinder, 
in dem sich die Sonde 

und der Absorber be- 


Anordnung I 


Aus Abb. 4 ist die ge- 
naue Lage von Quelle, 
Sonde und Absorber für 
beide Fälle zu ent- 
nehmen. Der Paraffin- 
zylinder war auf einem 


Metallteller gelagert, der 
gelag 
Anordnung I durch einen stark unter- 
setzten Elektromotor 
Ra-Be-Quelle 
Abb. 4. Schematische Darstellung der beiden in der Ver- 5 me 
suchsapparatur benutzten Anordnungen von Neutronen- mlauigeschwindigkel 
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gen pro Minute. Der rotierende Zylinder paßte genau in die Ausbohrung 
eines größeren Paraffinklotzes, dessen Höhe 500 mm und dessen Kantenlänge 
250 x 250 mm war. Nahe der Zylinderwand befand sich in einer passenden 
Bohrung die Ra—Be- Quelle. 

Über die Reinheit des verwendeten technischen Paraffins kann nur so viel 
gesagt werden, daß sich vor allem kein Bor und kein Kadmium in chemisch 
nachweisbaren Mengen darin befanden. 

Als Sondenmaterial wurde Dysprosium benutzt, das in Form von Dy,O, 
vorlag. Es hat einen besonders großen Einfangquerschnitt für thermische 
Neutronen und einen verhältnismäßig kleinen für schnelle Neutronen. Die 
Halbwertszeit des in Frage kommenden Isotops Dy!# beträgt 140 Min. Die 
Sonde wurde aus 15 mg Dy,0, hergestellt, das mittels Zaponlack auf eine Glim- 
merfolie geklebt wurde. 


Meßergebnisse 


Die Abb. 1 zeigt, daß die größten relativen Differenzen der beiden zur 
Diskussion stehenden Theorien zwischen 0,5 < y < 2 auftreten, die Kadmium- 
dicken von etwa 50 u bis 150 w entsprechen. Solche Folien wurden benutzt. 
Ihre Dicke wurde mit einem Orthotest bestimmt, das auf + 1 genau miBt. 
Der Meßfehler war also, abgesehen vielleicht von der dünnsten Folie, kleiner 
als der statistische Fehler bei der Aktivitätsbestimmung. 

Die Aktivitäten wurden auf Sättigungswerte umgerechnet. Diese Um- 
rechnung war wegen der unterschiedlichen Bestrahlungszeiten erforderlich. 
In einzelnen Fällen, wenn das Standardpräparat vor und nach der Aktivitäts- 
messung um mehr als ein Prozent vom üblichen Wert abwich, wurde am MeB- 
ergebnis eine Korrektur vorgenommen. 

Die Korrekturen hinsichtlich des Anteiles der Sondenstörung und der 
schnellen Neutronen wurde berücksichtigt und in der gleichen Größe von allen 
Meßwerten abgezogen, da die 
schnellen Neutronen in keiner 


oder jedenfalls in einer hier fü 

nicht ins Gewicht fallenden Li 

Weise von der Cd-Dicke be- i 

einflußt wurden. | 1 
Der Fehler setzt sich aus 06 Y% 


den statistischen Fehlern der 
ax 


Werte zusammen, die am Zu- 
standekommen von f(y) be- 
teiligt sind. Bei jeder Messung 
wurden mindestens 5000 Im- 
pulse registriert, so daß der 
statistische Fehler kleiner als 


1. Aus dem linearen Fermischen Modell 


2. Aus einer exakten Theorie 
(nach lyons) 


Z 3 ey 
1,6 ist. 5 
Prozent ist. In der Abb. 5 Abb. 5. Meßergebnisse einschließlich Fehlerbreiten, 


sind die Meßergebnisse mit eingetragen in die Kurven aus den zur Diskussion 
ihren Fehlerbreiten einge- stehenden Theorien 


tragen. 

Die Meßergebnisse bestätigen recht eindeutig die von D. Lyons benutzte 
exaktere Theorie. Da die Unterschiede zwischen den beiden Theorien maximal 
(bei y = 1) 10—11 Prozent betragen, die MeBfehler aber nur + 4 Prozent aus- 
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machen und die Meßpunkte zudem noch eindeutig in der Nachbarschaft de 


Kurve aus der exakten Theorie liegen, ist diese Behauptung ohne weiters 
gerechtfertigt. 


Für die Unterstützung bei dieser Arbeit sei auch an dieser Stelle Herm 
Professor Dr. G. Richter und Herrn Dr. D. Lyons gedankt. 


Zeuthen bei Berlin, Deutsche Akademie der Wissenschaften zu Berlin, 
Kernphysikalisches Institut. 


Bei der Redaktion eingegangen am 22. Dezember 1959. 


] 

stoff 
halt 

bem 

west 
übeı 

Aus: 

] 

Seit 

flac! 

Que 

gleic 

zu ¢ 

Med 

ist 1 

unt 
unn 

tror 

troy 

D.] 

schı 

wol 

der 


haft der 
weiteres 


> Herm 


Berlin, 


Der Dichteverlauf thermischer Neutronen 
in der Nähe des Randes eines Streumediums'‘) 


Von M. Richter 
Mit 6 Abbildungen 


Inhaltsübersicht 


Es wurde der Dichteverlauf thermischer Neutronen am Rande eines wasser- 
stoffhaltigen Moderators gemessen. Um die notwendige Genauigkeit zu er- 
halten, wurde der Reaktor in Rossendorf bei Dresden als Neutronenquelle 
benutzt. Das Ergebnis bestätigt, daß die exakte Theorie bei diesem Problem 
wesentlich genauer ist als die lineare Fermische Diffusionsvorstellung. Dar- 
über hinaus kann aus den Messungen mit einer gewissen Wahrscheinlichkeit eine 
Aussage über den Grad der Anisotropie gemacht werden. 


Problemstellung 


Es wird ein wasserstoffhaltiges Streumedium angenommen, das nach einer 
Seite unbegrenzt ist. Beiz — 0 hat es eine unendlich ausgedehnte Begrenzungs- 
fläche und in einiger Entfernung davon im Streumedium eine ausgedehnte 
Quelle thermischer Neutronen. Diese soll so beschaffen sein, daß die Flächen 
gleicher Neutronendichte in der Nähe der Grenzfläche eben sind und parallel 
zu dieser liegen. Von der anderen Seite her treten keine Neutronen in das 
Medium ein. Ferner wird vorausgesetzt, daß die Neutronenabsorption gering 
ist und vernachlässigt werden kann. 

In einer Reihe von theoretischen Überlegungen zu diesem Problem hatten 
unter anderen schon Mark?), Placzek*) und Seidel?) festgestellt, daß in 
unmittelbarer Randnähe das lineare Fermische Modell versagt und die Neu- 
tronendichte sicher nicht richtig wiedergibt. Unter der Voraussetzung iso- 
troper Streuung entwickelten sie verbesserte Theorien. Schließlich konnte 
D.Lyons®) die Dichte auch für den Fall einer anisotropen Streuung be- 
schreiben, wobei der Grad der Anisotropie sogar noch offen blieb. 

Es existiert für die Randdichte yw, (z) die Beziehung 


Yo (2) = 3 (2 + 20) — 2 Ye (2), (1) 
wobei 2, die Extrapolationslänge und z die Entfernung vom Rand in Einheiten 
der freien Weglänge ist. Die Größe 2 y, (z) ergibt sich anschaulich als die 


!) Auszug aus der Dissertation, eingereicht in Rostsock 1959. 
2) C. Mark, Physic. Rev. 72, 558 (1947). 

3) G. Placzek, Physic. Rev. 72, 556 (1947). 

4) G. Placzek u. W. Seidel, Physic. Rev. 72, 550 (1947). 

5) D. Lyons, Ann. Physik (7), 1, 400 (1958). 
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Differenz zwischen der Asymptote am Dichteverlauf im Innern des Streu. 
mediums und der wahren Dichte an der Stelle z. Dieser Wert 2 y, (z), in dem 
ein Parameter x noch unbestimmt und 
experimentell zu ermitteln ist, wurde 
Aus dem von Lyons berechnet. In der Abb. 1 
linearen Mode sind die Asymptote und die theo. 


Extrapotation aus ö Ausdem retisch möglichen Dichteverläufe dar. 
der exokten Model gestellt, in Abb. 2 die Abhängigkeit 
A der Größe 2, (z) sowohl von z ak 

a auch vom Parameter x. 
Da die Extrapolationslänge und 
(die Dichte der thermischen Neutronen 
I A 0 in Abhängigkeit von z experimentell 
Om, V3 =— untersucht und festgelegt werden 
Abb. 1. Dichteverlauf und Asymptote aus können, ist also 2y,(z) über die 
den Theorien Gl. (1) zu ermitteln. So kann einmal 


der von den exakten Rechnungen 
vorausgesagte Dichteverlauf nachgeprüft werden, zum anderen sind bei 
genügender Meßgenauigkeit Aussagen über den Parameter x zu erwarten. 


Experimente 

Der zu messende Effekt ist sehr klein, weswegen der Reaktor in Rossendorf 
bei Dresden als Neutronenquelle benutzt wurde. Nach einigen orientierenden 
Vorversuchen zeigte sich, daß eine thermische Leistung des Reaktors von 
300 bis 500kW für die 
Aktivierungen genügte. 
Dabei blieben die erforder- 
lichen biologischen Strah- 

lenschutzmaßnahmen 
noch in _erträglichen 
Grenzen. 

Wegen des hohen Neu- 
tronenstromes war 68 
nicht erforderlich, die Dy- 
Sonden bis zur Sättigung 
zu aktivieren. Im allge- 


a 02 —: 03 04 meinen reichte eine Halb- 
Abb. 2. Abhängigkeit der Größe 2 y, (z) von z und vom wertszeit aus, also etwa 
Parameter x 2,5 Stunden. Dieser Um- 


stand wiederum machtees 

erforderlich, mehrere Sonden zur Verfügung zu halten, denn eine einmal ak- 

tivierte Sonde konnte frühestens nach 20 Stunden ein zweites Mal benutzt 

werden. So wurden die sechs Sonden Nr. 11 bis 16 der nachfolgenden Tab.l 

aus Dy,O, angefertigt, das mittels Zaponlack auf 60 u starke Aluminiumfolien 

aufgespritzt wurde. Die Dysprosiumsonden wurden so dünn hergestellt, daß 

die durch sie zu erwartende Feldstörung wesentlich kleiner als ein Prozent 
blieb. 

Um die Messungen mit den verschiedenen Sonden miteinander vergleichen 

zu können, mußten ihre relativen Größen bekannt sein. Sie wurden durch 
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mehrere Aktivierungen in unmittelbarer Nähe einer Ra—Be-Quelle und 
anschließende Ausmessung der induzierten Aktivitäten ermittelt. Die MeB- 
ergebnisse sind in Spalte 2 der Tabelle eingetragen. In Spalte 3 ist ein Faktor 
F, aufgeführt, der der reziproken Größe der 


jeweiligen Sonde proportional ist. Durch Tabelle 1 
Multiplikation eines Meßwertes mit dem Nr. _Impulse/min | F, 
Faktor F, konnte dieser normiert und zu 
weiteren Vergleichen herangezogen werden. 11 548 + 0,7% | 0,597 
2 | 9 

Die Versuchsanordnung aus einem festen - | = + 10% ae 
Moderator ist in Abb. 3 dargestellt. Sie be- 14 | 332 40,9% 0.985 
steht aus einem Paraffinklotz mit einer 15 327 + 0,9% 1,000 
zylindrischen Bohrung von 100 mm Durch- 16 | 327+0,9% 1,000 


messer. Senkrecht zur Bohrung wird der 
Klotz auf einer Seite von einer Kadmiumwand abgeschlossen, auf der an- 
deren befindet sich die Eintrittsebene für die Neutronen. Das Reaktorbündel 
verläuft in der Verlängerung der im Klotz angebrachten Bohrung. In die 
Bohrung passend wird eine Anzahl von Zylindern und Scheiben aus Paraffin 
oder aus Polyäthylen eingeschoben, deren Stärke genau bekannt ist. Die 
Dieken der Scheiben lagen bei den Ver- 

suchen zwischen einem halben und drei Paraffinklofz 

Millimetern. Die Sonden waren zwischen 
den Scheiben angebracht. 


Es muß an dieser Stelle bemerkt 
werden, daß es mit dieser Anordnung | __---- 
möglich war, gleichzeitig mehrere Sonden it 
in verschiedenen Abständen von der Rück- 
wand zu bestrahlen. Der Vorteil ist ein- 
leuchtend, wenn man folgendes bedenkt: 


1. Wie schon gesagt, wurde nicht bis 
zur Sättigung aktiviert, sondern immer 
nur eine wesentlich kürzere Zeit. Da aber, ‚Auswechselbarer Zylinder Kadmiumwand 
wenn nur eine Sonde aktiviert wird, in 
der anschließenden Auswertung alle Er- Sonde (zwischen zwei Parattinscheiben) 
gebnisse einheitlich auf Sättigungsaktivi- Abb. 3. Versuchsanordnung aus 
täten umgerechnet werden, wäre es not- Paraffin und Polyäthylenscheiben 
wendig, die Zeit t, für die Aktivierung 
genau festzuhalten. Fehler bei der Messung von t, beeinflussen also die Ge- 
nauigkeit der zu ermittelnden Sättigungsaktivität. Dieser Fehler entfällt bei 
gleichzeitiger Aktivierung zweier Sonden, wenn eine davon als Monitor be- 
nutzt wird. 


2. Die Reaktorleistung ist im Verlaufe von Stunden nicht konstant. Das 
zur Zeit benutzte Verfahren zur Anzeige der Leistung gewährleistet höchstens 
eine Genauigkeit von einigen Prozent. Deswegen ist ein Monitor geradezu 
notwendig, d. h. in unserem Falle eine Sonde, die stets in der gleichen Lage 
aktiviert wird und deren Aktivität ein Maß für die integrierte Reaktorleistung 
ist. Eventuelle Abweichungen von einem Normalwert können so in der er- 
forderlichen Weise korrigiert werden. 
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MeBergebnisse 


In vorbereitenden Untersuchungen wurden der Dichteverlauf der thermi. 
schen Neutronen in Abhängigkeit von der Schichtdicke des Moderators bei 
Abwesenheit von Absorbern, die Eigenaktivität der Sondenträger und die 
Abhängigkeit der Aktivitäten von der Reaktorleistung ermittelt. Der Verlauf 
der Dichte der thermischen Neutronen in Abhängigkeit von der Moderator. 

dicke ist in Abb. 4 dargestellt. Daraus 

Dichte läßt sich folgern, daß es von Vorteil ist, 

die Moderatorschicht bis zur Kadmium. 

wand nicht dicker als 10cm zu machen, 

um noch optimale Intensitäten erwarten 
N zu können. 

N Bei den Auswertungen mußte not- 

N wendigerweise die Transportweglänge 4, 

a bekannt sein. Für Paraffin konnte aus der 

Literatur®) A, = 0,395 + 0,02 em ent- 

0 5 © & 2% 2 m nommen werden. Man kann annehmen, 

Abb. 4. Neutronendichte (in willkür- daß der Wert für Polyäthylen in der 

lichen Einheiten) in Abhängigkeit von gleichen Größenordnung liegt, da die 

der Moderatordicke Wasserstoffkonzentration pro Kubikzenti- 

meter der Substanzen etwa gleich ist. 

Allerdings kann die Dichte des Polyäthylens je nach der Kristallstruktur um 

10%, schwanken. Demzufolge werden auch die Diffusionskenngrößen nicht 

konstant sein. 

Zur experimentellen Ermittlung der Transportweglänge A, des verwendeten 
Polyäthylens wurde die Extrapolationslänge d festgelegt. Sie beträgt bekannt- 
lich 0,71 2,. Die Extrapolationslänge ergab sich aus der Asymptote am Dichte- 
verlauf, die für die weiteren Untersuchungen ohnehin notwendig war. Es 
wurden jeweils zwei Sonden gleichzeitig aktiviert, eine Sonde mit unterschied- 
lichem Abstand von der Kadmiumwand bei z = 0, die andere in der stets 
gleichen Entfernung von 9,6 mm als Monitor. Dabei konnte der erste Monitor- 
Wert auch als Meßpunkt benutzt werden, denn nur die beiden nächsten 
wurden darauf normiert. Zu bemerken ist noch, daß in einer Entfernung von 
3,15 mm von der Wand die tatsächliche Dichte schon etwa 0,3 Prozent unter- 
halb der Asymptote liegt. Bei der Auswertung wurde dieser Umstand beachtet. 

Mit Hilfe der Methode der kleinsten Quadrate ergab sich für die Gerade 
die Gleichung 


Y, = 5,46 + 0,58 (x, — 6,39). 
Dabei verstehen sich die Abszissengrößen x, in Millimeter des Polyäthylens. 
Daraus war 
d = 3,02 mm 
und durch die Beziehung d = 0,71 A, 


A, = 4,25 mm 
zu errechnen. 


®) C.W. Tittle, Physic. Rev. 80, 756 (1950). 
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Die Meßpunkte und die sich daraus ergebende Asymptote sind in Abb. 5 
graphisch dargestellt. Bei den folgenden Messungen des Dichteverlaufes 
wird die Abszisse z in Einheiten der freien Weglänge angegeben, um Vergleiche 
mit der Theorie und eventuellen 

Dichte 
späteren Messungen zu ermög- 
lichen. 

Für den Dichteverlauf ist 
das Intervall von 0,1<z< 0,4 
interessant. Die Messungen wur- 
den in Abständen von 0,530 mm 
und 1,055 mm vom Rande durch- 
geführt. Der Monitor befand sich 
jedesmal, wie schon bei den Mes- 
sungen des Asymptote, in 9,6 mm 


Entfernung von der Kadmium- 

wand. Die Monitorwerte und im Abb.5. Asymptote, die durch einige Meßpunkte 
3eichen Verhältnisdie Meßgröß im Inneren des Moderators festgelegt wurde. 
gleichen Verhattnisdle Mebgroven Auf der Abszisse ist die Dicke der Polyäthylen- 
der zugehörigen Sonden wurden scheiben aufgetragen 


normiert. Da die Sondenstörung 
die Intensität nur um etwa 0,3—0,4%, verkleinert, wurde sie bei der Auswer- 
tung nicht berücksichtigt. In Abb. 6 sind die Meßergebnisse eingetragen. Die 
ausgezogene Kurve ist die sich aus der Theorie für x = 5 (isotroper Fall) er- 
gebende. Der Fehler beträgt etwa + 1 Prozent, wobei der statistische Fehler 
eines einzelnen Meßwertes nicht größer als 0,35%, ist, denn es wurden jedes 
Mal mindestens 8—10 - 10% Impulse 
registriert. Ein Meßpunkt bei z = 0, 
der in der Abbildung eingetragen ist, 
wurde zu weiteren Aussagen nicht her- 
angezogen, da die Randeinflüsse theo- 
retisch schwer zu erfassen und experi- 
mentell nicht auszuschalten sind. 


An dieser Stelle soll bemerkt wer- 
den, daß die Bedingung o,=0 für 
den Moderator (Paraffin) erfüllt sein 
dürfte. Da nämlich die vier Meßpunkte, 
die die Asymptote definieren, innerhalb 0 i 
der Meßfehler eine Gerade bilden, ist u, 0 02 03 
die Absorption in Paraffin zu vernach- Abb.6. Experimentell ermittelter Dichte- 
lässigen. Bei merklicher Absorption verlauf in unmittelbarer Nähe des Ab- 
hätten die Meßpunkte auf einer Expo- sorbers (Abszisse in z-Einheiten) 
nentialkurve liegen müssen. 

Mit Hilfe der Beziehung (1) konnte jetzt 2 wy, (2) errechnet werden. Die 
sich ergebenden Werte sind in Abb. 2 in die theoretische Kurvenschar einge- 
tragen. 


Es dürfte stets mit Schwierigkeiten verbunden sein, den genauen Grad der 
Anisotropie aus derartigen Experimenten zu bestimmen. Die zu erwartenden 
Dichteunterschiede sind so klein, daß die Meßgenauigkeit unter 0,5 Prozent 
liegen müßte, sollten genaue Aussagen über ein gültiges x gemacht werden. 
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Nach den vorliegenden Messungen, der Fehler liegt bei + 1 Prozent, sche 
es möglich zu sein, daß für x ein Wert von 6 oder 7 richtig ist. 


Für nützliche Aussprachen und vielerlei Unterstützung sei auch an dien 
Stelle Herrn Professor Dr. G. Richter gedankt, außerdem aber auch Herm 
Professor Dr. H. Barwich und Herrn Dr. habil. K. F. Alexander vom 
Zentralinstitut in Rossendorf. 


Zeuthen bei Berlin, Deutsche Akademie der Wissenschaften zu Berlin, 
Kernphysikalisches Institut. 
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